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A line of sargassum weed running parallel to the inshore edge of the Gulf Stream at 


40° N., 63° W. There is usually an abrupt change of velocity of half a knot or so on 
both sides of a weed line such as this. Often there are a number of such lines, and 
corresponding velocity steps, along the inshore edge of the Stream. The vessel is about 
300 feet long. This aerial photograph was given to me by Commander William Kielhorn, 
United States Coast Guard Reserve. 
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PREFACE 


The purpose of this book is to describe and explain what is known about 
the Gulf Stream in a way which will interest physical scientists. The name 
“Gulf Stream’ is familiar enough to everyone, but few scientists have any 
knowledge of the nature of this grand natural phenomenon. I hope that 
by means of this book I shall be able to communicate the facts and theories 
concerning the Gulf Stream to a wide scientific audience. Those interested 
in an authoritative treatise which covers the entire field of oceanography 
will do best to refer to The Oceans: Their Physics, Chemistry and General 
Biology, by H. U. Sverdrup, Martin Johnson, and Richard Fleming; those 
who seek an appreciation of the theoretical framework of oceanography 
should read Dynamical Oceanography, by Joseph Proudman. So far as the 
restricted subject of the Gulf Stream itself is concerned, I believe the dis- 
cussion in this work is more comprehensive than that in any other source. 
At any rate, I have sought to make it so. 

I am very much indebted to all my colleagues at the Woods Hole 
Oceanographic Institution, many of whom have spent weary months at 
sea in our uncomfortable little exploring vessels gathering information 
about the western North Atlantic. In particular, I have profited from 
discussion of a descriptive nature with Messrs F. C. Fuglister, C. O’D. 
Iselin, W. S. von Arx, A. H. Woodcock, and L. V. Worthington. My 
collaboration with Mr Donald Parson, Jr., will always evoke memories 
of happy instrument-making. I have had the pleasure of theoretical 
discussions with Dr Willem Malkus, Dr Jule G. Charney, Dr George 
Veronis, and Dr George W. Morgan. Some of the ideas expressed in 
this book arose and took shape in these informal discussions. I do 
not claim them as my own. Wherever recollection permits I have 
indicated their source. And I want to express my gratitude to Dr R. 8S. 
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Arthur, of the Scripps Institution of Oceanography, for his advice and 
encouragement. 

In a field as small as dynamical oceanography it is inevitable that many 
exchanges of information take place by word of mouth. It is always a 
pleasure to discuss problems freely with investigators from other countries, 
and I have been able to talk with most of them. Owing to world conditions, 
I have never met the leading Russian dynamical oceanographer, Professor 
W. B. Stockmann, and hence my acquaintance with his ideas is limited to 
occasional translations of his papers. The Japanese studies of the Kuroshio 
are, for the most part, not treated in this book, partly on account of 
language difficulties, and partly because they are voluminous and deserve 
more careful study than I can give them. The Kuroshio is not unlike the 
Gulf Stream in many ways, and it is hoped that some day a thorough 
comparison can be made. 

Over the course of the years my researches have been generously sup- 
ported by the Office of Naval Research through contracts with the Woods 
Hole Oceanographic Institution. Without this sustenance I could certainly 
not have undertaken such experimental and observational studies as I 
have, and I should not have been afforded the opportunity of associating 
with the many keen minds and eager students of the sea who work at 
Woods Hole. But the conception of and responsibility for a book of this 
kind is, in the last analysis, an individual matter; therefore, it is also 
befitting to record here the fact that the writing of this book has been 
entirely a private undertaking at home and that preparation of the 
manuscript and of the original figures has been at my own expense. 


HENRY STOMMEL 
Woods Hole, Massachusetts 
June, 1955 
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Chapter One 
HISTORICAL INTRODUCTION 


From the time of the recorded discovery of the Gulf Stream to the present 
there have been many ideas about its cause. The historical references! 
I shall make will be restricted to those for which some documentary 
evidence exists—a plan that excludes a large body of speculation about 
early Norse, Arabian, and Portuguese navigators. 


EARLY IDEAS AND EXPLORATIONS 


Although the island of Cuba was first cireumnavigated in 1508, it was not 
until 1513 that the Gulf Stream (specifically, the Florida Current) was 
described by Ponce de Leén, who, sailing from Puerto Rico, crossed the 
stream north of Cape Canaveral and then sailed south to Tortugas. The 
current was so swift that his three ships were frequently unable to stem it 
(see Herrera y Tordesillas, 1601). 

By 1515 Peter Martyr of Anghiera reported various conjectures about 
the Gulf Stream (see the 1577 translation of his Decades). His arguments 
were based essentially upon the principle of the conservation of mass and 
upon the tacit assumption that the current velocity is independent of 
depth. Peter Martyr argued that the North Equatorial Current must 
either (i) pile up large masses of water at the Brazilian coast, or (ii) pass 
through some great straits or passages into the Pacific and thence round 
again into the Atlantic, or (iii) be deflected by the mainland so as to flow 
back into the ocean (Kohl, 1868). The first possibility was ruled out, he 

1 The material in this chapter is drawn from an article by the author in the 


Scientific Monthly for April, 1950. 
I SGS 
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claimed, because the explorers of the Brazilian coast had never noticed 
such a piling up of water. The second possibility was doubtful, because it 
was the consensus among the Spanish navigators that the mainland was 
not open, but presented a continuous barrier to the westward flow. Hence, 
by elimination, only the third possibility was left; and as an example of 
a deflection of the Equatorial Current by the mainland, the Gulf Stream, 
said Peter Martyr, was a case in point. 

The westward flow of the Equatorial Current itself was usually attri- 
buted to the primum mobile—in some manner not clearly understood the 
general westward motion of the celestial bodies across the sky drew the 
water and air of the equatorial regions along with it. 

By 1519 the Gulf Stream was so well known that Spanish ships bound 
for America came by way of the Equatorial Current but, on their return, 
passed through the Florida Straits, followed the Gulf Stream to about the 
latitude of Cape Hatteras, and then sailed eastward to Spain. In this way 
they had favorable winds and avoided contrary currents over the whole 
voyage. 

The sixteenth century marked the beginning of a period of intense 
activity in the western North Atlantic Ocean, of exploration of the coasts, 
and of the search for the Northwest Passage. Many of the early cruises in 
and about the Gulf Stream, chronicled in Kohl (1868), need not concern 
us here. Navigators of various nations investigated the geographic extent 
of the Gulf Stream System. Among these we may number Sir Humphrey 
Gilbert, who first suggested using a deep-sea anchor to determine surface 
drift, Martin Frobisher, Ribault, and Laudonniére. Frobisher and John 
Davis made numerous observations of the Labrador Current. 

Toward the end of the century André Thevet (1575) attributed the Gulf 
Stream to the great rivers that flow into the Gulf of Mexico. It was not 
until many years later, when the transports of mass through the mouth 
of the Mississippi and the Straits of Florida were measured and compared, 
that the utter inadequacy of such an explanation was completely revealed. 
The mass flux through the former is only one one-thousandth of that 
through the latter. 

The sharp line of demarcation between the warm- and cold-water masses 
was apparently first recorded by Lescarbot in 1609. His comment, as 
quoted (from 2d ed., 1612, 2:531) by Kohl (1868, p. 68), reads: 


I have found something remarkable upon which a natural 
philosopher should meditate. On the 18th of June, 1606, in 
latitude 45° at a distance of six times twenty leagues east of the 
Newfoundland Banks, we found ourselves in the midst of very warm 
water despite the fact that the air was cold. But on the 21st of 
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June all of a sudden we were in so cold a fog that it seemed like 
January and the sea was extremely cold too. [Translation.] 


In 1590 John White took a trip from Florida to Virginia. In order to 
stay within the Gulf Stream, he reported, one had to stand far out to sea, 
because along the coast there was a countercurrent—‘eddy currents setting 
to the south and southwest’ (Kohl, 1868). This was the first mention of 
countercurrents on the shoreward side of the Gulf Stream. 

The seventeenth century saw the colonization of the Atlantic coast of 
North America, and the Gulf Stream was of course traversed countless 
times at various latitudes. A number of studies of ocean currents were 
published at this time, works of most varied quality. Varenius (in 1671; 
2d ed., 1681) published a very comprehensive description of the surface 
currents then known, and Isaac Vossius (1663) postulated a complete 
circulation of the North Atlantic Ocean, turning in an ocean-wide clock- 
wise motion. The first chart showing the Gulf Stream was Kircher’s, 
published in 1665 (3d ed., 1678); and the next was a current chart by 
Happelius, in 1685 (see Kohl, 1868). In addition to showing certain true 
features, these charts display certain extraordinary phenomena—for 
example, two distinct surface currents which cross over each other; and 
a great whirlpool off the Lofoten Islands, the legendary Maelstrom. The 
reader can find some of these charts reproduced in Pillsbury’s (1891) 
account of his studies. On the whole, these charts were far superior to the 
theories advanced to explain them. Vossius contended that a great moun- 
tain of water was formed each day at the equator by the heat of the sun, 
and that this water mass was carried westward and broke upon the 
American shore, and then flowed along the coasts in the form of currents. 
Kircher, it is true, suggested that the trade winds contributed to the ocean 
circulation, but he also enumerated other, fantastic, causes. Even well- 
informed men like Kepler had curious ideas about the causes of ocean 
currents. Kepler believed that because the water is only loosely attached 
to the earth it could not keep up with the diurnal rotation and hence fell 
behind, the result being the westward drift of the Equatorial Current 
(Kohl, 1868, p. 87). 

These various theories were, at least, honest attempts to explain physical 
phenomena by an as yet poorly developed physics. In addition, there were 
advanced fantastic theories that enjoyed a certain popularity. An example 
is Merula’s statement (Kohl, 1868, p. 63): 


At the North Pole one finds four large islands. . . between which 
are four deep and broad channels. The water flows together near 
the Pole, but at the Pole itself is a great Black Rock, 33 leagues 
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in circumference. Ships which once enter one of these channels 
never return, not even with the most favourable winds, and next 
to the Black Rock all the water is engulfed into the bowels of the 
earth, whence it flows through springs and river sources once 
again into the light of day. ['Translation.] 


From the Louisiane in 1702 Laval observed what he thought were 
variations in the strength of the Florida Current associated with the north 
component of wind; and in his book about the voyage published some 
years later (1728), he stated that he believed this phenomenon to be 
common knowledge among sailors. 


THE PERIOD 1700-1850 


During the early 1700’s the great American whale fishery sent ships all 
over the world, and the names of such patches of sand as Nantucket be- 
came known in every land touching the sea. The practical knowledge 
gained by these seafarers was not widely published in technical journals; 
instead, it was handed down by a system of apprenticeship and by word of 
mouth. And meanwhile, the basic understanding of fluid mechanics was 
growing. Daniel Bernoulli’s Hydrodynamica was published in 1738. In the 
last third of the century the study of the theoretical aspects of oceanic 
tides was brought to a high point in the contributions of Laplace. The 
influence of these developments in theoretical mechanics, and the general 
intellectual atmosphere of this Age of Enlightenment, discouraged further 
extraphysical and purely imaginative theories of ocean currents. 

In 1770 the Board of Customs at Boston complained to the Lords of the 
Treasury at London that the mail packets usually required two weeks 
longer to make the trip from England to New England than did the 
merchant ships. Benjamin Franklin was Postmaster General at the time 
and happened to discuss the matter with a Nantucket sea captain, Timothy 
Folger. The captain said he believed that charge to be true (Franklin, 
1786, p. 314): 


“We are well acquainted with the stream because in our pursuit 
of whales, which keep to the sides of it but are not met within it, 
we run along the side and frequently cross it to change our side, 
and in crossing it have sometimes met and spoke with those 
packets who were in the middle of it and stemming it. We have 
informed them that they were stemming a current that was 
against them to the value of three miles an hour and advised them 
to cross it, but they were too wise to be councelled [sic] by simple 
American fishermen.’ 
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Franklin had Folger plot the course of the Gulf Stream for him and then 
had a chart engraved and printed by the General Post Office. 
Franklin (1786, p. 315) believed that the Gulf Stream was caused 


by the accumulation of water on the eastern coast of America 
between the tropics by the trade winds. It is known that a large 
piece of water 10 miles broad and generally only 3 feet deep, has, 
by a strong wind, had its water driven to one side and sustained 
so as to become 6 feet deep while the windward side was laid dry. 
This may give some idea of the quantity heaped up by the 
American coast, and the reason of its running down in a strong 
current through the islands into the Gulf of Mexico and from 
thence proceeding along the coasts and banks of Newfoundland 
where it turns off towards and runs down through the Western 
Islands. 


By the time of Maury, in the middle of the nineteenth century (see 
Maury, 1859), Franklin’s estimates of the velocities of the Stream were 
regarded as excessive, but more recent studies tend to confirm them. 
Franklin did not give any details concerning the edge of the Stream. 
Starting in 1775, both Franklin and Charles Blagden (1782), independently, 
conceived the idea of using the thermometer as an instrument of navi- 
gation, and each made a series of surface temperature measurements while 
crossing the Atlantic. On Franklin’s last voyage in 1785 he even attempted 
to measure subsurface temperatures to a depth of about 100 ft., first with 
a bottle and later with a cask fitted with a valve at each end. 

A number of subsequent investigators made use of the surface thermo- 
meter, among them Governor Pownall (1787) and Captain Strickland 
(1802). It was in this fashion that Captain Strickland discovered a north- 
easterly extension of the Gulf Stream toward England and Scandinavia. 
These temperature measurements were not made with any idea of deter- 
mining the pressure field and geostrophic current, as is done today, but 
were simply regarded as an indication of the type of water through which 
the ship was sailing. 

In passing, the invention of the marine chronometer by John Harrison 
and its perfection by Thomas Earnshaw should be mentioned. By 1785, 
accurate chronometers were generally available to ships; this made the 
determination of longitude at sea at last a possibility and the deter- 
mination of the set of a current much more exact. Another important 
oceanographic tool, the drift bottle, was probably first used in 1802, when 
such a bottle was cast from the Rainbow. The use of drift bottles continued 
for more than a century, and finally received great impetus at the hands 
of the Prince of Monaco. 
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One important consequence of these early thermometric measurements 
was the discovery of pockets of cold water in the Gulf Stream. These 
pockets were first observed on the noteworthy cruise of the packet Eliza, 
en route from Halifax to England in April, 1810 (see Kohl, 1868), In the 
middle of the warm water of the Gulf Stream a mass of water colder than 
the surrounding water by about 10-15° Fahrenheit, and some 200 miles 
in diameter, was discovered. The explanation offered at the time was that 
these cold spots were due to the melting of icebergs entrapped by the Gulf 
Stream. This occurrence suggests a phenomenon very much like the large 
cold eddy reported in recent times by Iselin and Fuglister (1948), but its 
real nature and cause must remain in the realm of surmise. Another view, 
much in accord with modern synoptic experience, was that held by John 
Hamilton, who made serial air and water temperature measurements 
during twenty-six voyages to and from Europe, and asserted that the 
Gulf Stream was so unsteady and shifted its position so frequently that it 
was impossible to define its limits. Both Humboldt (1814) and Sabine 
(1825) were convinced that changes in the strength of the trades affected 
the Gulf Stream, and Sabine even suggested the use of weather ships 
which, he said, should observe the Florida Current and then sail quickly to 
Europe with the news of how strongly it was flowing, so that weather 
predictions could be made. 

In 1832 the results of an extensive compilation by James Rennell of data 
from the British Admiralty Office were published posthumously. Rennell 
distinguished clearly between ‘drift currents’, which are produced by the 
direct stress of the wind, and what he called ‘stream currents’, which are 
produced by a horizontal pressure gradient in the direction of flow. 
Rennell, in accord with Franklin’s earlier view, regarded the Gulf Stream 
as a current of the second kind and decided from his investigation that: 
(i) the breadth of the Stream changes from time to time; (ii) the breadth 
can vary as much as twofold even within so short a period as ten weeks; 
(iii) the variations are not seasonal; (iv) the north side of the Stream is 
more permanent than the south side; (v) temperature alone does not prove 
the existence of the Stream, for even warm countercurrents may exist; and 
(vi) cold-water inclusions occur within the body of the warm water. Also, 
Rennell proposed a special nomenclature for various parts of the Gulf 
Stream System. 

The work of Rennell was so authoritative and exhaustive at the time 
that it must have seemed to his contemporaries that the major features of 
the Gulf Stream had been delineated and charted. Few seriously enter- 
tained doubts that the Gulf Stream was a result of the downhill flow of 
water piled up by the trade winds along the American coast and in the 
Gulf of Mexico. 
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At this point a discovery was made which for a time completely dis- 
credited the notion that the winds are the cause of the ocean surface cur- 
rents. Arago drew attention in 1836 to the results of a leveling survey 
across Florida which showed a difference of not more than 7}in., and 
probably less. This difference of level seemed much too small to produce 
the Gulf Stream; hence Arago advanced the idea that the cause of the 
ocean currents is simply the density differences at the equator and poles 
due to unequal solar heating. ‘We should use the same theory for ocean 
currents which we use to explain the Trade Winds’, he stated. Seafarers 
were still inclined to the theory of wind-driven currents. Our present ideas 
of the magnitude of frictional forces in the Gulf Stream suggest, however, 
that a head of 74 in. is adequate to drive the Stream. 


THE LAST HALF OF THE NINETEENTH CENTURY 


In 1844 an intensive study of ocean circulation was begun by Matthew 
Fontaine Maury (1859), who was partial to Arago’s theory. He dismissed 
Rennell’s theory of wind-driven circulation on the following grounds: 
(i) assuming that the transport through the Florida Straits is the same as 
that of the current off Cape Hatteras, he (Maury) deduced the depth at 
Hatteras, completely ignoring the possibility that the Stream may mix 
with its environment; (ii) since the depth deduced in this manner is less 
than that at the Straits of Florida, Maury made the naive claim that the 
Gulf Stream would have to run uphill, and this reduced the wind theory to 
absurdity. He did not realize that the horizontal pressure gradient along 
the axis of the Stream is not controlled by the slope of the bottom boundary, 
so that if the top surface slopes down toward the north, the current can 
very well flow northward despite its diminishing depth. Maury was very 
much confused concerning fluid mechanics, even though that science was 
being rapidly advanced in his time; Stokes’s famous paper on the flow of 
viscous fluids was published in 1845, and the significant work by Coriolis 
had appeared in 1835. Maury uses such examples as the accumulation of 
sawdust in the center of a basin filled with water to explain the accumu- 
lation of sargassum in the Sargasso Sea, and, to illustrate the Coriolis 
force, the fact (?) that railroad trains usually run off the rails to the right— 
examples having no more than extremely doubtful applicability to the 
ocean. He goes on to ask why, if the Atlantic circulation is a completely 
closed circuit, there should be a piling up at any one place. This shows that 
Maury did not understand Rennell’s distinction between drift currents and 
stream currents. There seems to be little point in enumerating Maury’s 
misconceptions about the effect of the earth’s rotation on ocean currents, 
his roof-shaped current theory, his strange idea that waters of different 
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salinities are immiscible, and his suggestion that there is a kind of peri- 
staltic driving force acting upon the Gulf Stream. Maury says that the 
difference of density at the poles and the equator must produce the ocean 
currents, but he had no quantitative idea about how salinity and tem- 
perature are related to density. Accurate hydrographic tables had not 
yet been constructed, and the salinity distribution in the sea was only 
imperfectly known. 

The encouragement which Maury gave to the collection of ships’ data, 
however, his calling of the Maritime Conference at Brussels in 1853, and his 
dissemination of hydrographic information, in the form of good pilot charts 
and revised sailing directions, were of great practical and commercial 
importance. 

Modern surveying of the Gulf Stream began in 1844 with the work done 
by the United States Coast and Geodetic Survey under the direction 
of Franklin’s great-grandson, A. D. Bache (1860): fourteen temperature 
sections between Tortugas and Nantucket. These early survey cruises 
confirmed the existence of cold veins within the Stream. Bache supposed 
that the cold veins were a result of the bottom configuration which diverted 
the Gulf Stream in separate bands—irregularities which are now regarded 
not as straight bands, but as meanders. In Bache’s time the ‘bands’ of 
cold water were supposed to be invariable in number and position, and the 
USCGS Gulf Stream Chart of 1860 shows them so. The currents were 
observed only at the surface by the drift of the ship. 

During the Civil War there was a lull in Gulf Stream investigations, but 
by 1867 they were resumed by Henry Mitchell, who attempted to measure 
subsurface currents by means of two floats attached by a fine wire, one 
of which remained at the surface while the other sank slowly. Mitchell 
concluded that the velocity of the Gulf Stream off Fort Chorrera, Cuba, 
extended undiminished to a depth of at least 600 fathoms. For the next 
ten years the Bibb, the Bache, and the Blake continued to survey the 
Stream. Early soundings had been made by rope, but in 1881 Bartlett 
made use of the Thomson sounding machine, in which piano wire was 
used. Bartlett did not detect the cold and warm bands shown on Bache’s 
chart. According to the modern view of the temporary nature of cold 
inclusions in the Gulf Stream, this is not surprising. The Challenger, on her 
world-wide cruise, visited the Gulf Stream in 1873. 

John Elhott Pillsbury (1891) commanded the Blake during a remarkable 
series of observations beginning in 1885. The ship was anchored along 
several cross sections in the Florida Straits, and the current vector and the 
temperature at various depths were determined. The anchoring gear and 
the current meter were of Pillsbury’s own design. Pillsbury went about his 
task in a careful and painstaking manner. For example, he took two years 


HISTORICAL INTRODUCTION 9 


to occupy the section from Fowey Rocks to Gun Cay, ‘the time actually 
employed in observations at Section A being over 1,100 hours’. Obser- 
vations were also made at four other sections across the Florida Straits, 
at a station off Cape Hatteras, and at a number of stations in the passages 
of the Windward Islands. 

Pillsbury noted the westward intensification of streamlines within the 
Straits themselves, as earlier reported by Fremont. By sounding across 
the Gulf Stream all the way from the Straits to Cape Hatteras and showing 
that the bottom in this region was comparatively smooth, he disproved 
Bache’s hypothesis that mountain ranges under the sea cause ‘bands’ of 
cold water. His observations of temperature and current in the Florida 
Straits are unique and are so valuable that they are still used (Wiist, 1924) 
as the classical example of the accuracy of the geostrophic current deter- 
mination. He was a firm believer in the theory of wind-driven ocean 
currents. He felt that the time variations of temperature, horizontally 
and vertically, along the Gulf Stream were abundantly proved, and he 
attributed them to random local winds and tidal influences along the 
coast. He established that the Antilles Current is a tributary to the Gulf 
Stream. He emphasized the importance of the prevailing westerlies in 
maintaining the North Atlantic Drift. He surmised that the Stream 
becomes increasingly meandering as it flows northeastward. 

During the 1870’s there was a good deal of discussion on a qualitative 
plane between adherents and antagonists of the wind theory. The reader 
can form a fair idea of the level at which these discussions were carried on 
by reading the views of Aitken (1877), Carpenter (1874), and Sir C. Wyville 
Thomson (1874). A mathematical physicist, Z6ppritz (1878), attempted to 
show that the wind stress could cause appreciable ocean currents only after 
hundreds of thousands of years of constantly acting upon the water, 
because the molecular viscosity of water is so small. The importance of the 
role of turbulence in the ocean as an agency in the transfer of momentum, 
as well as of other properties, was not at that time appreciated. In 1883, 
five years after the work by Z6ppritz, Reynolds’ famous paper on his 
experiments on turbulent flow in pipes appeared. The realization of the 
fact that the ocean is essentially a turbulent regime showed the error in 
Zoppritz’ reasoning. 

It is important to realize that thus far in the nineteenth century the 
influence of the earth’s rotation upon ocean currents was only imperfectly 
understood by oceanographers, although the hydrodynamical equations of 
a perfect barotropic fluid relative to a rotating sphere had been written 
down in the previous century by Laplace (1778). The extremely important 
fact that the Coriolis force is almost everywhere balanced by the horizontal 
pressure gradients associated with the distribution of mass was not realized, 
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despite the fact that these terms are several orders of magnitude greater 
than the terms of the driving force and dissipative force. William Ferrel 
(1882) was one of the first to understand the role of the Coriolis force in the 
distribution of ocean currents caused by the wind. He derived the relation 
between the barometric gradient and the velocity of the wind (which is the 
counterpart of the geostrophic current relationship of oceanography). For 
some reason these remarkable achievements of Ferrel were overlooked by 
oceanographers, and therefore his work did not have a direct influence on 
the development of physical oceanography. The formula for computing 
ocean currents from the slope of isobaric surfaces was derived by Henrik 
Mohn (1885). 

The first careful study of the possibility of obtaining an approximation 
to the velocity field in the ocean from an exact knowledge of the mass field 
was made by Sandstrém and Helland-Hansen (1903), on the basis of 
Vilhelm Bjerknes’ (1898) circulation theorem. It became evident that the 
Coriolis force acting upon the Gulf Stream is counterbalanced by a hori- 
zontal pressure gradient associated with the mass of lighter water in the 
Sargasso Sea, and that the Stream is not so much a warm current as a 
boundary phenomenon associated with, but not caused by, the sudden 
downward slope of the isotherms toward the center of the ocean. It was 
recognized that the differences in density across the Stream have nothing 
to do with the driving of the Stream, but are simply part of an equilibrium 
brought about indirectly by the stress of the wind. The great step forward 
represented by this method of dynamic computation—as it came to be 
called—was that at last oceanographers had begun to use the hydro- 
dynamical equations of motion in their study of the sea, even though only 
in a piecemeal way. 

By the turn of the century the dependence of density of sea water on its 
temperature, salinity, and pressure had been determined. The ocean-wide 
distribution of salinity had been charted and plotted, mostly as a result of 
the pioneering work of Forchhammer (1865). This development came 
surprisingly late in synoptic oceanography, especially when one considers 
how much earlier the surface temperatures in most parts of the North 
Atlantic had been known. 


THE TWENTIETH CENTURY 


We now enter the modern period of physical oceanography in which the 
hydrodynamical equations of mean motion are used, although never in 
their complete form. There are several reasons for using only a few of the 
terms. First, the equations are nonlinear in their complete form, and 
therefore the mathematical methods for solving them are unavailable. 
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Secondly, the parameters depending upon turbulence are very poorly 
known. Thirdly, even today we do not have nearly as much accurate 
synoptic information about the Gulf Stream as we desire. 

V. W. Ekman (1905) developed a wind-drift theory, taking into account 
the effect of the earth’s rotation, which for the first time made possible 
computations of the effect of wind stress in transporting water. In Ekman’s 
treatment the ocean is regarded as homogeneous and of infinite extent 
horizontally. A constant wind blows over this ocean. In deep water the 
surface-current vector is 45° cwm sole to the wind and changes direction 
and diminishes with depth according to the Ekman spiral. 

Later, Ekman (1923, 1932, 1939) carried out further investigations of 
ocean currents in which he discovered certain deflecting effects of the 
inequalities in the depth of a homogeneous ocean. A qualitative study of 
the effect of bottom topography upon currents in a heterogeneous ocean 
has been advanced by Sverdrup, Johnson, and Fleming (1942). Certain of 
these deflections have been apparently observed (Neumann, 1940) in the 
northern parts of the Gulf Stream System associated with the Altair cone. 

The first intensive surveys of the Gulf Stream System by the research 
ketch Atlantis began in the early 1930’s, soon after the establishment of 
the Woods Hole Oceanographic Institution. This program was carried out 
under the leadership of C. O’D. Iselin. 

From this point on, I shall dispense with the historical narrative in 
favor of a topical exploration of the present knowledge concerning the 
Gulf Stream. It is not surprising that even now, after many years of effort, 
our conception of the Gulf Stream is incomplete. I shall try to show in 
what respects this knowledge is deficient, as well as those features which 
we do understand, and also endeavor to explain some of the greatest 
difficulties encountered in planning and executing surveys. 


Chapter Two 
METHODS OF OBSERVATION 


The obstacles encountered by the collector of hydrographic data pertaining 
to the Gulf Stream are numerous and formidable. Although many diffi- 
culties remain, modern technology has been of great assistance in oceano- 
graphic research. 


INSTRUMENTS USED IN GULF STREAM WORK 


The instruments utilized are not, for the most part, produced exclusively 
for the study of the Gulf Stream. The reversing bottle, the thermometer, 
and the bathythermograph are employed in every ocean area. However, 
the use of these instruments in the Gulf Stream does require special 
techniques. For example, anchor stations are out of the question in the 
strongest part of the Stream, and hydrographic stations are usually made 
by steaming against the surface current to reduce the wire angle in the low- 
speed water below. Because most readers are not likely to be familiar 
with these instruments, this chapter offers a brief description of them and 
of the techniques for their use at sea. 

Reversing bottle —The reversing bottle is made of metal and is fitted 
with a valve at each end. When lowered to the desired depth it is 
tripped by a messenger (a weight which slides along the wire). As it 
turns over, the valves close and a little more than 1 liter of water is 
sampled. Normally a number of bottles are lowered in series on the 
same cable. 

Salinity is determined by chemical titration, usually on shore. The 
standard of accuracy is about +0-02%, (%, means parts per thousand). 
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Reversing thermometers.—The reversing thermometer is mounted on the 
sampling bottle. The thermometer is so contrived that when it is turned 
upside down the mercury column breaks. When the thermometer is brought 
to the surface the reading must be corrected, because the temperature at 
which the mercury is broken is usually different from that at which it is 
read on deck. 

These thermometers are of two distinct types, which are used in pairs; 
that is, two thermometers, one of each type, are mounted on the same 
bottle. One type is protected from the hydrostatic pressure by an exterior 
glass shell (thermal contact with the outside is maintained by a pool of 
mercury); the other type is unprotected from the pressure, and hence 
registers depth because the mercury bulb is squeezed. The deviation 
amounts to about 0°01 C. per meter. The standard of accuracy of each type 
is + 0°01 C. Reversing thermometers and reversing bottles are designed for 
use with the ship hove to or steaming very slowly against the current or 
into the wind. They are necessary for making geostrophic current calcu- 
lations. 

Bathythermograph.—The bathythermograph is an instrument designed 
for use from a ship under way. A smoked-glass slide is mounted on a 
spring-loaded bellows so that it moves under variations in hydrostatic 
pressure. A bourdon thermal element moves a needle over the slide at 
a right angle to the direction of the pressure-induced motion. The scratch 
on the smoked slide is therefore a plot of temperature versus depth. The 
standard of accuracy in a model which works to a depth of 900 ft. is 10 ft., 
and +0°2C. Models for other depth ranges are also available. 

Bathypitotmeter.—A promising device for towing at low speeds, to make 
simultaneous recordings of the velocity of water relative to the instrument, 
the temperature, and depth, was devised by Malkus (1953). The record is 
made on a moving waxed-paper strip. This instrument is capable of making 
a determination of a complete vertical sounding in about two hours. The long 
length of cable, which is payed out when the ship is under way, apparently 
isolates the instrument from the irregular pitch and roll motions of the 
ship, which normally cause spurious readings on velocity meters suspended 
vertically from a drifting ship. 

Propeller-type current meters.—There are various kinds of current meters 
of the propeller type, from the early Ekman meter to the recent Watson 
meter. They record velocity and direction, but they take so long a time to 
operate that they are not often used for Gulf Stream studies. 

All velocity meters measure water velocity relative to the ship, but the 
ship’s velocity relative to the earth is seldom known adequately from 
navigation, even in areas with good Loran coverage (where radio fixes good 
to 1/2 mile are available at any time). Therefore there is always an 
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ambiguity concerning the zero point of all velocity measurements made at 
sea. Even anchored ships or buoys move about at the end of their cable 
so much as to be unreliable reference points in deep water. 

Towed electrodes—The component of a ship’s velocity relative to the 
earth and perpendicular to the ship’s heading may be determined by 
measuring the electromotive force induced by the earth’s magnetic field 
in a 100 m. length of wire towed astern. The complete vector velocity may 
be obtained by occasionally altering the ship’s course. As developed by 
William von Arx, this method has proved very valuable for practical 
survey purposes in indicating the presence and direction of surface currents 
(see Longuet-Higgins et al., 1954; and von Arx, 1950). The time required 
is about 5 min. per measurement, because the ship must be maneuvered in 
directions other than the regular heading. 

The chief difficulty of this method is that it does not give uniformly 
accurate surface velocities. The reason for this deficiency can be easily 
made clear by a concrete example. Assume that a ship and the cable 
behind it are swept sideways by the component of surface current velocity 
perpendicular to the ship’s heading. An electromotive force equal to vH,1 
is induced in the cable, where v is the sideways velocity of the ship and 
cable, H, is the vertical component of the earth’s magnetic field, and / is 
the length of the cable. If electrodes making good contact with the sea are 
fitted at both ends of the cable, and were there no potential induced in the 
sea by virtue of its motion, the electromotive force induced in the cable 
could be measured by a potentiometer and the velocity determined by the 
relation expressed above. In general, however, there is a potential gradient 
component developed in the ocean in the same direction as that in the 
cable; and since there is no way of knowing exactly what this ocean 
potential is except by knowing a great deal about the ocean velocity field, 
the method cannot give true surface velocity. In ocean current systems 
that are broad compared to the depth of the ocean, the potential gradient 
developed in the ocean is approximately 0H,, where 0 is the velocity 
averaged over the total depth. Since in the Gulf Stream the strongest 
currents extend to only a fraction of the total depth, <v, the towed 
electrodes give a good approximation to surface velocities, the deficiency 
varying between 2 and 40 per cent. 

Air-borne radiation thermometer.—Stommel and Parson (see Stommel et 
al., 1953) conducted some field experiments to determine the feasibility of 
making rapid synoptic surveys of the surface temperature over the Gulf 
Stream by airplane, using an infrared bolometer. The bolometer is exposed 
alternately to the sea surface and to a black body in the plane. The chopped 
radiation, falling on the bolometer, is restricted to the band of maximum 
water-vapor transmission (8—13 ~) which contains about 30 per cent of the 
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energy of the ocean’s emitted radiation. The airplane is flown below an 
altitude of 1000 ft., and preliminary tests show that readings accurate to 
1° F. are obtainable. The contrast in surface temperature across the Stream 
is about 4° F. in late summer, 20° F. in late winter. It is estimated that 
under favorable weather conditions for flying, the plane could make about 
six zigzag crossings of the Stream from Cape Hatteras to longitude 65° W. 
in a single 10 hr. flight. Such rapid surveys, which might be made twice 
a month, are still in the planning stage. 


LORAN 


Loran is a form of radio navigation developed during the Second World 
War by the Radiation Laboratory of the Massachusetts Institute of 
Technology. Its coverage is not world-wide, but fortunately there is a 
network of stations covering the central sector of the Gulf Stream System. 
A hyperbola of position is determined by accurately measuring (to 1 psec.) 
the difference in time between pulse signals sent at a constant rate from 
two shore stations which emit in perfect synchronization. A Loran fix is 
obtained by the intersection of several lines of position. Under the best 
circumstances the fix is accurate to within 1/4 nautical mile, and fixes can 
be made successively as often as every 10 min. This is a very great advance 
over the method of celestial navigation, which is usually limited to 
obtaining fixes at sunrise and sunset, and a latitude sight at noon. 


Chapter Three 


THE GEOSTROPHIC 
RELATIONSHIP 


It is probable that the reader unfamiliar with the literature and conven- 
tions of meteorology and oceanography will find that the manner in which the 
hydrodynamical equations are introduced and abbreviated in this and later 
chapters is neither sufficiently explained nor fully justified. Because of the 
specialized nature of the topic of this book it did not seem advisable to me to 
try to develop the formulation from first principles. The unsatisfied reader 
is therefore advised to refer to a good treatise on dynamical meteorology, or 
to Proudman’s (1953) textbook on dynamical oceanography, where an 
adequately detailed exposition can be found. Lamb (1932) gives a concise 
derivation of the equations of motion referred to a rotating reference frame. 


THE BALANCE OF FORCES IN THE GULF STREAM 


In the ocean there are two sets of forces which are nearly always almost 
balanced. The first of the near-balances, and perhaps the more obvious, is 
that in the vertical between gravity force! and the vertical pressure gradient. 
We introduce a system of rectangular codrdinates in the ocean, with the 
x-axis directed toward the east, the y-axis directed toward the north, and 
the z-axis directed vertically upward. The expression for this hydrostatic 


balance of forces is 
op a 1) 
pope (iD ( 


1 Gravity is defined as the vector sum of the gravitational force due to the earth’s 
mass and the centrifugal force due to its rotation. 
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where p is pressure, g is the acceleration due to gravity, and p is the 
density of the ocean water. In general, the order of magnitude of each of 
these two terms is 103 dynes/em.*. Under surface waves, where there are 
vertical accelerations, terms up to the order of magnitude of 10? dynes/cm.? 
may be expected, but these are important only in the upper 100 m. The 
acceleration term for ocean tides is much less, perhaps of the order of 
magnitude of 10-* dynes/cm.?. The vertical Coriolis force due to horizontal 
motions is not likely to exceed 10-* dynes/cm.*? anywhere in the ocean; 
and, because large-scale vertical motions are so small, this force may be 
neglected for practical purposes, such as in the computation of depth in 
terms of measured pressures. The balance of forces expressed by this 
equation is therefore correct to a very high order of approximation for 
most large-scale oceanic phenomena. 

The horizontal equations of motion for large-scale oceanic features 
exhibit a similar balance between two large forces compared to which most 
other terms in the equations of motion are often small. In order to illus- 
trate this, we shall first write the two equations of motion including only 
these two forces (the so-called geostrophic equations) : 


_ op 
ppo=s (2) 
Bee 
mila age (3) 


The symbols wu and v are the x and y components of velocity, and f is a 
quantity called the Coriolis parameter f=2w sin ¢, where w is the angular 
velocity of the earth and ¢ is the geographic latitude. In mid-latitudes the 
value of f is approximately 10-*/sec. The maximum current velocities in 
the Gulf Stream range from 100 to 250 cm./sec., and hence the Coriolis 
force, acting at right angles to the Stream, is of the order of 10-* dynes/cm.°. 
This Coriolis force is nearly balanced by horizontal pressure gradients 
due to the density distribution in the ocean, as shown in equations (2) 
and (3). 

It is interesting to compare the order of magnitude of these terms with 
that of other terms in the equations of motion. If there is a local ac- 
celeration of the ocean currents at a particular locality such that a 
250 cm./sec. current is diminished to zero in the course of a week, the local 
acceleration term is of the order of magnitude of 4 x 10-4 dynes/cm.’. 

If the current happens to be flowing in a curvilinear path, with a radius 
of curvature of &, the inertial terms will be of the order of magnitude of 
u2/R. For example, if the radius of curvature of the streamlines is 200 km., 
and the current velocity is 200 cm./sec., the inertial terms in the equations 
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of motion will be of the order of magnitude of 2 x 10-3 dynes/cm.3. This 
term sometimes approaches the Coriolis force and is not always negligible, 
especially in meanders of the Gulf Stream. As will be shown in Chapter 
VIII, the inertial terms in the direction of flow of the Gulf Stream are 
quite likely to be important, even though the cross-stream force com- 
ponents are essentially geostrophic. 

So little is known about the nature of the turbulent shearing stresses 
in the ocean that it is dangerous to attempt to estimate their order of 
magnitude. The Laplacian of the horizontal velocity in the Gulf Stream 
has been observed to reach values as high as 2 x 10-"'/cm./sec. Information 
concerning the order of the magnitude of the horizontal-eddy viscosity is 
scarce. We might take as a conceivable maximum value 5 x 10? cm?./sec. 
Under these circumstances the maximum viscous force in the Gulf Stream 
would be of the order of magnitude of 10-* dynes/cm.*. Slightly larger 
eddy viscosity could make the Stream appreciably nongeostrophic. 

Under normal conditions, the stress of the wind on the ocean surface is 
of the order of 1 dyne/cm.?. If this stress is distributed evenly over a 
layer of water 50 m. deep the effect of wind stress by vertical turbulence 
is of the order of 2 x 10-4 dynes/cm.’. Thus, all these terms are small com- 
pared to the Coriolis force and to the horizontal pressure gradients in the 
Gulf Stream. 

One sees, therefore, that in the horizontal equations of motion there is 
an approximate balance between the term of the Coriolis force and the 
term of the horizontal pressure gradient. This relation, often called the 
geostrophic equation, has been of practical use in estimating the velocities, 
and, more especially, the transports, of the Gulf Stream. The actual 
numerical process is somewhat involved, but is essentially based on the 
equations obtained by cross-differentiation of the geostrophic equations 
and the hydrostatic equation; elimination of pressure results in the 
following pair: 


O(rfp)_ _— ap 
O(ufp) op 


Thus, in an ocean current flowing toward the north (positive y-direction), 
dv/dz at each level is associated with a decrease of density toward the east 
at that level (in the Northern Hemisphere), and there is no appreciable 
horizontal gradient of density in the direction of the stream: du/éz=u=0. 
The density of sea water as a function of depth and position along a 
vertical section through the ocean thus provides us with a basis for com- 
puting the vertical gradient of horizontal velocity normal to the section. 
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Were the velocity known for any one value of z it would be possible to 
determine the velocity at all z by numerical integration with respect to z. 
In this way one calculates the geostrophic currents (often called, in oceano- 
graphic parlance, by the unfortunate solecism ‘dynamic currents’) from 
the density field. Meteorologists use the same principles for computing 
winds aloft from observed density fields, but they have the inestimable 
advantage of knowing the true pressure distribution at the earth’s surface, 
and thus there is no uncertainty in the constant of integration. Every- 
thing hinges on a proper determination of the constant of integration. In 
oceanography we speak of a ‘level of no motion’, that is, a value of z at 
which we can believe the velocity is zero, and from which we can carry 
out the vertical integration of differential equations (4) and (5) for the 
velocity field. The fact that this all-important level of no motion, or 
reference level, is still, to all intents and purposes, undetermined, is one of 
the most disconcerting features of physical oceanography. Some methods 
for estimating this reference level are given in Sverdrup’s textbook 
(Sverdrup et al., 1942, pp. 456-457). One school of thought simply relies on 
placing the reference level sufficiently deep to be below the most intense 
horizontal gradients of density. Thus Iselin’s (1940, p. 24) transport calcu- 
lations of the Gulf Stream are based upon an arbitrary choice of 2000 m. 
as the reference level, or level of no motion.? Were the chosen level signifi- 
cantly lower, say at the bottom, the resulting transports would be at least 
50 per cent greater; however, the velocities in the very surface layers of 
the Gulf Stream would not be changed much. Fortunately, the choice of 
the reference level has less effect upon velocities of surface water than on 
those of deep water. It is very important to remember that the 2000 m. 
reference level is completely arbitrary. Defant (1941) has proposed that 
the level of no motion coincides with the level of no vertical gradient of 
geostrophic velocity. Since he has been able to find levels of no vertical 
gradient in most of the ocean, he has drawn up charts of velocity based on 
this completely intuitive criterion. Finally, the use of continuity of mass 
and of conservation of various properties such as salt and heat content. 
has been proposed in special cases in which a section covers all possible 
entrances and exits for water in a closed arm of the sea. Fuglister and I 
have tried such calculations. They place a very severe load upon the 
accuracy of the observations and upon the assumption that there is no 
time variability or mixing in the structure of the deep water. They involve 


2 The depth of an observation level, or reference level, is sometimes given in terms 
of hydrostatic pressure rather than in linear measure. The decibar is nearly equivalent 
to the meter. I have used meters even when referring to works which use the decibar 
unit. Also, in the past oceanographers usually neglected the transports below the 
reference level. 
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very small differences of very large numbers. It is difficult to place much 
confidence in them. 

Thus the choice of the reference level for geostrophic calculations be- 
comes mostly a matter of taste, and we should admit that that is ultimately 
intolerable. The determination of the level of no motion is not a matter 
for debate, but for direct measurement—a subject to which I shall return 
in a polemical section at the end of this book. 


Fig. 1. Schematic cross section of the Gulf Stream. The arrows indicate 
direction of horizontal pressure gradient. The little circles indicate vanishing 
horizontal pressure gradient. The current is flowing perpendicular to the plane 
of the page. The dashed lines are contours of equal water density. The line 1 
is a level surface; the line s is the actual sea surface. 


SCHEMATIC DENSITY AND PRESSURE FIELDS 
ACROSS THE GULF STREAM 


Fig. 1 is a schematic diagram showing a hypothetical cross section of the 
Gulf Stream. A level surface (such as the ocean at rest would assume) is 
shown by a heavy solid line. The hypothetical sea surface across the Gulf 
Stream is shown by the thin solid line. The difference in elevation between 
one side of the Stream and the other is supposed to be about 1 m. The 
vertical scale of the diagram is exaggerated very greatly to show the shape 
of the free surface. The broken lines beneath the surface represent contours 


THE GEOSTROPHIC RELATIONSHIP P| 


of equal density, increasing downward. The spacing of these lines is also 
vertically exaggerated, but not so much as the free surface. The difference 
in level of the broken lines is about 700 m. across the Stream. 

The diagram extends to a depth of about 1000 m., the lower 3000 m. 
being omitted. The arrows show the hypothetical magnitude and direction 
of the horizontal pressure gradient across the Stream. Near the surface 
the pressure gradient is controlled entirely by the shape of the free surface. 
At depth, the slope of the density surfaces takes effect until at some depth 
(in the diagram, 1000 m., but usually assumed* to be 2000 m.) the hori- 
zontal pressure gradients vanish. 

In order that such a distribution of density can long persist, these 
pressure gradients must be opposed by an opposite and equal force. The 
essence of the geostrophic method is that we suppose that these pressure 
gradients are opposed by Coriolis forces acting to the right of the direction 
of motion of the water. This implies velocities in a direction perpendicular 
to the plane of the figure, and directed into the page where the arrows 
point left, and outward where they point right. Since these forces are 
perpendicular to the direction of motion they neither drive nor brake the 
motion. 

The Gulf Stream is not a river of hot water flowing through the ocean, 
but a narrow ribbon of high-velocity water acting as a boundary that pre- 
vents the warm water on the Sargasso Sea (right-hand) side from over- 
flowing the colder, denser waters on the inshore (left-hand) side. 


3 See footnote 2, on p. 163. 


Chapter Four 


LARGE-SCALE FEATURES OF THE 
NORTH ATLANTIC CIRCULATION 


For more than a century the surface features of the oceans have been 
crudely determined from literally millions of ship reports sent in to the 
various naval hydrographic services throughout the world. Since no 
special effort is made to navigate with precision in the deep ocean, the 
reported surface currents are not determined with very high precision, but 
the charts drawn up from this large source of data give a good over-all 
view of surface currents. An example of such a compilation is the United 
States Navy Hydrographic Office Current Atlas (1946). 


SURFACE FEATURES 


A schematic chart of surface currents is shown in fig. 2. The part of the 
Atlantic on which our attention is focused in this book is the remarkably 
intense set of currents on the western side of the ocean, along the open 
coast of North America. People commonly speak of this whole system of 
currents as the Gulf Stream, even though very little water from the Gulf 
of Mexico is actually in the Stream. 

Iselin (1936, pp. 73-75) attempted to introduce a well-defined nomen- 
clature for various parts of the current system of the western North 
Atlantic. The entire set of western currents was to be called the Gulf 
Stream System, and the current from Tortugas, in the Florida Straits, to 
Cape Hatteras was to be called the Florida Current. The old term Gulf 
Stream was retained for the section of the current between Cape Hatteras 
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and the tail of the Grand Banks. Extensions of the Gulf Stream System to 
the eastward were to be spoken of as the North Atlantic Current. The North 
Atlantic Current, which appears to be made up of a number of separate 
streams, eddies, or branches (exactly which we do not yet know) is often 
obscured by a shallow, wind-driven surface movement called the North 
Atlantic Drift, which varies from time to time, depending upon the winds. 

In spite of the advantages of this nomenclature, no one has strictly 
adhered to it. In this book I often use the term Gulf Stream in a more 


Fig. 2. Chart showing the chief features of the surface-water circulation of 
the North Atlantic circulation, according to Sverdrup, Johnson, and Fleming 
(1942, fig. 187). In general, the chart is much oversimplified, and it should be 
regarded as essentially schematic. 


general sense than that proposed by Iselin; and I do not speak of the 
Florida Current as extending to Cape Hatteras, but restrict the use of this 
term to mean the current actually within the Florida Straits. Unfor- 
tunately, the naming of things is more a matter of common usage than of 
good sense. 

The names of other currents referred to in this book are also shown in 
fig. 2. 
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In order to discuss the Gulf Stream in relation to the broad features of 
the North Atlantic Ocean it is necessary to introduce at this point a 
number of charts and graphical presentations of data which will be useful 
in later descriptions and for reference. Thus there are presented, in figs. 3 


SURFACE TEMPERATURE 


DEGREES FAHRENHEIT 


FEBRUARY 


Fig. 3. Contours of surface temperature in the western North Atlantic for 
February, according to Fuglister (1947, pl. 2). 


and 4, the surface-temperature charts prepared by Fuglister (1947, pls. 2 
and 8) for the months of February and August. The currents and hori- 
zontal temperature gradients are actually much more pronounced than 
these average charts show, because the averaging processes used in their 
construction tend to blur the fine and shifting detail of the instantaneous 
velocity and temperature fields. 
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It was the close correspondence of the current charts (see fig. 2) and 
the surface-temperature charts (see figs. 3 and 4) which stimulated the 
interest of early seafarers in ‘thermometric’ navigation so many years 
ago (see Chapter I). 


SURFACE TEMPERATURE 


DEGREES FAHRENHEIT 


AUGUST 


Fig. 4. Contours of surface temperature in the western North Atlantic for 
August, according to Fuglister (1947, pl. 8). 


THE DISTRIBUTION OF PROPERTIES WITH DEPTH 


An average temperature sounding (with depth) for two seasons of the year, 
and on both sides of the Gulf Stream, is shown in fig. 5,a. The curve marked 
‘slope’ is the sounding made in the slope water between the continental 
shelf along the coast and the Gulf Stream near Chesapeake Bay; the curve 
marked ‘central’ corresponds to a sounding in the North Atlantic Central 
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Water on the offshore side of the Stream. Fig. 5, b, shows the corresponding 
values of salinity. One sees that there are three clearly defined thermal 
regions in the ocean: (i) a surface layer several hundred meters in depth 
which is subject to seasonal thermal variations; (ii) a zone of marked verti- 


DEPTH IN KILOMETERS 
DEPTH IN KILOMETERS 


a b 


Fig. 5. Temperature and salinity soundings. a, temperature soundings at 
two different seasons of the year in slope and central water. The depth scale is 
in kilometers. The temperature is in degrees Centigrade. b, salinity soundings 
at two different seasons of the year in slope and central water. The depth scale 
is in kilometers. The salinity unit %, means parts per thousand. 


cal temperature gradient, called the main thermocline; and (iii) a very large 
mass of cold deep water below 1500 m. 

Fig. 6 shows the average volume transports of the various surface cur- 
rents in the North Atlantic Ocean in millions of cubic meters per second 
according to Iselin (1936). Fig. 7 isa block diagram of the gross features of 
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the North Atlantic Ocean. The observer is standing high above the Gulf of 
Mexico looking toward the northeast, that is, toward England. The ocean 
is dissected to show its structure with depth as well as along the surface. 

Block 1 represents the westernmost part of the North Atlantic, thus 
including the Gulf Stream. Block 2 represents the Sargasso Sea. Block 3 
represents the area of the North Atlantic Current, its eddies, and its 


Fig. 6. Iselin’s sketch showing sources (broken lines) and pattern (solid 
lines) of the Gulf Stream System. In the western half of the ocean each trans- 
port line represents about 12 x 10° m.3/sec. From Iselin (1936, fig. 48). 


multiple currents. The deep water of the ocean, which is indicated by the 
darkest shading, apparently does not circulate as rapidly as the surface 
waters. With an average age of some several centuries, it is only slowly 
renewed by sinking, in very limited areas such as that indicated by the 
crooked arrow in Block 4. Water between 5 and 16° C. lies mostly in the 
region of the main thermocline and is shown by intermediate shading in 
all the blocks. We have no direct information about its sinking rate, but 
there are some indications that this intermediate water is formed in the 
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northern half of Block 3, probably mostly in the wintertime. The actual 
transports of water at various locations in the Atlantic, as inferred from 
geostrophic calculations based on actual data, are discussed in more detail 
in Chapter XI. The thickness of the intermediate layer is remarkably con- 
stant throughout the whole North Atlantic Ocean; this is in marked con- 
trast to the thickness of the surface layer (water at a temperature higher 
than 17° C.), which varies widely from place to place. It seems reasonable 
to suppose that the water sinking into this intermediate layer in northern 
latitudes (where it reaches the surface in the winter) is eventually mixed 
upward into the surface layer in subtropical latitudes. This mixing occurs 


Fig. 7. A dissected-block diagram of the thermal structure and circulation 
of the North Atlantic Ocean, as viewed from a great height over, the Gulf of 
Mexico. The unshaded part is the warmest water; the lightly shaded part in- 
dicates the water of the thermocline; and the heavily shaded part represents 
cold deep water. The smooth curved arrows indicate the direction of flow of the 
horizontal currents of the surface and thermocline; the zigzag arrows indicate 
hypothetical slow vertical flows. A description of the individual blocks is 
given in the text. 


mostly in the fall and winter of the year, when cooling and wind stirring of 
the surface layer are at a maximum. To complete the cycle one must sup- 
pose that part of the surface water eventually finds its way back into the 
northern half of Block 3 through eddies and multiple streams and that 
it is there reconverted into intermediate water. A large fraction of the 
surface water, however, circulates without transformation in the hori- 
zontal wind-driven surface gyre. 

Water-mass analysis —Sverdrup, Johnson, and Fleming (1942) have dis- 
cussed, in very concise form, the nature and probable origin of various 
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water masses present in the North Atlantic Ocean. A more comprehensive 
and detailed study of the area can be found in the earlier monograph by 
Iselin (1936). Inasmuch as the subject of the present book is limited to the 
Gulf Stream System itself, it would probably be too much of a digression 
to try to discuss in detail the features of the entire North Atlantic Ocean. 
Therefore, the description given here is very brief. The reader who is 
interested in further detail is referred to Iselin (1936) and to the magni- 
ficent oceanographic atlases published in the scientific reports of the 
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Fig. 8. Plot of temperature versus salinity of selected station data from the 
North Atlantic, according to Sverdrup, Johnson, and Fleming (1942, fig. 183). 
A=Atlantis; AH=Armauer Hansen; G=General Greene; M= Meteor; 
N.A.=North Atlantic Central Region; M.W.= Mediterranean Water. 


German Meteor Expedition. Albert Defant’s study (1941) is one of the 
reports in this series. The two great atlases are Band 5—Atlas (Bohnecke, 
1936) and Band 6—Atlas (Wiist and Defant, 1936). 

The standard procedure used in water-mass analysis is to plot the hydro- 
graphic data from a number of representative soundings on a graph in 
which the coérdinates are temperature and salinity. This graph is custo- 
marily called a T-S diagram. Fig. 8, from Sverdrup ef al. (1942, fig. 183), 
shows the T-S diagram for a number of specially selected stations in the 
North Atlantic Ocean. Observations from the upper 100m. have been 
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omitted because they usually show a great scatter, being variable for 
different seasons. The depths of the shallowest points in the diagram are 
indicated by numerals (meters), and the insert chart shows the positions . 
of the various stations. 

The T-S diagram is used to identify different kinds of water in the 
ocean, to note similarities of properties between one region and another, 
and to obtain a qualitative picture of the amount of mixing between water 
masses of different geographic location. Four considerations must be kept 
in mind in the use of T-S diagrams. 

First, potential temperature and salinity are conservative properties in 
the ocean, except in a shallow homogeneous layer near the surface, where 
evaporation, solar radiation, cooling by back radiation, an exchange of 
heat with the atmosphere, and similar changes occur. Sea water is com- 
pressible; therefore, true temperature is not strictly conservative. The 
change in temperature of a water parcel that is subjected to an adiabatic 
vertical displacement of 1000 m. is roughly 0°1 C., and therefore, for many 
practical purposes of graphing, the potential temperature and the true 
temperature on a T-S diagram are interchangeable, although potential 
temperature is preferable. 

Secondly, mixing processes in the interior of the ocean mix both tem- 
perature and salinity in the same way. Thus the mixture of two masses 
of water represented by two points on a T-S diagram lies along a straight 
line joining them. 

The third consideration is that the oceans of tropical and temperate 
regions are stratified and vertically stable at all depths, with the possible 
exceptions of the shallow homogeneous surface layer and of certain deep 
isolated basins on the ocean bottom. Because of the damping of vertical 
turbulence by the stability, there is a tendency for most major flows to 
occur along surfaces of equal potential density. The potential density is 
the density that a sample of water would have if it were brought adia- 
batically to the surface at atmospheric pressure. It is important to make 
the distinction between density and potential density, because of the con- 
siderable compressibility of sea water. At a pressure corresponding to that 
present at the average bottom depth of the ocean, 4000 m., water is com- 
pressed about 2 per cent. Deep flows tend to occur along surfaces of equal 
potential density; hence it is important that the potential density be 
represented on the T-S diagram, to indicate the directions of preferred 
flow and mixing. The family of smooth curved lines on the T-S plane in 
fig. 8 represent loci of constant o;, a quantity which, unfortunately, does 
not have any special name, but which is very nearly the same as potential 
density. Because of its importance to physical oceanography it must be 
defined carefully. 


LARGE-SCALE FEATURES 31 


The density of sea water p is always a little greater than unity in the 
c.g.s. system; therefore, in order to avoid having to write long decimals all 
the time, the quantity o has been introduced, and is defined in the follow- 
eel o=(p—1) 1000. (1) 
Thus, a density of 1-:02721 is written as a o of 27-21. 

The quantity o, is defined as the o of sea water at atmospheric pressure 
and at the temperature at which it was collected. Therefore, it differs 
slightly from the o which would be computed from the potential density 
of the same specimen, in that no account is taken of the adiabatic tem- 
perature change of the specimen during reduction of the pressure to 
standard pressure. This difference is very small. The essential effect of 
compressibility of sea water is taken into account by the use of o, surfaces 
on a T-S diagram, instead of the rigorously correct equal-potential-density 
surfaces. 

The relation of density to temperature, salinity, and pressure has been 
the subject of elaborate and highly precise laboratory measurement. 
Accurate tables have been prepared. They are discussed in more detail by 
Sverdrup e¢ al. (1942). A very much abbreviated table for o;, is given in 
Appendix III of the present study. 

Finally, it is necessary to make one further remark about mixing along 
a surface of equal potential density. Since mixture on a T-S diagram 
occurs along straight lines, whereas the lines of equal potential density (or 
equal o,) are curved convexly upward, the mixture of two water masses of 
the same potential density tends to increase the potential density of the 
mixture. Although this effect, called cabbeling, should be borne in mind, 
it is not of great importance in the crude qualitative analysis that is usually 
done with T-S diagrams. 

North Atlantic water masses.—The water in the North Atlantic Ocean is 
made up essentially of two water masses: one, the so-called North Atlantic 
Central Water, and the other, the North Atlantic Deep Water. 

The North Atlantic Central Water is that water with a temperature 
between 8 and 19°C. and a salinity of between 35-10 and 36-70%. (The 
symbol %, is read ‘parts per thousand’.) This water mass is delineated by 
the two sloping broken lines on the T-S diagram (fig. 8). As can be seen, 
most of the soundings for temperature higher than 8° lie between these 
narrow limits. Iselin (1936, figs. 22 and 25) has drawn T-S diagrams for 
water masses in the western North Atlantic which show even sharper 
definition, indicating that over much of the Sargasso Sea the salinity of 
water of a given temperature does not vary much more than the accepted 
limits of error of salinity determination accommodate. The sounding 
which lies farthest from the mean in fig. 8 is that for Atlantis station 
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1175, which appears to represent South Atlantic conditions (low average 
salinity). 

The North Atlantic Deep Water is characterized by temperatures be- 
tween 3-5 and 2°2 C. and salinities between 34-97 and 34-90%,. This body 
of water composes most of the North Atlantic by volume. It is extremely 
homogeneous over the entire North Atlantic. It does not take part in the 
wind-driven surface circulation, but does exhibit large transports in a 
narrow stream near the western coasts of both the North and the South 
Atlantic—presumably as a result of the thermodynamically driven part 
of the circulation (see Chapter XI). It apparently originates in the most 
northerly parts of the wintertime North Atlantic, but the exact way in 
which it is formed is something of a mystery (see Worthington, 1954a). 
Near the very bottom it is probably mixed somewhat with bottom water 
of Antarctic origin. 

Between these two water masses, which may be regarded as the principal 
masses of the North Atlantic, there are other, smaller, amounts of water 
which are produced by mixing at intermediate depths along surfaces of 
equal potential density (07, range: 27-2—-27-8) with water from the South 
Atlantic and with an outflow from the Mediterranean Sea. The highly 
saline water from the Mediterranean is shown by the point M.W. in fig. 8. 
As the Mediterranean Water leaves the Straits of Gibraltar, it sinks along 
the 27-6 o, surface and mixes with water over most of the eastern North 
Atlantic. Its presence is obvious in both the two southern Armauer Hansen 
stations. Similarly, there is a mixture at mid-depths with water from the 
South Atlantic Ocean, called Antarctic Intermediate Water. This water 
mass, of low salinity, forms at the surface of the South Atlantic in a broad 
band extending from the Cape of Good Hope to Cape Horn. It then sinks 
beneath the surface and flows northward; and eventually it crosses the 
equator and mixes with waters in the North Atlantic. The particular 
stations chosen for fig. 8 do not show this minimum-salinity water very 
well except at Atlantis station 1175, which is really characteristic of South 
Atlantic conditions. A very small amount of low-salinity Arctic Inter- 
mediate Water is formed in the northern regions of the Atlantic Ocean; 
the effect of this mass shows up in General Greene station 1990. 

So far as the study of the Gulf Stream System is concerned, the most 
important water mass is the North Atlantic Central Water, since it 
occupies the upper 1000 m. of the central regions of the North Atlantic and 
takes part in the wind-driven surface circulation which gives rise to the 
Gulf Stream. In particular, there is a very large mass of nearly homo- 
geneous water at temperatures between 17 and 18° C., as shown in fig. 5, a. 
It is not possible to explain the origin of this water satisfactorily on a 
quantitative basis. Iselin (1939) calls attention to the fact that most of 
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the o, surfaces in the North Atlantic Central Water intersect the ocean 
surface at points where the surface T-S relation is similar to that at greater 
depths. He suggests that the surface waters sink along a; surfaces, without 
mixing, preserving the T-S characteristics acquired at the surface. Sverdrup 
(Sverdrup ef al., 1942, p. 145) has suggested that vertical mixing is likely 
to be of considerable importance in the area. It is possible (see Worthing- 
ton, 1954a) that the entire treatment of the origin of the North Atlantic 
Central Water to date has suffered from attempts to imagine it as a 
stationary phenomenon, in which every water-producing process acts all 
of the time. Upon closer scrutiny every scheme of circulation that has 
been suggested exhibits some shortcoming. In addition to the subtleties 
involved in the thermodynamical processes of a continuously stratified 
fluid, slow density flows, and imperfectly understood vertical and horizontal 
turbulent processes, there are a number of constraints of a dynamical 
sort (p. 122), associated with the wind-driven horizontal circulation of the 
North Atlantic, which probably play very significant roles in the forma- 
tion of the North Atlantic Central Water. 

A fair idea of the exchange of heat and water across the surface can be 
computed by semiempirical laws from a knowledge of average air and water 
temperatures (Jacobs, 1942), but even with this detailed information of 
net heat energy and water flux at the surface, the problem of the formation 
of the North Atlantic Central Water defies quantitative explanation. 

Water masses of the Gulf Stream.—As is shown in the current chart 
(fig. 2), most of the water which enters the Gulf Stream System is water 
previously driven westward by the trade winds. The westward flow is a 
broad band of moving water, called the North Equatorial Current, which 
moves slowly, and, especially on its southern side, is very shallow (200 m.). 
In reaching longitude 60° W. it divides its flow into two parts: one part 
flows through the Caribbean, in a series of gradually narrowing channels 
and straits, and eventually finds its way out into the North Atlantic 
through the Florida Straits; the other, northern, part flows north of the 
West Indies, then joins the Florida Current over about 8° of latitude. The 
combined flow, the mature Gulf Stream, leaves the coast at Cape Hatteras. 
The total transport of the Gulf Stream off Chesapeake Bay (e.g., in April, 
1932, 82 x 10° m.3/sec., assuming a 2000 m. reference level) also includes 
some water which recirculates in a long quasi-elliptical orbit to the south- 
east of the Stream, but which can hardly be properly called North 
Equatorial Current water. The transports of these various currents and 
parts of currents are shown in fig. 6. It is clear that since the various 
portions of water which make up the Gulf Stream proper come from a 
wide range of different localities in the North Atlantic, they bring to the 
Stream rather distinctive T-S relations, and that these may be used, to 
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some extent, to trace the origin of various kinds of water in the Stream. 
It would be premature to explore this subject further at this point of the 
book; it is therefore deferred until after the nature of the Stream itself 
has been discussed in greater detail. A recent section made by Worthing- 
ton from Nova Scotia to South America along approximately 64° W. 
longitude, sketched in fig. 9, illustrates how narrow the current of the Gulf 
Stream is compared to the Equatorial and Caribbean currents that feed it.1 
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Fig. 9. Vertical section along approximately 64° W. longitude, extending 
from Nova Scotia, on the left, through Bermuda and then through Mona 
Passage in the West Indies, and through the Caribbean to the coast of South 
America. Drawn from deep stations made by Worthington on the Atlantis 
and Caryn in 1954. This sketch shows particularly well the contrast between 
the sharp Gulf Stream, in the northern half of the section, and the broad North 
Equatorial Current in the southern half. Mr Worthington very kindly let me 
use some of his observations to construct this figure in advance of his publi- 
cation of the full details and final interpretation of his work. In deference to his 
privilege of prior publication this figure is only a sketch. It is not definitive. 


The density field of the Atlantic—Because the three-dimensional density 
field is fundamental in describing the geostrophic currents of the ocean, 
and will later be idealized in theoretical models, and also because equal 
o, surfaces are preferred directions for mixing and movement of water 
masses, it is worth while to reproduce here the charts of the depths of 
various o, surfaces drawn by Montgomery and Pollak (1942) from the 
Meteor data. These are given in figs. 10-15. 

1 Tn deference to Mr Worthington’s privilege of first publication of the data for this 


new section which he has obtained, I present in fig. 9 only a rough schematic sketch of 
the section. 
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Fig. 10. Depth in hundreds of meters of the 26-5 o, surface, according to 
Montgomery and Pollak (1942, fig. 12). 
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Fig. 12. Depth in hundreds of meters of the 27:2 7, surface, according to 
Montgomery and Pollak (1942, fig. 14.) 
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Fig. 13. Depth in hundreds of meters of the 27:4 o, surface, according to 
Montgomery and Pollak (1942, fig. 15). 
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Fig. 14. Depth in hundreds of meters of the 27-6 o; surface, according to 
Montgomery and Pollak (1942, fig. 16). 
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Fig. 15. Depth in hundreds of n 
Montgomery and Pollak (1942, fig. 17). 
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Bathymetric features along the Stream.—The course of the currents in the 
Gulf Stream System seems to be determined in part by the submarine 
topography of the western North Atlantic. Fig. 16 is a bathymetric chart 
of the region according to Tolstoy (1951, pl. 1). Fig. 17 shows the bathy- 
metry through the Straits of Florida, obtained from Hydrographic Office 
charts. The depths in these two figures are given in fathoms (1 fathom = 6 ft.). 
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Fig. 16. Bathymetric chart of the western North Atlantic in contours of 
500 fathoms, according to Tolstoy (1951, pl. 1). Sea mounts are shown in the 
upper right quarter of the chart. 


The Stream continues on a shelf of about 800 m. depth along the Blake 
Plateau to about 33° N., where it leaves the shelf. From Cape Hatteras 
northeast, the Stream flows through a region that is about 4000-5000 m. 
deep. Since the high-velocity part of the Stream does not penetrate much 
below 1500 m., it would be difficult to see how bottom topography could 
influence the Stream here, were it not for the recent discovery (by echo 
sounder) of numerous sea mounts in the area. Fig. 16 shows the position 
of some of these. 

Farther east the North Atlantic Current encounters the Mid-Atlantic 
Ridge, which, over large areas, comes to within 2000 m. of the surface. 

In order to aid the reader in locating the various sections and positions 
referred to throughout the text, fig. 18, an index chart has been prepared. 
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Fig. 17. Bathymetric chart of the Straits of Florida in 500-fathom intervals, 
but also showing 10- and 100-fathom curves. 
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BERMUDA 


Fig. 18. Index chart, showing the positions of various sections referred to in 
figures in this book. The numerals on the lines are the figure numbers of the 
corresponding sections. 


Chapter Five 


THE HYDROGRAPHY OF 
THE GULF STREAM 


The first detailed series of hydrographic stations across the Gulf Stream 
was that begun in 1931 by the Atlantis, then the only seagoing research 
vessel of the newly founded Woods Hole Oceanographic Institution, and 
consisting of soundings repeated quarterly for a number of years. 


‘ATLANTIS’ SECTIONS OF THE GULF STREAM, 1931-1939 


Many early sections were made between Bermuda and Chesapeake Bay 
(Iselin, 1936) and between Bermuda and Montau’ Point, New York 
(Iselin, 1940). Temperature sections were made at four different seasons, 
as shown in figs. 19-22; and salinity sections at the same seasons, as shown 
in figs. 23-26. In these figures the vertical scale above the depth of 
2000 m. is much exaggerated (a distortion of 1:370); but below 2000 m. it 
is less so (1:148). The warmer the water, the less dense. The most striking 
feature of all these sections is the pronounced change in level of the 
isotherms in a narrow region. According to the geostrophic relation 
[Chapter III, equations (4) and (5)], this narrow zone is where the high 
current velocities, perpendicular to the plane of the page (in figs. 19-26), 
occur. The surface of the waters to the left of the Stream, called slope water, 
is subject to wider seasonal fluctuations than the Sargasso water to the 
right of the Stream, where the primary seasonal change is the appearance 
of a shallow thermocline in the summer. 
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BATHYTHERMOGRAPH SECTIONS 


Fig. 27 shows the kind of detail available in a bathythermograph crossing 
of the Stream. In discussing the results of a multiple-ship cruise called 
‘Operation Cabot’, Fuglister and Worthington (1951, p. 3) have suggested 
the following definition of terms: 


Early in the planning and operational stages of Operation 
Cabot it became evident that precise definitions were needed for 
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Fig. 19. Temperature section across the Gulf Stream, Chesapeake Bay to 
Bermuda, February 11—18, 1932, according to Iselin (1936, fig. 3). 
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the various terms used in association with the Gulf Stream. The 
frequent references to the ‘cold wall’, ‘edge of the Stream’, 
‘warm core’ and ‘front’ led to a certain amount of confusion and 
misunderstanding. The term ‘inner edge’ was most frequently 
used and most variously interpreted. This confusion is caused 
primarily because, although the words ‘Gulf Stream’ denote a 
current, they also imply a distinct water mass, and secondarily 
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because water masses that may be motionless are included as part 
of the Stream because they lie below a surface current. 

Since the Gulf Stream is a boundary or front in the western 
North Atlantic between the slope water and the Sargasso Sea we 
may define it as follows: it is a continuous band stretching from 
the continental shelf off Cape Hatteras to the 50th meridian of 
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Fig. 20. Temperature section across the Gulf Stream, Chesapeake Bay to 
Bermuda, April 17—23, 1932, according to Iselin (1936, fig. 5). 


longitude, south of the Grand Banks of Newfoundland. This band 
consists of a pronounced pressure gradient between the warm, 
highly saline water to the south, and the colder, fresher water to 
the north. Using this definition then, the inner and outer limits 
or edges of the Gulf Stream can be defined as the points where 
this pressure gradient becomes zero. These points can be located 
only if deep, closely spaced temperature and salinity data are 
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Fig. 21. Temperature section across the Gulf Stream, Chesapeake Bay to 
Bermuda, August 28-September 3, 1932, according to Iselin (1936, fig. 7). 
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Fig. 22. Temperature section across the Gulf Stream, Chesapeake Bay to 
Bermuda, November 30—December 5, 1932, according to Iselin (1936, fig. 10). 


HYDROGRAPHY OF THE STREAM 


obtained, and the cross-current pressure gradients calculated. 
Also, because of the large eddies found both north and south of 
the Stream, any section made across the area must be long enough 
to ascertain whether or not more than one pronounced pressure 
gradient exists. If only one is found it defines the Gulf Stream, 
but if more than one are located then the position of the Stream 
cannot be determined by that single section. 
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Fig. 23. Salinity section across the Gulf Stream, Chesapeake Bay to 
Bermuda, February 11-18, 1932, according to Iselin (1936, fig. 4). 


Not to be confused with the inner or left-hand edge of the 
Stream is the temperature-salinity boundary at the surface. This 
generally abrupt change that occurs to the left of the ‘warm core’ 
may or may not coincide with the left-hand edge of the Gulf 
Stream as defined above. This applies also to the color boundary 
and the long thick lines of Sargassum frequently seen on the 
surface ; all of these surface phenomena are apparently associated 
with shear zones to the left of the ‘warm core’ but they are not 
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necessarily coincident with the left-hand or inner edge of the Gulf 
Stream. 

The ‘warm core’ is defined here as that part of the Gulf Stream 
containing water warmer than the water at the same depth to the 
right, facing down stream, of the current. This ‘warm core’ is 
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Fig. 24. Salinity section across the Gulf Stream, Chesapeake Bay to 
Bermuda, April 17-23, 1932, according to Iselin (1936, fig. 6). 


generally 300 to 400 meters deep with the maximum temperature 
anomalies at a depth of about 100 meters. 

The word ‘front’ is considered synonymous with the pro- 
nounced pressure gradient and therefore with the Gulf Stream 
itself. 

The term ‘cold wall’ dates back to 1845 and is still frequently 
used to denote the ‘inner edge of the Stream’ or, according to 
Church (1937), ‘the temperature gradient between the slope 
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Fig. 25. Salinity section across the Gulf Stream, Chesapeake Bay to 
Bermuda, August 28—September 3, 1932, according to Iselin (1936, fig. 8). 
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Fig. 26. Salinity section across the Gulf Stream, Chesapeake Bay to 
Bermuda, November 30—December 5, 1932, according to Iselin (1936, fig. 11). 
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water and the Gulf Stream’. According to these definitions it 
could equally well be called the ‘warm wall’ though in neither 
case do we have anything resembling a wall. This temperature 
gradient exists at different depths across the entire width of the 
Gulf Stream and therefore cannot be considered as something 
separate or adjoining the Stream. Because of the misleading con- 
notations of the term ‘cold wall’ it will not be used in this paper. 


800 
Temperature (°F.) 
900 


Fig. 27. Sample of temperature sections, in shallow surface layers (upper 
900 ft.) across the Gulf Stream, made by means of the bathythermograph. This 
figure is a simplification of a section published by Iselin and Fuglister (1948, 
fig. 2). The simplification consists of drawing 5° intervals for the isotherms. 
Temperatures are in degrees Fahrenheit. The section extends from the con- 
tinental shelf at the left, across the Gulf Stream, where the 65° F. isotherm 
drops abruptly, well into the Sargasso Sea on the right. Even in this simplified 
drawing there is a great deal of fine detail. Just what is the cause of these fine 
variations we do not know, although internal waves are often mentioned as an 
explanation. The cold filament of water on the left-hand side of the Stream 
is clearly shown at a depth of between 200 and 300 ft., just left of center. 


MEANDERS AND EDDIES 


The position of the Gulf Stream is not always the same, nor is its path even 
approximately straight. Church (1937) was able to demonstrate the truth 
of the former statement conclusively on the basis of 1200 thermograph 
records from ships crossing the Stream. The number of crossings per week 
was about three or four, and hence it was impossible to develop any detail 
about the presence of wavelike disturbances along the Stream, or of eddies 
on either side. 

Fuglister and Worthington (1951) have prepared a chart showing the 
positions of the maximum cross-stream temperature gradients at a depth 
of 100 m., from all bathythermograph surveys made in the five years 1946 
through 1950 (see fig. 28 of the present study). These lines, of course, show 
a great deal more detail than Church’s. They confirm Church’s deduction 
that the position of the Stream varies from time to time, but, more 
important, they indicate that the Stream does not shift position bodily, 
but in wavelike patterns which have come to be spoken of as meanders. 
The first information on the way in which these meanders move was 
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obtained from the Multiple Ship Survey of 1950. This was the most de- 
tailed cruise ever made. Seven ships were employed. The time sequence of 
meander patterns was observed. 

According to the analysis by Fuglister and Worthington, two meanders 
in the western half of the surveyed region moved eastward at a rate of 
about 11 nautical miles a day. The water in the swiftest part of the Stream 
itself moves more than a hundred miles a day. The amplitude of the 
meanders nearly doubled in two weeks. Fig. 29 shows the mean tem- 


Fig. 28. Positions of the maximum cross-current temperature gradients at 
a depth of 100m., from all surveys in the period 1946-1950, according to 
Fuglister and Worthington (1951, fig. 4). 


perature, in degrees Fahrenheit, of the upper 200 m. layer of the western 
part of the survey area at the beginning of the period of observation. 

In the eastern half of the surveyed area a very much distorted meander 
was observed to break off into a clearly defined eddy. Fig. 30 shows the 
position of this eddy on June 17. The existence of eddies had been inferred 
before, but this one was very closely studied because so many ships were 
in the area. 

Fig. 31 shows a survey of the ‘edge’ of the Gulf Stream as determined 
from a single airplane flight, by means of the air-borne radiation thermo- 
meter (Stommel ef al., 1953). The black dots indicate the position of 
crossing of a strong surface temperature discontinuity as observed from 
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Fig. 30. A continuation of fig. 29 toward the east, showing the large eddy 
which developed on June 17, 1950 (Fuglister and Worthington, 1951, fig. 7). 
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Fig. 31. Edge of the Gulf Stream, as determined by air-borne radiation 
thermometer (see Stommel et al., 1953, fig. 1). The dots show the position of a 
strong horizontal temperature contrast. The curved-line segments indicate 
regions where there was clear visual evidence of the inshore edge of the Stream. 
The position of this area can be easily ascertained by reference to the index 
chart, fig. 18. 
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the airplane. The continuous curves indicate parts of the ‘edge’ where 
there was a visible indication of an edge: a sharp change in color, a change 
in number of whitecaps, and so forth (see frontispiece). The irregularities 
of the ‘edge’ in this survey are on a scale different from that of fig. 29.1 
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Fig. 32. Closely spaced measurements, by Worthington (19546, fig. 7), of 
surface current across the Gulf Stream. Surface velocities as computed in three 
ways are indicated: by geostrophic equation—solid black squares; by towed 
electrode, with correction applied—solid black circles; by successive Loran 
fixes—open circles. 


VELOCITY DETERMINATIONS ACROSS THE STREAM 


Worthington (19546) made three very closely spaced hydrographic 
traverses across the Stream in October and November, 1950. These are the 
best ever made, and for this reason it is worth considering one of them in 
some detail. Fig. 32 shows the surface currents across the Stream as deter- 
mined by the geostrophic equation, by towed electrodes with a correction 
factor applied, and by set of the ship as determined by successive Loran 


1 A more recent and detailed survey of the variability of the edge of the Gulf 
Stream is given by W. 8. von Arx, D. F. Bumpus, and W. 8. Richardson, in Deep-Sea 
Research, 3 (1955): 46-65, under the title ‘On the Fine Structure of the Gulf Stream 
Front’. 
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fixes. The outstanding features of this velocity profile are the narrowness of 
the current, the high velocities, and the countercurrent on the right of the 
Stream. The streaky, banded nature of the Stream, as shown by the curve 
connecting black squares, is probably not real. Neither the results obtained 
from the towed electrodes nor those from the Loran fixes show a similar 
streakiness of velocity across the Stream. It seems likely that the streaki- 
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Fig. 33. Closely spaced temperature measurements across the Stream, in 
degrees Centigrade. By Worthington (19540, fig. 4). 


ness is a false effect arising from internal inertial gravity or tidal waves 
acting upon the density field. 

Fig. 33 shows the thermal field across the same section of the Stream. 
Only the upper half of the ocean is shown, the total depth in that area 
being nearly 4000 m. The warm core is wider than the core of high velocity 
and extends toward the right into a weak countercurrent. This is contrary 
to the first natural conjecture that high downstream velocities and high 
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temperature are directly correlated. The main thermocline (10° C. iso- 
therm) drops from 200 to 900 m. across the Stream in less than 70 nautical 
miles. The temperature-versus-depth curve of the water below the 16° 
isotherm is nearly the same on the two sides of the Stream, except for a 
vertical displacement of about 700 m. The chief dissimilarity between the 
water above the main thermocline on the one side of the Stream and that 
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Fig. 34. The geostrophic current velocity (em./sec.), according to Worthing- 
ton’s (19546) closely spaced section (his fig. 10). 


on the other consists in differences of temperature of the upper 200 m. and 
the presence of a large body of nearly isothermal water (ca. 18° C.) between 
200 and 600 m. on the Sargasso Sea (right-hand) side. 

From a computation of transport, obtained from fifteen crossings, Iselin 
(1940) has shown that the total transport of the Gulf Stream above 2000 m. 
in this area is between 76 x 10° and 93 x 10® cm.8/sec. Fig. 34 shows the 
geostrophic velocity across the Stream in centimeters per second. The 
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current is sensible even at a depth of 1000 m. More recent crossings by 
Worthington, in 1954 and 1955, indicate that the current may be appre- 
ciable even at the bottom. 


VELOCITY SECTIONS BY DIRECT MEASUREMENT 


The only direct measurements of velocity at various depths thus far re- 
corded are those made on a six-day cruise of the Bear and the Caryn in 
July, 1952. The measurements were taken in a crest of a meander near Cape 
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Fig. 35. Bathypitotmeter measurements of current at two stations in the 
Gulf Stream. The solid black circles indicate velocity; the open circles, tem- 
perature. Both series of bathypitotmeter soundings, BP V and BP VI, were 
made in the left-hand edge of the warm core, in the part of the Stream having 
the highest velocity, with the ship hove to and drifting with the current, and 
were obtained by Dr Willem Malkus in June, 1952. 


Hatteras. Vertical velocity profiles made by the Malkus bathypitotmeter 
are shown in fig. 35. The northernmost station was made on the northern 
edge of the Stream; the other station shown was very close to the center of 
the Stream. The Watson propeller-type meter was also used and showed 
little variation of direction with depth. In general, there does not appear 
to be any great discrepancy between these velocities and the computed 
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geostrophic velocities, but the data are too irregularly spaced for quanti- 
tative comparison. Future survey work will, it is hoped, make possible 
such comparisons. 


PROBLEMS OF CONTOURING 


Problems of contouring are encountered in two ways: (i) in the con- 
struction of vertical temperature and salinity sections; and (ii) in the 


ae 


b 


Fig. 36,a and b. Different ways of contouring a property the distribution of 
which is measured only along the vertical dashed lines. See text. 


construction of charts showing horizontal distributions. There is never any 
difficulty in drawing isotherms in a vertical section, for example, when the 
temperature is a monotonic function of depth, but when the temperature 
passes through some maxima and minima, the interpretation is ambiguous, 
Fig. 36 illustrates how a single set of data can be interpreted by contouring 
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in two different ways. The essential decision to make in interpretation is 
whether there is an isolated central core or whether there is a tongue con- 
nected to the lower fluid. If the chart is regarded as a vertical temperature 
section the likeliest interpretation is that there is no connection and that 
the cold central mass is fresher; hence, there would be no violation of 
vertical stability. There is no proof, however, that the water cannot be 
momentarily unstable (from intense turbulence, the breaking of internal 
waves, etc.), and hence the alternative interpretation is not automatically 
ruled out. Simultaneous salinity measurements help to settle the issue, but 
they are not usually available in bathythermograph sections. Such tem- 
perature inversions are common in the side of the Stream toward the coast 
and have been discussed by Ford, Longard, and Banks (1952), who find 
them frequently made up of freshened water (30-5-34-5%, as opposed to 
36-5%, in the Sargasso water in the central North Atlantic). 

Similar ambiguities of contouring arise in drawing horizontal charts, and 
these ambiguities are likely to be even more perplexing because the lines 
of data on the chart are fewer. It is often impossible to decide whether to 
draw a stream with an eddy to one side or the other, a sharp S-shaped 
curve in a single stream, or three separate streams. The only successful and 
practical way to resolve this difficulty is to meet it at the time when it 
arises, by proper planning on board the survey vessel. This is a good 
example to demonstrate why cruises cannot be entirely planned ahead of 
time and why the uncomfortable job of taking oceanographic data cannot 
be left entirely to untrained personnel. 


THE STRATEGY OF EXPLORATION 


In the present state of physical oceanography most of the effort must be 
devoted to field work, in an attempt to explore and describe the chief 
physical features of oceanic phenomena such as the Gulf Stream. This does 
not mean that theoretical studies are not worth while, but it does mean 
that the chief features are so poorly known that a serious theoretical 
analysis of the Gulf Stream cannot be made at present. We have seen how 
the oceanographic surveys of the 1930’s have given a broad, general view 
of the ocean surface circulation and have delineated the major, average 
thermal and salinity structure of the waters. These early surveys are all 
characterized by a common strategic practice: the courses and spacing of 
stations were planned in advance; and the sections actually made corres- 
ponded as closely to these preliminary plans as the exigencies of current 
set, weather, and instrumental failure would permit. The early Meteor 
studies of the South Atlantic, the present-day West Coast survey, and the 
first ten years of Atlantis cruises in the Gulf Stream, are all examples of 
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this type of strategy. So far as the object of exploration is a broad, coarse, 
climatological-mean picture, these premeditated cruises are successful. 

The bathythermograph introduced a new strategy into the game. Its 
great value lies perhaps less in the continuous nature of its trace with 
depth than in the fact that it can be used from a ship under way and that 
the data obtained from it are quickly and easily interpreted. The oceano- 
grapher studying the Gulf Stream nowadays keeps a running plot of all 
bathythermograph information as it comes in. On the basis of his inter- 
pretation of this information he frequently alters the course of the ship 
to explore whatever feature he is particularly interested in. The result is 
that he directs the ship in a series of traverses which he hopes will cross 
the Stream again and again, in order to obtain a three-dimensional image 
of its thermal structure and currents. This new strategy has endless 
possibilities and ramifications, only a few of which have been explored. 
For example, the chief tactic has been to follow the Stream by a series of 
zigzag legs. Fig. 28 shows the positions of the Stream as determined by 
such zigzag bathythermograph cruises during the period 1946-1950, ac- 
cording to Fuglister and Worthington (1951, fig. 4). 

More than a single ship may be employed to advantage in this type of 
operation: the Multiple Ship Survey of 1950 employed the zigzag tactic 
and was able to delineate for the first time a meander pattern showing 
several distinct waves and an eddy breaking off. The expense and effort 
involved in portraying this picture were considerable. 

Another type of strategy was developed during the short 1952 cruise of 
the Bear and Caryn; one ship was employed to make the rapid bathy- 
thermograph traverses of the Stream, and the other was used in making 
the slow hydrographic stations and velocity measurements. In this way 
it was possible to keep a constant check on the position of the slower vessel 
relative to the chief thermal features of the Stream. This strategy requires 
at least two ships, one of which should be capable of reasonable speed; 
they should never be farther apart than a few hundred miles, otherwise the 
cruise loses the effectiveness of the multiple-ship type of operation. Even 
two ships face a formidable situation if the ocean feature under obser- 
vation is changing rapidly. The 1952 Bear—Caryn cruise is an example. 
The evidence of six crossings of the Stream showed that the ships were 
working on the crest of a meander in which the deep velocity and thermal 
fields were essentially stationary during the six days spent in the area. The 
surface salinities, as indicated by a recording conductometric cell, were 
apparently changing rapidly with time, in the course of a widespread in- 
vasion of freshened water from coastal areas (Chesapeake Bay, Delaware 
Bay, and Hudson River) which was spreading along and across the Stream. 
The surface salinities plotted in the course of the cruise cannot, therefore, 
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be interpreted as an instantaneous picture at any time. It is possible, of 
course, to reduce all observations to their approximate positions at any 
epoch, if the surface velocities are well known and steady, but the practical 
difficulties are formidable and the final distribution is questionable. Thus 
the two ships were able to delineate the quasi-stationary thermal and 
velocity fields fairly well, but were completely unable to do the same for 
the rapidly changing surface salinity field. It will be very interesting to 
see what new strategic plan will meet such situations. 

There are a number of possibilities. For example, more use could be 
made of drifting and anchored recording and radiotelemetering instru- 
ments. Instruments left on the bottom for long periods might also be used. 
Electrical potential gradient and pressure seem to be the best variables to 
measure on the bottom, where fluctuations in temperature and velocity 
may be insignificant. The airplane can be used in measuring surface tem- 
perature rapidly over large areas. It is too early to foresee the new tactics 
these will demand. 

A leg crossing the Gulf Stream may be estimated as 80-100 nautical 
miles in length. To take twelve hydrographic stations, each of two casts,” 
requires about 70 hr., weather permitting. A leg with a similar number of 
bathypitotmeter velocity measurements might consume 56 hr. If both 
hydrographic stations and velocity stations are made, the total time in 
crossing is 100 hr. These periods may be contrasted to 10 hr. when the 
bathythermograph alone is used, and 12 hr. for measuring with towed 
electrodes. An airplane can traverse the same leg in less than 40 min., but 
the data it can obtain are pretty well limited to surface temperature and 
photographs. 


TRANSFER PROCESSES ACROSS THE STREAM 


In 1936 Rossby (1936b) showed that if the Gulf Stream were regarded as 
a wake stream in the sense used by Tollmien (1926) the observed down- 
stream increase of mass transport could be accounted for. Rossby also 
suggested that large lateral exchange processes might be at work carrying 
Sargasso Sea water across the Stream and into the slope water on the 
left-hand side. It now appears that the increase of transport in the down- 
stream direction can be accounted for by other means (see Chapters VII 
and VIII), and that the turbulent-jet analogy is not necessary; but the 
question still remains: what kind of transfer occurs across the Stream? 


2 Usually, at any deep-sea hydrographic station not more than twelve reversing 
bottles and twelve pairs of reversing thermometers are attached to the cable. Three 
or four separate lowerings must therefore be made, in order to obtain samples at 
all depths. Each lowering is called a ‘cast’. 
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Indeed, the T-S diagrams (see fig. 37) of water masses below 500 m. on 
opposite sides of the Stream are so similar that it seems very likely that 
in deep water such a transfer occurs. The large eddies which detach from 
the Stream, as exemplified by the eddy observed on the Multiple Ship 
cruise of June, 1950, certainly effect a transfer of water across the Stream. 
It should be noted that they are rather solitary phenomena in the Gulf 
Stream, at least so far as we now know. For example, in June, 1950, there 
was only one eddy in the process of detachment in a distance of more than 
1200 miles along the Stream. The wavelike meanders may transfer momen- 


---- SLOPE WATER 
—— CENTRAL ATLANTIC WATER 


Fig. 37. Comparison between the mean temperature-salinity correlation 
curve for the slope water and that for the Sargasso water, on the Atlantis 
sections from Chesapeake Bay to Bermuda. From Iselin (1936, fig. 23). 


tum, but not heat and salt. A transfer of properties across the Stream on 
a scale smaller than that of these large eddies does not appear to be 
indicated. Ford, Longard, and Banks (1952) have called attention to a 
narrow band of fresh cold water extending from the surface to a depth of 
about 400 ft., and located within the left-hand side of the Stream. Fig. 38 
shows a bathythermograph section through such a filament of cold fresh 
water. Ford et al. point out that of the ninety-seven crossings of the 
Stream made in June, 1950, fifty-seven show the cold layer. It is quite 
possible that on the other crossings the bathythermograph lowerings 
missed the cold layer altogether—it appears to be a very slender ribbon, 
not more than about 5 miles wide. Fig. 39, from Ford et al., shows the 
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Fig. 38. Salinity, temperature, and velocity across the Gulf Stream (above), 
showing the fresh filament on the inshore edge; and (below) a bathythermo- 
graphic section, in degrees Fahrenheit, based on data collected at the same 
time. From Ford, Longard, and Banks (1952, fig. 3). 


Fig. 39. Chart showing the filament of cold water along the edge of the 
Stream. The contour lines represent surface temperatures obtained by a con- 
tinuously recording thermograph on board H.M.C.S. New Liskeard. The ship’s 
course is shown by arrowed lines. The current parallels the isotherms. The 
contours are drawn in by eye, and, of course, are subject to the ambiguities of 
contouring (see p. 59). From Ford, Longard, and Banks (1952, fig. 2). 
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position of the cold water in June, 1950. Fig. 40, from the same article, 
shows the positions of the temperature inversions for the entire Multiple 
Ship operation, by shaded areas. It is quite evident from the temperature 
and salinity of this water that it does not come from depth, but must 
originate from somewhere along the shelf near Cape Hatteras. Were this 
filament of fresh water a permanent feature, the supply of fresh water 
required to maintain it would be of the order of magnitude 104 m.3/sec., 


Fig. 40. Positions of temperature inversions along the inshore side of the 
Gulf Stream, determined on the 1950 Multiple Ship Survey and shown by dark 
areas. The areas of shelf and slope are indicated. The Gulf Stream current is 
shown by arrows between the curved solid line and the dotted line, which serve 
to mark the left- and right-hand sides of the Stream. From Ford, Longard, 
and Banks (1952, fig. 5). 


which could be supplied by river discharge along the coast. The very fact 
that such a slender filament can preserve its integrity along at least 
1200 miles of the Gulf Stream is an indication that smali-scale turbulent 
processes tending to transfer properties across the Stream in the upper 
layer are inconsiderable. 
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SOURCES OF THE GULF STREAM WATER MASSES 


As was explained in Chapter IV, the water which flows out of the Florida 
Straits originally comes in large part from the southern half of the North 
Equatorial Current, and in part from a branch of the South Equatorial 
Current which is apparently split off from the South Atlantic by the 
peculiarly wedgelike coastline of Brazil. This water flows through the 
Caribbean, and then, without mixing with the waters endemic to the Gulf 
of Mexico, emerges from the Florida Straits in very much the same state 
as when it entered the Caribbean. Because of the large admixture of 
Antarctic Intermediate Water at mid-depths which this water has acquired 
from the South Atlantic, there is a distinct salinity minimum (between 600 
and 800 m. depth) in the water coming out of the Florida Straits. This 
salinity minimum is not so marked in the water flowing westward north of 
the West Indies. Therefore, the difference in the intensity of the salinity 
minimum serves as an indicator, or tracer, for distinguishing water in the 
Stream which has come through the Caribbean from that which has joined 
the Stream north of the Indies. Iselin (1936) has used the intensity of the 
salinity minimum in this way to trace the source regions of various 
portions of the Gulf Stream in various areas of the North Atlantic. He has 
plotted charts of salinity anomalies which give a rough indication of the 
sources and degree of mixing of mid-depth waters. 

In a detailed analysis of a Gulf Stream section made in April, 1932, off 
Chesapeake Bay, Iselin (1936) shows that almost all the water colder than 
8° C. must have joined the Stream north of the Florida Straits, an amount 
equal to 15 per cent of the total transport of 82 x 10° m.3/sec. which he 
computed for the Stream using a 2000 m. reference level. This is consistent 
with the fact that most of the water at a temperature below 8° C. is blocked 
from flowing through the Florida Straits by the shallowness of the channel. 
The water between 8 and 20° C. flowing in the Gulf Stream off Chesapeake 
Bay is made up of nearly equal parts of water which has come through the 
Florida Straits and of water from the North Atlantic north of the Indies; 
together these components amount to about 71 per cent of the total 
transport. An additional 13 per cent of the total transport occurs in the 
warm core of water consisting of surface waters warmer than 20° C., but 
these waters are so subject to wind mixing and to atmospheric cooling 
that for them the T-S relationship means very little as a tracer. Less than 
1 per cent of the total Gulf Stream transport off Chesapeake Bay involves 
entrained masses of coastal, shelf, and slope waters on the inshore edge of 
the Stream. The reader is reminded that these figures are based on an 
assumed reference level. 

The T-S method actually gives us a rather broad picture of the possible 
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origins of various parts of the Stream, but in order to refine the method 
it would be necessary to have some quantitative information about mixing, 
and it would be necessary to make more widespread and frequent hydro- 
graphic sections than are now available. This would require great expense 
and effort with present techniques. In addition, the uncertainty about the 
choice of the reference level makes transport calculations rather fuzzy. 
Thus, in a few words, sparsity of data, ignorance of mixing, and uncertain 
transport computations conspire to maintain the technique of water-mass 
analysis at a qualitative rather than a quantitative level. 

Although the hope that we shall ever have really frequent and wide- 
spread hydrographic data is very dim, there is a need for repeating certain 
hydrographic stations made in the past. With the waning popularity of the 
hydrographic station there has been a tendency toward slovenliness in the 
work which makes some of the new data worthless for water-mass analysis. 
In order to offset this tendency, as well as to gather new precise data on 
possible secular changes in the properties of the deep waters, Worthington 
began, in 1954, to rerun all the old Atlantis sections of twenty-odd years 
ago. This tedious and very demanding work, requiring precision measure- 
ments under difficult environmental conditions, is now half finished. When 
completed, it will be a useful step toward a firm and secure beginning for 
a long-term study of slow climatic changes of the ocean. 

Just as Iselin (1936) used the salinity minimum at the depth of 700 m. 
to trace water masses in the Gulf Stream, Richards and Redfield (1955) 
have recently used dissolved-oxygen deficiency to attempt to trace water 
at about 200 m. It is known that the waters from the Straits of Florida 
(at 200m.) are more deficient in dissolved oxygen than are the waters 
that join the Stream north of the Indies. A diagram analogous to the T-S 
diagram is used, but the codrdinates are dissolved-oxygen concentration 
and o;. The results of Richards and Redfield’s analysis of eight available 
Gulf Stream cross sections made within the period 1950-1953 were highly 
variable; that is, the amount and position of Florida Straits water at 
200 m. present in the Gulf Stream varied greatly from section to section. 
Of course, eight sections are too few from which to establish the nature of 
the variability and to relate it to any continuously measured quantity 
such as variation in difference of sea level across the Straits; but it is 
certainly worth while to demonstrate the fact of variability itself. The 
study also showed that there are separate filaments of the Florida Straits 
type of water in the countercurrent to the right of the current. Possibly 
this has some bearing on the question whether the countercurrent is 
frictionally (p. 97) or advectively (p. 123) produced. 
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THE FLORIDA CURRENT 


Most of the water in the Florida Current comes through the Yucatan 
Channel from the Caribbean Sea, rather than from the Gulf of Mexico 
(Iselin, 1936). The water is forced through the long channel between the 
Florida peninsula on one side and the islands of Cuba and the Bahamas 
on the other by a head of water of about 19 cm., as inferred from a leveling 
survey across Florida. The Straits between Key West and Havana are 
140 km. wide, and their greatest depth is 1500 m. This channel becomes 
narrower and shallower downstream, takes an abrupt 90° turn to the left, 
and reaches a minimum cross section off Miami, where the width is 80 km. 
and the greatest depth is 800 m. The stream then flows north along the 
continental shelf to about 33° N. latitude, where it leaves the shelf and 
flows into deep water south of Cape Hatteras. From Cape Canaveral 
north, the Florida Current increases in mass transport. This is particularly 
obvious at the place where it leaves the continental shelf and begins to 
carry water at temperatures lower than 8° C. along with it. There is also 
supposed to be an increase of volume transport by confluence with the 
Antilles Current immediately north of the Bahamas (Iselin, 1936), but 
evidence concerning the transport of the Antilles Current is conflicting. 

The mass transport of the Florida Current is estimated to be about 
26 x 10° m.3/sec. The computed values vary considerably, as is shown by 
computations (Montgomery, 1941b) made from four Aflantis sections 
across the Straits at Havana: March 4, 1934; February 19-20, 1935; 
April 12-13, 1935; and March 24-25, 1938. The transports computed are 
26-0, 30-3, 29-0, and 26-0 x 10° m.3/sec., respectively. There are numerous 
sources of error in such computations: first, the level of vanishing hori- 
zontal pressure gradient must be assumed; secondly, the geostrophic 
equilibrium may not hold strictly; thirdly, there are gaps between the end 
stations and the shore where transport must be estimated; and, fourthly, 
tidal effects on the density field may upset the computation. Woodcock 
(see Parr, 19376) made a series of anchor stations across the Straits at 
Miami which show marked semidiurnal variation in the density structure. 
Parr (1937) has described a highly speculative process of cross-stream flow 
based on Rossby’s (19366) wake-stream analogy. It seems to me that the 
tidal influence on the density structure in the Straits of Florida might be 
of the nature of an internal seiche. The dimensions of the channel at 
Miami and the density’ structure appear to favor a resonance in a cross- 
stream semidiurnal internal seiche. 

Since the tides in the Gulf of Mexico are very small, one would suppose 
that there is a tendency for a semidiurnal progressive tidal wave to move 
upstream from Miami. A study of the tidal constants at various points 


HYDROGRAPHY OF THE STREAM 69 


along the Straits of Florida suggests that this is in fact true, the wave 
being of the Kelvin type (Thomson, 1871), and the range of tide being 
greater along the United States coast than along the coast of Cuba or the 
Bahamas. 

Electromagnetic measurement.—An interesting series of observations is 
being accumulated through the codperation of the Western Union Tele- 
graph Company and the Woods Hole Oceanographic Institution. Elec- 
trodes have been buried in the beaches at Key West and Havana, and 
connected by a submarine cable. The potential developed by the water 
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Fig. 41. Mass transport of the Florida Current, as determined from electrical 
measurements made on the Western Union Telegraph Company cable between 
Key West and Havana by P. J. Moore and analyzed by Wertheim. 


moving through the earth’s magnetic field is about 1 volt across the Straits, 
and, under the assumption of a nonconducting bottom, indicates a volume 
transport of approximately 26 x 10® m.3/sec. The usefulness and feasibility 
of this kind of measurement will appear after more years of records have 
been obtained. 

Data for about two years have now been obtained from the submarine 
cable between Key West and Havana, across the Florida Current. 
Electrical potential measurements covering several more years will prob- 
ably be required before any definite conclusions can be drawn concerning 
fluctuations in the actual volume transport. Wertheim (1954) has ex- 
amined the first year’s readings and has plotted the transport as a function 
of time (fig. 41). Each point is a 24 hr. average of readings of potential 
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(thus, tidal variations and magnetic storm signals are averaged out). 
Perhaps the most striking feature of these fluctuations is the extreme 
rapidity with which major changes in transport can occur. 


THE NORTH ATLANTIC CURRENT 


There are not a great many observations of the North Atlantic Current 
(Defant, 1939; Neumann, 1940; Soule, 1950). The Atlantis made a section 
from 37° N. to 52° N. along 30° W. longitude in 1931 (Iselin, 1936), here 
shown as fig. 42, and the International Gulf Stream Expedition of 1938 
(the Altair and the Armauer Hansen) made some sections in the area. What 
data there are suggest that the North Atlantic Current is not a single, 
narrow current like the Gulf Stream off the United States coast, but is 
made up of several distinct broad currents. Just how permanent these 
currents are has not been determined. Iselin (1936) and Sverdrup, Johnson, 
and Fleming (1942) have supposed that the Gulf Stream System splits into 
a number of branches just east of the tail of the Grand Banks. Fuglister 
(19516) has proposed that there is no real branching at all. He has drawn 
a schematic set of streamlines which fit the data well, and calls his schema 
the Multiple Current Hypothesis. He suggests that an instantaneous chart 
would show not a continuous stream, but a number of disconnected fila- 
ments of current. This pattern changes from time to time. The possibility 
that the current system of the North Atlantic Current is irregular and 
varying and actually discontinuous is very disconcerting. It is even more 
disturbing to find that Fuglister is able to draw these alternative inter- 
pretations even about the water as far west as 65° W. longitude, where we 
have always thought of the Gulf Stream as unambiguously identifiable. 
Figs. 43-45 represent three interpretations of the nature of the north- 
eastern parts of the Gulf Stream, according to Fuglister (1955, charts 
3a—3c, all based on the same data). These illustrate clearly the difficulties 
of adequate description of the Gulf Stream System. In my opinion, the 
interpretation in fig. 45 is a bit forced, from an attempt to spread iso- 
therms apart wherever possible. On the other hand, fig. 43 is forced in 
the opposite way. Fig. 44 corresponds most nearly to Fuglister’s picture of 
multiple streams and seems intuitively, to me, to be the most natural 
interpretation. 
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Chapter Six 


THE WIND SYSTEM OVER THE 
NORTH ATLANTIC 


It is important to describe briefly the winds over the North Atlantic 
Ocean, because it is generally supposed that the surface circulation of the 
ocean is caused by them. Many of the older textbooks of meteorology 
emphasize the mean wind distribution so heavily that one is likely to over- 
look the considerable fluctuations which occur in wind distribution from 
day to day. At present, theories such as those described in Chapter VII 
take into account only the mean wind distribution. Just what the effects 
of the widespread interruptions and irregularities of atmospheric flow on 
the ocean circulation are, is still to be told. 


NORMAL CIRCULATION 


Chase (MS, 1951) has made an extensive study of the surface pressures in 
the last half-century on the North Atlantic, using the historical weather 
maps for the Northern Hemisphere (U.S. Weather Bureau, for the period 
beginning in 1899), and has since incorporated information from the normal 
weather maps, Northern Hemisphere (U.S. Weather Bureau, 1952). The 
normal sea-level pressure (yearly average) is shown in fig. 46. By far the 
largest area of the ocean is covered by a single high-pressure cell with winds 
traveling in a clockwise direction about its center. The normal sea-level 
pressures at two different seasons of the year are shown in figs. 47 and 48. 

The center of the high is at its most northeasterly position in January ; 
by March it has moved almost 1200 miles to the southwest and is at its 
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Fig. 46. Map showing annual mean sea-lev i o Joseph 
Chase (1951, MS). Pressures are giv 
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Fig. 47. Map showing mean sea-level pressure for January (Chase, 1951, 
MS). Pressures are given in the excess of millibars over 1000 mb. 
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most southerly position. During the spring the center of the high gradient 
moves due north, and from July on, due east, until in November the 
center lies almost exactly at the center of the yearly average chart (Chase, 
1951). Chase described the changes in pressure distribution as follows 
(his p. 4): 


The general appearance of the mean monthly pressure distri- 
butions is similar to that of the annual average. Although the 


ca 


Fig. 48. Map showing mean sea-level pressure for July (Chase, 1951, MS). 
Pressures are given in the excess of millibars over 1000 mb. 


Westerlies are normally stronger in winter than in summer, the 
circulation around the Bermuda—Azores High, when taken as a 
whole, is greater in summer. The mean maps for January and 
July...illustrate this fact. In July the Trades are stronger in 
proportion to the increased pressure gradient and the east and 
west ends of the high are much better developed. The Westerlies 
show a decrease in width of band and in pressure gradient. The 
decrease is greatest in the area north of 50° N. South of that 
parallel, the average Westerlies are nearly as strong as in winter. 
This is the result of more steadiness of direction in summer rather 
than of high velocities. 
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INTERRUPTIONS IN THE NORMAL CIRCULATION 


Chase (1951, p. 6) also pointed out that only about half of the daily surface- 
pressure maps look at all like the monthly normal maps: ‘A large pro- 
portion of the daily maps resemble the average maps in position and 
orientation of the Bermuda—Azores High and in the circulation about it. 
These situations are most common in the summer (frequency estimated at 
some 60% of the time, and least common in winter [about 30%]). The 
frequency for fall is slightly higher than spring.’ The rest of the daily 
maps are quite irregular, showing a number of interruptions which he has 
classified as: fronts; stagnant lows; linkages of the Bermuda—Azores High 
to continental highs; and hurricanes. 

Some fronts pass to the north of the Bermuda—Azores High without 
particularly disturbing it. Others, such as the front shown in fig. 49, seem 
to force it toward the south. Normally these frontal interruptions occur 
about twice as often in the winter as in the summer, the duration of each 
interruption being about a week. 

The stagnant lows generally persist over the North Atlantic for a 
longer time. Such lows occur about half as often as the frontal 
passages. 

Frequently a continental high coalesces with the Bermuda—Azores high, 
causing interruptions of the normal circulation, and these in turn result in 
widespread irregularities and shifts in wind. 

Hurricanes occur in summer and autumn, but primarily in September. 
They are by far the least frequent interruption to the mean circulation. 
Fig. 50 represents a very rare situation, that of four hurricanes occurring 
simultaneously. 


THE WIND STRESS ON THE OCEAN 


Even though we have a good idea of average wind conditions over the 
ocean, we can only make crude approximations to the stress produced by 
these winds on the ocean surface (Rossby, 19365). 

The first approximate values of wind stress were computed (Ekman, 
1905) from the mean slope of the water surface in the Baltic Sea, as 
observed by tide gauges. According to a first-order theory, the local slope 
of a water surface depends only upon three factors, namely, the surface 
stress, the bottom stress, and the depth of the water (if the water is 
assumed to be vertically homogeneous). The vertical internal turbulent- 
eddy viscosity does not appear in this expression. In a natural body of 
water where neither the depth nor the wind distribution is uniform, and 
very little is known concerning the bottom stress, this method is very 
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Fig. 49. A daily map showing the disturbance of the Bermuda—Azores high 
caused by a front pushing in from the northwest (Chase, 1951, MS). 
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Fig. 50. A rare case of four simultaneous hurricanes over the North Atlantic 
(Chase, 1951, MS). 
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uncertain. Recently (Van Dorn, 1953), very precise measurements have 
been made of the slope of the water surface in a large regular artificial 
basin. Under such ideal conditions the method gives better results, but 
the fetch is so short that large waves do not develop. 

Experiments to determine the drag of wind on a free water surface have 
also been made by Keulegan (1951) and Francis (1951) in wind tunnels. 
There is evidence that the presence of large waves may not have an 
important effect on the drag coefficient. Francis (ibid.) suggests that only 
the very smallest wavelets and ripples act as roughness elements. 

Another method of determining the wind stress stems from laboratory 
results in the field of aerodynamics. In these experiments it has been 
found that equations for the stress in a boundary layer can be computed 
from a knowledge of the variation of wind velocity with altitude over the 
first 10 m. above the sea. Rossby and Montgomery (1935) applied these 
aerodynamic equations to observed wind profiles. Table 1, adapted from 
Sverdrup, Johnson, and Fleming (1942, p. 67), gives drag computed from 
these equations. 


TABLE 1 


CORRESPONDENCE OF WIND STRESS (IN DYNES PER SQUARE CENTIMETER) 
To WIND VELOCITIES MEASURED FIFTEEN METERS ABOVE SEA LEVEL 


Wind velocity (m./sec.).... 2 4 6 8 10 12 14 16 18 
Wind stress (dynes/em.?) .. 0:04 0-16 0-34 1-81 2-83 4:09 5-56 7-25 9-20 


Source: Sverdrup, Johnson, and Fleming (1942, table 67). 


No account of the thermal instability or stability of the air relative to 
the water is included in these studies. There is need for ingenious and 
careful experiments and observations on this problem. 

A measure of the wind stress on the sea can be obtained (Sheppard 
and Omar, 1952) from purely meteorological data, in the following 
way: 

Let w and v be the horizontal components of the wind along axes x and y, 
where « is in the direction of the surface wind uw). One then supposes the 
following equilibrium to hold in the x-direction: 


v—-—+-—=0. (1) 
Let us suppose that v, is the y-component of the geostrophic wind; then 


OT, ; 
Be Pits) ? (2) 


WIND SYSTEM OVER THE NorTH ATLANTIC 81 


and we now integrate this equation from the surface z=0 to a height 
z=h, where T,=0; hence 79, the stress at the surface, is given by 


h 
=I, P(¥g—2) dz. (3) 


The method has been used successfully to measure wind stress over the 
land by a number of investigators, and over tropical seas by Sheppard and 


10> X SHEAR-STRESS COEFFICIENT 
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WIND SPEED (m./sec.) AT 10 m. ABOVE WATER LEVEL 


Fig. 51. Plot showing correlation of previous field values of shear stress 
with model results. X=wind-tunnel results; S=Shoulejkin; P= Palmén; 
E=Ekman; M=Montgomery; R=Roll; B=Bruch; H=Hela; J=N. K. 
Johnson; D=Durst; C=Corkan; L=Sutcliffe. Redrawn from Francis (1951, 
fig. 9). 


Omar (1952). Sheppard, Charnock, and Francis (1952) have shown, from 
a series of observations at the Scilly Islands, that this method does not 
work well in the westerlies. 


Francis (1951) has given an interesting summary of the determination 
6 SGSs 
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of the drag coefficient y? (y?=79/p,u*, where p, is the density of air, and u 
is the wind velocity at 6 m.) by these different methods. To date, values of 
the drag coefficient range from 2 x 10-* to 10-*. The largest of these were 
determined in strong winds by the slope-of-the-sea-surface method, the 
smallest values under light winds from shipboard. Fig. 51 shows various 
values of the drag coefficient. In spite of the intensive effort that has been 
made to measure the stress of the wind on the sea, there is so much scatter 
in the various determinations that the stress is still not well known. It is 
hard to imagine an object of study more important to the physical oceano- 
graphy of ocean currents (see Montgomery, 19360). 

Finally, mention should be made of a fundamental difficulty in pre- 
paring charts of the mean wind stress over the ocean: because of the non- 
linear relation between surface wind and surface wind stress, the use of 
mean wind charts for the computation of mean stress charts is not satis- 
factory, especially in temperate and polar latitudes, where there are major 
fluctuations in the observed wind field every few days. The proper way to 
produce a mean stress chart would be the laborious one of first constructing 
daily stress charts, and then obtaining vectorial mean stress from them. 
Actually, the present state of the theory of wind-driven currents does not 
justify such refinement. An approximate method, proposed by Reid 
(1948), is based upon the use of charts of mean wind plus charts of mean 
variability. 


Chapter Seven 


LINEAR THEORIES OF 
THE GULF STREAM 


Perhaps the most striking feature of the large-scale horizontal surface 
circulation in the North Atlantic Ocean is its east-west asymmetry. 
Although the mean winds are very broad and diffuse over the entire ocean, 
the currents along the western shores of the Atlantic are very narrow and 
intense. This same westward intensification of surface currents is observed 
in other oceans, and has been the object of study by the writer (Stommel, 
1948), who proposed that the latitudinal variation of the Coriolis parameter 
was the cause of the asymmetry, and by Munk (1950), who first developed 
a thorough dynamical theory which yields the main features of the ocean 
circulation and the correct order of magnitude of the transports of ocean 
currents from the computed wind-stress distribution. 

The Kuroshio is the counterpart of the Gulf Stream in the North Pacific. 
In the Indian Ocean the Agulhas Current hugs the coast of Africa. In the 
South Atlantic there is the Brazil Current. The coastlines and areas of these 
oceans are quite different, so we are tempted to think that local bathy- 
graphic peculiarities are not an essential influence in the western intensifi- 
cation of ocean currents; thus the existence of the Straits of Florida is not 
essential to the formation of the Gulf Stream. Ifthe Antilles were excavated 
the Gulf Stream would still exist. However, the South Pacific Ocean offers a 
somewhat embarrassing exception. There does not appear to be any current 
of great intensity off Australia; in fact, the Humboldt Current off Peru is the 
strongest South Pacific current, and it lies in the eastern part of the South 
Pacific. With this important exception the following rule does seem to hold: 

6-2 
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On the western edge of most of the world’s oceans there is a system of strong 
currents. Now, we may ask ourselves, why is this so? If the wind system, a 
broad, widespread phenomenon, drives the ocean currents, how is it that 
the resulting ocean-current system should be so asymmetrical; in particular, 
why should the strongest currents be squeezed into a narrow belt on the 
western edge of the oceans? It is this question which the theories of Stommel 
(1948), Munk (1950), and Hidaka (19496) attempt to answer. 
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Fig. 52. Schematic wind system (broken lines) and currents (solid lines) 
which would occur were there no asymmetry in the circulation. 


THE THEORY OF WIND-DRIVEN OCEAN CURRENTS 


The physical situation is as follows. Suppose we consider a large ocean 
basin, such as is depicted in fig. 52, with an anticyclonic wind circulating 
over it. In both figs. 52 and 53 the streamlines of ocean current are shown as 
solid lines; those of the wind, as broken lines. Our intuition tells us that this 
wind will produce an ocean surface current in an anticyclonic sense. From 
what little we know about the distribution of properties in the deeper parts 
of the ocean it seems reasonable to suppose that the current which the 
wind induces does not extend to the bottom. Our intuition also tells us that 
if the ocean-current system has arrived at a state of steady motion under the 
stress of the wind, the free surface and isopleths of density in the ocean will 
have adjusted themselves in such a way as to produce the horizontal 
pressure gradients necessary to balance the Coriolis forces acting upon the 
moving water. But that is about as far as our intuition takes us, and it 
gives no hint of the necessary asymmetry in the current system. To pro- 
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gress further it is necessary to consider the processes tending to increase or 
decrease the vorticity of the ocean water. First we shall consider these pro- 
cesses qualitatively, and then proceed to a mathematical demonstration. 
The vorticity of a vertical column of water! is taken as positive if counter- 
clockwise, and negative if clockwise. If the spin is measured in its relation 
to the earth, we speak of relative vorticity; by adding the Coriolis parameter 
to the relative vorticity we obtain the absolute vorticity of the column. 
Unlike the velocity distribution, the vorticity distribution in the ocean may 


Fig. 53. Schematic relation of wind (broken lines) and current (solid lines), 
with asymmetry. 


be discussed without reference to the pressure field—a fact which greatly 
simplifies consideration of the qualitative factors. (Formally, the vorticity 
equation is obtained by elimination of the pressure from the equations of 
motion by cross-differentiation.) The wind system in fig. 52 is one that 
would tend to decrease (make negative) the relative vorticity of all the water 
columns in the ocean. In the steady state of motion which must exist under 
a prolonged exposure of the ocean to this wind system, the relative vorticity 
has a fixed value at each point of the ocean independent of time. This 


1 The vorticity of an element of fluid is twice its angular velocity. An element of 
fluid in a shearing motion is subjected to an instantaneous spin. The qualitative 
reasoning is done entirely in terms of the total vorticity of a vertical water column— 
that is, the sum of the vertical components of vorticity of all the elements of a 
vertical water column. The mathematical definitions of several different kinds of 
vorticity are given in Chapter VIII, equations (3) and (5). 
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means that other processes must act to counteract the negative-vorticity 
tendency due to the wind alone. 

The linear theories involve two processes, the first of which is friction. 
Since, as we have already seen, the ocean surface circulation does not seem 
to be frictionally bound to the bottom, we inquire whether it is frictionally 
bound to the ocean shores by horizontal eddies. Such a horizontal viscosity 
would provide a positive-vorticity tendency over the ocean shown in fig. 52. 
A numerical check, made using (i) values of lateral-eddy viscosity inferred 
from the distribution of conservative properties, and (ii) a horizontal 
oceanic circulation which looks like the wind system, without any evidence 
of asymmetry, requires a circulation many times as fast as the real ocean 
circulation for lateral friction to produce a positive-vorticity tendency 
strong enough to balance the wind-stress-vorticity tendency. 

The second process is the tendency of planetary vorticity. Regardless of 
hemisphere, columns of water moving northward without convergence or 
divergence have a negative-vorticity tendency, and those moving south- 
ward have a positive-vorticity tendency. This follows from the conservation 
of angular momentum, or, to put it in other words, the variation of the 
Coriolis parameter with latitude. Since the net meridional transport of 
water across a parallel of latitude is zero (as much water moves north as 
south), the planetary-vorticity tendency is positive for water in the eastern 
part of the ocean presented in fig. 52, and negative for water in the western 
part. Therefore the planetary-vorticity tendency alone is incapable of 
balancing the wind-stress-vorticity tendency. In the steady state we must 
have azero over-all vorticity tendency, by definition. That is, at every point 
in the ocean the wind-stress-, frictional-, and planetary-vorticity tendencies 
must cancel out one another. 

The distribution of the wind-stress-vorticity tendency may be regarded 
as fixed, let us say of an order of magnitude —1. If there were only a broad 
current system without the asymmetry that is actually observed but with 
a transport of water similar to the observed transport, the frictional- 
vorticity tendency would be of a smaller order of magnitude, say +0-1, and 
the planetary-vorticity tendency would be of the order —1 in the western 
and +1 in the eastern part of the ocean shown in fig. 52. Thus there would 
be an approximate balance of tendencies in the eastern part of the ocean, but 
the western part would not be in equilibrium; hence, for such an ocean as 
this in a steady state of motion, a symmetrical current system is not 
physically possible. The state of affairs for a symmetrical circulation is 
summarized in table 2. 

If we let the current system be strongly asymmetrical, as it is in fig. 58, 
we do not seriously affect the balance between the wind-stress-vorticity 
tendency and the planetary-vorticity tendency in the eastern part of the 
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TABLE 2 


Vorticiry TENDENCIES IN A SYMMETRICAL CIRCULATION 


North-flowing South-flowing 
Vorticity tendency currents in western currents in eastern 
side of ocean side of ocean 
VV IEICE SPN OSS Soy eue is).o1 5 jo-cuevrae, woe els: one's — 1:0 — 1:0 
HEC ELOMA Nets ar Ae eos, sale ois eves sere + 0-1 + 0-1 
PIGmGiCiny pas eboo oH Oem onl OOo — 1:0 + 1-0 
AR Gear hee karcnrcucis laa els 6 — 1-9 + 0-1 


ocean, where the frictional-vorticity tendency is very small. However, on 
the western side, where the velocity and shear are great because of the 
narrowness of the western meridional current, the frictional- and planetary- 
vorticity tendencies are greatly enhanced, to orders of magnitude greater 
than that of the tendency of wind-stress vorticity, say +10 and —10 
respectively. In this way we may achieve a balance between the tendencies 
in the western side of the ocean as well, and so obtain a steady state. Table 3 
shows the balance of terms in an asymmetric circulation. 


TABLE 3 


Vorticiry TENDENCIES IN AN ASYMMETRIC CIRCULATION 


North-flowing South-flowing 
Vorticity tendency currents in the currents over remainder 
western edge of ocean 
WimGestrose (is 5 oie dec lee ed ole — 1:0 — 1:0 
SI INTL eA tae ay opty oss ebeks tis) o veveicyelors + 10-0 + 0-1 
LISTED "S505 ed ogObde DoneaneE — 9-0 + 0-9 
pea uell eetetayefa crave sl exe vale evel wiaveiiehe 0-0 0-0 


If one analyzes ocean currents in this qualitative manner for anticyclonic 
and cyclonic wind systems in both hemispheres, he will find that the neces- 
sary intensification of the ocean circulation is always on the western side 
of the ocean. This simplified analysis applies only to the viscous theories 
of the Gulf Stream discussed in this chapter, not to the nonlinear theory 
presented in Chapter VIII. 


WESTWARD INTENSIFICATION 


In 1946, when I first began to study the North Atlantic circulation, Dr R. B. 

Montgomery suggested to me that the marked east-west asymmetry of the 

surface circulation was one of the main features requiring explanation. 
Soon afterward, I proposed a simple model of a wind-driven ocean 
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(Stommel, 1948), to show the important effect on the transport lines which 
results from the variation of the Coriolis parameter with latitude. 

A rectangular ocean is envisaged, with the origin of a Cartesian-coérdinate 
system at the southwest corner (see fig. 54). The y-axis points northward, 
the x-axis eastward. The shores of the ocean are at x=0,7, and y=0, b. The 
ocean is considered as a homogeneous layer of constant depth D when at 
rest. When currents occur, as in the real oceans, the depth differs from 


D everywhere by a small variable 
amount h. The quantity A is much 
smaller than D. The total depth of 
the water column is therefore D+A, 
D being everywhere constant, and h 
a variable yet to be determined. 

The winds over the ocean are the 
trades over the equatorial half of 
the rectangular basin, and prevailing 
westerlies over the poleward half. 
An expression for the wind stress 
acting upon a column of unit hori- 
zontal area and depth D+/ must 
include this dependence upon y. 
A simple functional form of the wind 


y 


x 

0 

Fig. 54. Codrdinates and boundaries 
of the rectangular ocean basin used by 
Stommel (1948, fig. 1), but with notation 
changes. The x-axis points toward the 
east; the y-axis to the north. The dimen- 
sions of the rectangular basin are r and b. 


stress is taken as — F' cos (77y/b). 

To keep the equations of motion as simple as possible, the component 
frictional forces are taken as — Ru and — Rv, where R is the coefficient of 
friction, and w and v are the x and y components of the velocity vector, 
respectively. The Coriolis parameter f is also introduced. In general, it is 
a function of y. 

The vertically integrated steady-state equations of motion, with the 
inertial terms omitted, are written in the form 


oh 
0=f(D+h)v—F cos" — Ru-g(D+h) = (1) 


Ox’ 


oh 

Dm ee lee ne OU ihe (2) 
The quantities w and v are taken to be independent of depth, an assump- 
tion which simplifies the analysis, but requires regarding the wind as 
essentially a body force instead of a surface stress. To these, the equation of 


continuity must be added: 


al(D+h)u]  A(D+h)e] 


Oa dy lt 


(3) 
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Cross-differentiation of the first two equations and use of the third result 
in the following equation: 
of Fn sin dv du 
pia ie R eae ; 
u(D+h), so ee art (=- =)= =0 (4) 
In the actual oceans, h is so much smaller than D that to a first degree of 
approximation this may be rewritten as 


S78 
photysns +a —5 =0; (5) 


where the following definitions have been made: / =0f/éy; and y= F7/Rb. 
This equation is called the vorticity equation. To the same degree of 
approximation, the equation of continuity may be replaced by 


ae Ov 
a (hy 6 

Ox oy Sak (6) 
A stream function 7 is introduced now by the relations: w=0y/dy; and 
v= —dy/éx. The vorticity equation is now rewritten in terms of the stream 
function: D oy 

5 UY 
CHIP Neca 2a i winiby 7 


The boundary conditions are that the shore of the ocean be a streamline 


¥(0,y)=W(r,y) = W(x, 0) = (a, 6) =0. (8) 


If f is a linear function of y, then (D/R)f is a constant. The general 
solution is 


(Nets) 
ae ery | ea 9 
i XY, (2) sin= > (9) 
where 
Y=) (c; sinn,y +d; cos n,y) , (10) 
X=) (p,e4i* + 9,ei") . (11) 


The constants A; and B; have been defined thus: 


Ay=—5h+ (se) +nt}, and peu (55) nj}. (12) 


The quantities c,, d;, p;, q; are undetermined constants. This solution is very 
general, but reduces to a simple closed form when the boundary conditions 
are imposed. First of all, the d; and c; vanish, except c, corresponding to 


90 LINEAR THEORIES—VISCOUS 


n,=7/b. This constant c, may be absorbed into p, and qg,. When subscripts 
are dropped, the stream function has the form 


2 
= (5) sin“? (pet#+-ge"*—1), (13) 
where 
1—eBr 
Par ore and g=1-—p. (14) 


The curves (77 =const.) are the streamlines of the ocean currents. 

In order to help visualize the meaning of this solution, it is advisable to 
compute some numerical examples that will show what role the various 
parameters play. Three cases are discussed. All involve the same effects of 
wind stress, bottom friction, and horizontal pressure gradients caused by 
variations of surface height. The role of the Coriolis force is different in each 
case. First it is assumed that the Coriolis parameter vanishes everywhere— 
the case of the nonrotating ocean. Secondly, it is assumed that the Coriolis 
parameter is constant everywhere—the case of the uniformly rotating 
ocean. In the third case it is assumed that the Coriolis parameter is a linear 
function of latitude. Of the three cases, the last one most nearly approxi- 
mates the state of affairs in the real ocean. 

For convenience of the numerical computations the dimensions of the 
ocean are taken as follows: 


r =10° cem.=10,000 km., 
b =27 x 108 cm. =6,249 km., 
D=2.=x 10* cm.=200 m. 


The maximum wind stress F' is taken to be 1 dyne/cm.?. 

The coefficient of friction R is the only quantity for which a value must 
be devised. If a value of R =0-02 is assumed, the velocities in the resulting 
systems approach those observed in nature. 

The case of the nonrotating ocean.—In the nonrotating ocean the con- 
stants p and q are fairly simple. Within 1 per cent, or as closely as graphs 
may be drawn, p and q are given by 


pe g— (15) 
The equation for the stream function is therefore 


2 


The east-west and north-south symmetry of the streamlines is im- 
mediately evident from this equation. The actual streamlines computed 
from it are exhibited in fig. 55. 
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The height contours are computed by integration of the primitive 
equations and are plotted in fig. 56. The general features of the nonrotating, 
wind-driven system constitute a broad circulation exhibiting absolutely no 
tendency toward crowding of the streamlines. 


— 


1000 km 


Fig. 55. Streamlines of a nonrotating ocean (Stommel, 1948, fig. 2). 


Fig. 56. Contours of surface height in the nonrotating ocean (Stommel, 1948, 
fig. 3). 


The case of the uniformly rotating ocean.—If the Coriolis parameter is a 
constant 0-25 x 10-4, the streamline diagram does not differ from that of the 
nonrotating basin. When the height contours are computed, however, a 
difference between the two cases becomes apparent, as shown in fig. 57. 
The large elevation in the central part of the ocean provides horizontal 
pressure gradients that nearly counterbalance the Coriolis forces. The 
height contours are not strictly parallel to the streamlines, but nearly so. 

Coriolis parameter a function of latitude.—In the real ocean the Coriolis 
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force is a function of latitude. In low latitudes this function is nearly a 
linear one: f= fy, ? = 10-/em. sec. The inequality in the absolute values 
of the quantities A and B that occurs in this case immediately makes clear 
the complete lack of east-west symmetry. The streamlines drawn from this 
formula are shown in fig. 58. The most striking feature of this figure is the 


—oo 


1000 km 
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Fig. 57. Contours of surface height in a uniformly rotating ocean (Stommel, 
1948, fig. 4). 


+ — 


1000 km 


Fig. 58. Streamlines for the case in which the Coriolis parameter is a linear 
function of latitude (Stommel, 1948, fig. 5). 


intense crowding of streamlines toward the western border of the ocean. 
The rest of the streamline picture is broad and diffuse. The resemblance that 
the velocity field of this simple case bears to that of the actual Gulf Stream 
suggests that the westward concentration of streamlines in the wind-driven 
oceanic circulation is a result of the variation of the Coriolis parameter 
with latitude. 
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The height contours computed from this example are shown in fig. 59. To 
extend the results of this study to the Southern Hemisphere, the reader will 
notice that since # is unaffected by crossing the equator and F simply 
changes sign, all the diagrams may be transformed to below the equator by 
simple reflection across the x-axis. The crowding of streamlines is therefore 
toward the western border of each ocean, irrespective of hemisphere. 

The artificial nature of this theoretical model should be emphasized, 
particularly the form of the dissipative term. 


— 


1000 km 


Fig. 59. Contours of surface height which are associated with flow shown in 
fig. 58 (Stommel, 1948, fig. 6). 


MUNK’S THEORY OF THE WIND-DRIVEN 
OCEAN CIRCULATION 


By far the most distinguished theoretical investigation into the wind-driven 
ocean circulation is that conducted by Munk (1950), who has succeeded in 
deducing many features of the mean ocean circulation from the wind stress 
alone. In the preceding section of this chapter we have discussed qualita- 
tively the reason for the great intensification of western currents; in this 
section we shall follow Munk’s quantitative development of the theory. 

Let us suppose the ocean surface to be a plane surface at rest: x is eastward, 
y is northward, and zis upward. The plane z=0 lies at the mean surface; the 
actual surface is at z=Zp. 

The equations of steady motion in the horizontal plane are: 


Op a C* 0 ou 
= Ke 17 
pie B+ Bal satis) Ye ( re) 9) 


op Pam eaNS |. / Lee 
ayo ea 18 
pfu et Kast) +5 rs} oe 
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where uw and v are the x and y velocity components, : is pressure, f is the 
Coriolis parameter, p is density, and Ky, and Ky, are the horizontal and 
vertical coefficients of eddy viscosity. We assume K, to be constant and 
uniform. 

These equations may be integrated from the surface z=z, to a depth 


z=-—h beneath which both the currents and the horizontal pressure 
gradients vanish: 2) 
P= | pdz, (19) 
M,= | " pudz, (20) 
—h 
u,= | " pede. (21) 
—h 
Now the integrals of the pressure gradients are given as follows: 
% Op oP 02 
| * dea F pled, (22) 
% Op oP C2 
— dz=—— p(%)—. 23 
[° pam GPS (23) 


There seems to be little reason to make use of the complete integral of the 
horizontal shearing-stress terms. We shall assume that 


Zo C2 C2 C2 C2 
J" alge a) = Ug) Me = 
and that 
Zo g2 2 gz gf 
J Bal setae) =A (Get ge) Mo 
where 
A = Ky . 
The vertical shearing-stress term integrates very simply: 
#0 ou 
[" g(&rg) ear (26) 
% 0 ou 
Be) iy eee | 
ee ( 5) a2 70, (27) 


where 7, and 7, are the x and y components of the wind stress applied at the 
surface, The integrated equations of motion are: 
a. a 


oP 
ea) |e 7 t4(gat dy? 


si) Mette, (28) 


oP (@ a 
+M,f= 5 t4(sat sya) My +Ty. (29) 
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We may now introduce a transport function wy, defined by M,= — dy/ey, 
M,=0y/0x, and eliminate P by cross-differentiation. We shall assume that 
f=Of/oy is a constant, and thus obtain a simple linear equation for yy: 


Gi eek | BN 8 Or, or 
A Ze p= 2. | 
(Fata ait a) e =| > Bose ie) 


The boundary conditions are that both 7 and its derivative normal to the 
boundary shall vanish. If the boundaries (coastlines) are taken as forming 
a simple rectangle x=0,7, and y= +8, and only an east-west wind-stress 
system is assumed, 7,=0, then an approximate solution is of the form: 


OT 
ly = ee II 
p=rXp By’ (31) 


where 


/3 (3 7 1 ee 
— — Bo-(1/2) ka ~ —— (ka —e-hr-2) — ] 
Xx Be cos (5 kato 6 +1 Jy (Oe € ); (82) 


where B= (2/,/3) —(./3/kr), and k=3,/(f/A). 

The w-field which Munk (1950, fig. 2) computed for a rectangular basin 
of Pacific Ocean dimensions, using mean wind data in the fashion described 
by Reid (1948a), is shown in fig. 60. If we consider mean annual zonal 
winds only, we see that the integrated oceanic wind-driven circulation is 
divided into closed circulatory systems—or ‘gyres’, as Munk calls them— 
bounded at latitudes ¢,, where curl, 7=0, and with latitudinal axes at 
latitudes ¢,, of extreme values of curl, 7, which of course do not necessarily 
correspond to latitudes where 7 =0. The Sargasso Sea in the Atlantic Ocean 
is thus situated at the inflection point of the mean wind stress between 
westerlies and Northern Hemisphere trades. 

The function X(x,) along the latitude circles ¢, can be interpreted as the 
total northward transport of the current between x=0 and «=, whereas 
the quantity X’(x) is the transport per unit width. Fig. 61 shows these 
functions drawn to an arbitrary scale. The region from «=0 to x=4/k is 
supposed to correspond to the Gulf Stream. 

When X and X’ are computed it is found that the equations fall naturally 
into three parts, each of which dominates in a given sector. At the western 
edge of the ocean x<r, and 


a oe kt eos (5 ka — 4 +1, (33) 


w —(1/2) ka ai v 
k => 3B é sin 5) ka 5 ( ) 
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representing slightly wnder-damped oscillations of wavelength 


47 47 


Wie Be apo (35) 


Table 4 gives the locations and values of the first few extrema. East of the 
main current there is a countercurrent of the magnitude exp (—7/,/3), or 
17 per cent of that of the main current. From the analyses of actual data 
by Wiist (1936) and Iselin (1936) it appears that this theoretical counter- 
current is very nearly equal to that observed. In Chapter VIII it is sug- 


CENTRAL DRIFT 
Exe 
Ck) k 


15 10 heron) 5 ie} 


Fig. 61. Plot of equations for X and X’/k giving west-east variation in 
transport (~X) and transport velocity (~X’) from the western shore (cx=0) 
to the eastern shore (l1—a2=0). The scale of the central and eastern solutions is 
exaggerated relative to the scale for the western solution. From Munk (1950, 
fig. 3). 


gested that the Gulf Stream is not essentially a frictional phenomenon. 
Frictionless models are described, but none give the remarkable similarity 
to the countercurrent which the Munk theory yields. 


TABLE 4 


EXTREMA OF Xw AND Xiy 


Location# 
Measurement Western Western Counter- Counter- 
current current current current 
axis limit axis limit 
hl oi ee 1/6 3/6 4/6 1 
Dy ests tic SRO OR RRO 0-45 1:17 1-09 0:97 
w/k. Biete te wicre sche 0°55 0-00 — 0-09 0-00 


® See fig. 61. 
Source: Munk (1950, table 1). 
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The shape of the X curve for the western current shown in fig. 61 may be 
compared with the transport function y computed geostrophically from 
hydrographic data for both the Gulf Stream and Kuroshio (fig. 62). The 
east-west scale of the theoretical curves depends upon the particular choice 
of k involved, hence of lateral-eddy viscosity. In order to produce a Gulf 


10°"? ¥ in g./sec. 


100 200 300 km. 


Fig. 62. The mass-transport stream function y computed from oceanographic 
data across the Kuroshio and Gulf Stream. Compare these with the left-hand 
side of fig. 61. From Munk (1950, fig. 4). 


Stream 200 km. wide, Munk chooses A =5 x 10? em.?/sec. A smaller value 
of A would give a narrower Gulf Stream. 

The more recent cross sections of the Gulf Stream (e.g., fig. 33) suggest 
that the width of the main current used by Munk is three or four times as 
great as it should be. Thus, for a purely frictional theory the coefficient of 
lateral-eddy viscosity must be reduced to about 10% em.?/sec. Although this 
adjustment of the parameter A is perfectly admissible, there being no in- 
dependent measure at present, the inertial terms are no longer negligible, 


— 
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and, as we shall see in the next chapter, there is therefore much reason to 
formulate the boundary layer by using inertial terms rather than friction. 
For the moment, however, we follow Munk’s development. 

The total transport of the western current, ~,,, is independent of A, 
and, using the value of X,, given in table 4, is simply 


Veg —l-Lirp + curl.7. (36) 


Table 5 contains a comparison of the computed transports of various 
western currents: those computed by Munk (1950, table 2) from the wind 
stress, and those computed from oceanographic observations. 


TABLE 5 


THE Mass TRANSPORT OF SOME WESTERN CURRENTS, DETERMINED FROM 
THE WIND STRESS AND FROM OCEANOGRAPHIC OBSERVATIONS 


103? y 
1018 f/ aa aa 
Garvont Lati- | om.) an aiieeedl Wind | Computed | Oceano- 
esl lc sae) facts | ates | epsceentions 
g./sec. g./sec. 

Gulf Stream ..| 35° N. 1-9 6,500 70 1-30 36 74 

Kuroshio Ext..| 35° N. 1-9 |10,000 50 1-25 39 65 

Oyashio Cur...| 50° N. 1:5 5,500 — 15 — — 6:5 —7 

Brazil’Cur... ..| 20° S. 2-2 5,500 — 20 — — 5:8 |— 5to — 10 


Source: Munk (1950, table 2). 


The effect of the meridional wind component, which supports the zonal 
winds, is introduced by a meridional wind factor 


el 
—curl,7 (7) : (37) 
dy 


The computed transports of the Gulf Stream and Kuroshio are too small 
by about a factor of two. For example, Munk’s computed value of the trans- 
port of the Gulf Stream, using £=1-9 x 10-3/em./sec., r=6500 km., and 
dr,/dy=70 x 10-" g./sec./sec., is 36 x 1012 g./sec., as compared to an actual 
value of about 74 x 10” g./sec., if the reference level be taken as the depth 
of 2000 m. 

The source of this discrepancy is not to be sought in the physics of the 
western current itself. The results shown in table 5 apply equally well to 
the nonlinear theories set forth in Chapter VIII. Munk (1950, p. 92) 
ascribes the discrepancy to an underestimate of wind stress for low wind 
speeds: 


Finally we may examine the question of why the computed 
transports of the Gulf Stream and the Kuroshio current amount to 


9-2 
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only about one half the observed values [see table 5]. It does not 
seem reasonable that the oceanographic observations should be 
off by more than, say, 20 per cent, nor that the theoretical ex- 
pression. . .for the transport should account for the discrepancy, 
since it is almost independent of the eddy viscosity and the shape 
of the ocean basin. 

Maury and others have ascribed the North Atlantic circulation, 
in particular the Gulf Stream, to differential heating between 
equator and pole, to the freezing of ice, and to other processes that 
make up the thermohaline circulation. If we assume that the circu- 
lation were half wind-driven, half thermohaline, it would be a 
strange coincidence that the general pattern of the circulation, such 
as the boundaries of the gyres, should conform so closely to the 
general atmospheric circulation. Furthermore, Fuglister [Munk 
refers to an unpublished manuscript, later published, 1951a] has 
found a high correlation between variations in the current with 
variations in the wind. It should also be noted that the thermo- 
haline circulation insofar as it is related to the outflow of river water 
along the Atlantic seaboard would tend to reduce rather than to 
strengthen the Gulf Stream. It would seem therefore that the sub- 
tropical gyre, and probably also the subpolar gyre, are predomin- 
antly wind-driven. 

Methods for computing wind stress from the observed wind 
speeds according to the equation 


T= Cp Pair U2 [38] 


are discussed by Reid [1948a]...Underestimates of 7 may first 
of all result from underestimates in the wind speeds on the climato- 
logical charts from which the appropriate averages were taken. 
The preponderance of coastal stations, and the tendency of ships 
at sea to avoid regions of high wind, would lead to consistent 
errors, but these cannot account for more than a fraction of the dis- 
crepancies. The weakest link is the drag coefficient Cp which is 
based on measurements of Baltic storm tides, and a few other 
measurements. . . In accordance with the views presently accepted 
we have assumed C,)=-0026 at high wind speeds, Cp, ~-008 at 
low speeds, with the discontinuity occurring at Beaufort 4 [Munk, 
1947]. In the trade-wind belt of the eastern Pacific, where the 
winds are predominantly Beaufort 4 and above, Sverdrup [1947] 
and Reid [1948a] have obtained satisfactory agreement between 
computed and observed transports. Supposing a value of -0026 
were applicable at all wind speeds, then the effect of the south- 
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westerly winds over the eastern Atlantic would make itself felt in 
a much higher meridional wind factor [table 5], and it can be de- 
monstrated that the discrepancy between computed and observed 
transports could be largely accounted for. We are therefore led to 
propose a higher value of C', at low wind speeds. The transports of 
the Gulf Stream and Kuroshio current are, after all, probably as 
good an indicator of the overall stress exerted by the winds on the 
ocean as any of the measurements on which the value of Cy is now 
based. 


There have been several attempts to obtain better agreement between 
computed and observed transports. I have expressed my own views in 
Chapter XI. Hidaka (1949) has performed the analysis of the linear theory 
of wind-driven ocean currents in many different forms. In one study, using 
spherical codrdinates, he has obtained a value of transport of the Kuroshio 
more nearly equal to that observed than has Munk. It is difficult to see 
how a change of codrdinates can make so great a difference in the transport. 
Recently, Sarkisyan (1954) has carried out a numerical study dynamically 
similar to Munk’s, but in an ocean shaped very much like the real North 
Atlantic. He obtained transports for the Gulf Stream of between 70 and 
90 x 10!" g./sec., but as I do not know the details of the wind-stress 
distribution which he used, the significance of his close agreement is not 
clear to me. 

Owing to the existence of a biharmonic operator in the viscous term of 
the governing equation, there is an exponentially decaying line of vortices 
to the east of the countercurrent, centered along the axis ¢, (fig. 63). There 
is no good evidence that such a line of vortices actually exists. By some 
stretch of the imagination one might find confirmation in the charts of 
Felber (1934) and Defant (1941), but to my mind this is rather special 
pleading. 

Beyond the vortices, over most of the central regions of the ocean the 
solution reduces to 


(39) A,=1——; —=——, (40) 


This corresponds, for zonal winds, to the solution of the equation given by 
Sverdrup (1947, p. 322): 
4) > 
ae i= carl, (41) 
which does not contain lateral-stress terms. This important relation is dis- 
cussed further in Chapter XI. 
Munk and Carrier (1950) have also investigated the effects of variously 
shaped ocean basins—in particular, a triangular one, which fits the North 


102 LINEAR THEORIES—VISCOUS 


Pacific Ocean more properly than the rectangular model. Munk, Groves, 
and Carrier (1950) have studied the effect of the nonlinear inertial terms, 
assuming they are small, using a method of successive approximations, and 
using Reid’s (19486) model of vertical density structure, which consists of 
an exponential decrease of density upward to the thermocline, and a 
homogeneous upper layer. They found a slight downstream displacement 
of the region of maximum currents. The two important observed features 
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Fig. 63. The volume-transport stream function, in units of 10° m.%/sec., near 
the western boundary for mean annual zonal winds over the Atlantic. The 
center line is at 31° N. For comparison with the Sargasso Sea circulation, the 
figure should be distorted by maintaining the west—east orientation of the 
x-axis and rotating the y-axis clockwise until it coincides with the coast. From 
Munk (1950, fig. 5). 


which were not accounted for in the linear theory, (i) the inshore counter- 
current and (ii) the continued sharpness of the Gulf Stream long after it 
leaves the coast, were not found in the higher-order solutions. Sub- 
sequently, Miyazaki (1952) found that an inshore countercurrent (i) can 
be obtained formally by letting the coefficient of eddy viscosity decrease 
toward the coast, and, as we shall see in the next chapter, there is no 
difficulty in obtaining feature ii if the inertial terms are large enough. 

So much of the quantitative aspects of the theory of wind-driven ocean 
currents depends upon the actual pattern of wind stress used that it is 
advisable to reproduce for reference the actual distribution of zonal wind 
stress which Munk has adopted (fig. 64). 
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Fig. 64. Distribution of zonal wind stress with latitude. (Zonal wind stress 
T, in dynes/em.?.) Note the longitudinal asymmetry over the North Atlantic. 
In the region south of Greenland, where observations are scarce, the stress 
averages have been corrected for this asymmetry. Dotted lines indicate curves 


plotted from scanty data. From Munk (1950, fig. 9). 


Chapter Eight 


NONLINEAR THEORIES OF 
THE GULF STREAM 


This chapter records attempts to include the effect of the nonlinear inertial 
terms in theoretical explanations of the Gulf Stream. From a historical 
point of view it is important that mention be made of the wake-stream 
theory advanced by Rossby. Rossby’s theory of the Gulf Stream assigns 
a very important role to lateral friction, just as the linear theories dis- 
cussed in the preceding chapter do, but does not take into account the 
effect of the variation of the Coriolis parameter with latitude. Next, we 
shall consider some of the more recent nonlinear theories in which lateral 
friction is assumed to play only a minor role. We therefore reéxamine 
the question how large lateral mixing and viscosity really are, and we 
find that although a definite answer cannot be given, there are grounds 
for supposing that lateral friction may actually be small as compared to 
inertial terms. 

We shall derive a principle of conservation of a quantity called ‘potential 
vorticity’ in order to incorporate the inertial terms in a convenient manner, 
and then make a vorticity analysis of a Gulf Stream section. After a few 
preliminary remarks about certain very simple potential-vorticity models, 
we shall proceed to a study of the recent ‘inertial boundary-layer’ theories 
of Morgan and Charney. Finally, reference will be made to Rossby’s 
concept of ‘critical’ flow as applied to the left-hand edge of the Gulf 
Stream. 
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ROSSBY’S WAKE-STREAM THEORY 


Fourteen years before the development of the linear theories discussed in 
the preceding chapter, Rossby (1936a) worked out an interesting nonlinear 
model of the Gulf Stream in which large-scale lateral mixing was involved. 
The geostrophic relationship alone does not prescribe a transverse velocity 
profile for the Stream; and hence Rossby made use of the ideas of a turbulent 
jet stream from the field of experimental fluid mechanics to define the 
profile. The Straits of Florida were supposed by Rossby to act as a nozzle, 
or jet; the Gulf Stream was therefore regarded as a purely inertial stream, 
interacting with its environment by mixing. 

Rossby reduced the problem to one on a nonrotating system, by ex- 
tracting from the equations the Coriolis-force terms and their associated 
pressure gradients, and by assuming that the motion is essentially non- 
divergent. The problem was thus reduced to the one that had already been 
discussed by Tollmien (1926), who was able to obtain good agreement with 
experimental results by assuming (i) that the pressure gradient along the 
jet is small, and therefore that the total momentum transport of the jet is 
the same at all sections downstream; and (ii) that the shearing stresses 
acting upon the jet may be described in terms of a mixing length propor- 
tional to the distance from the nozzle. 

Tollmien found that the mass transport and width of the jet increase 
downstream. This requires that there be an inflow of water from the sur- 
rounding medium into the jet. 

Although the mass transport of the Gulf Stream does increase after 
leaving the Florida Straits, it does not increase after passing Cape Hatteras. 
The angular spread of experimentally produced jets (Férthmann, 1934) 
varies greatly, but Peters and Bicknell (1936) obtained spreads of between 
8 and 14°. The spread of the actual Gulf Stream in its instantaneous form is 
less than 1°, although the mean spread averaged over many different sets of 
observations is greater (Stommel, 1951). 

After discussing the analogy to Tollmien’s jet in a homogeneous ocean, 
Rossby (1936a) made some qualitative studies of the effect of the stratifi- 
cation of the sea, by considering an ocean made up of two layers of slightly 
different density, the lower layer being at rest. By the geostrophic relation- 
ship, the difference in level of the interface at the two sides of the jet stream 
is a measure of the total mass transport of the stream; therefore, if the mass 
transport increases downstream, it soon becomes impossible to join the 
interface inside the stream to the level of the undisturbed water masses 
outside the stream, unless the interface is bent back (up on the right, down 
on the left) to these levels. This requires, geostrophically, a countercurrent 
on each side of the jet. 
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To carry the argument much further becomes difficult, because one has to 
allow for the possibility that changes in the stratification affect the mixing 
length in some way about which we can only guess. Since we do not know 
what these laws actually are, we may easily be led astray. Rossby (1936a) 
is led to argue that there should be a countercurrent on the inshore side of 
the Gulf Stream, but none on the Sargasso Sea side. In fact, the counter- 
current observed on the inshore side is very much weaker than the pro- 
nounced countercurrent on the Sargasso Sea side. Rossby’s treatment of 
the dynamics of the countercurrents, in the wake-stream model, has been 
very stimulating. We shall use the same dynamical equations later in this 
chapter, when we come to discuss the Stream as a stream of uniform 
potential vorticity. 


LATERAL MIXING 


The Rossby wake-stream theory and the linear theories discussed in the pre- 
ceding chapter all depend heavily upon the existence of large-scale turbu- 
lence: lateral mixing and the associated lateral-eddy viscosity. Whenever 
we attempt to frame a hydrodynamical problem in terms of eddy viscosity 
we must prescribe a certain scale of motion, to divide motions into what we 
willregard as mean motions and turbulent motions. Our dynamical equations 
are expressed explicitly in terms of the mean motions only; the turbulent 
motions are lumped statistically into one parameter: the lateral-eddy 
viscosity. It is clear that in a medium in which all scales of motion are 
present the distinction between mean motion and turbulence is purely 
arbitrary, but that it is of utmost importance to have a clear mental picture 
of the scale of motion which divides the two forms of motion in every 
theoretical discussion. Thus, as we have seen, the actual Gulf Stream 
meanders; a space and time average over several months would doubt- 
less look rather like the gradual fanning out past Cape Hatteras that both 
Munk and Rossby seem to have in mind. Individual eddies and meanders 
would be statistically averaged out as turbulence; we may speak of the 
mean motion in this situation as the climatological-mean Gulf Stream. The 
pertinent scale of motion which separates mean motion from turbulent 
motion in this case is about 100 km., the half wavelength of meanders. As 
I have shown (Stommel, 1951, 1953), it seems likely that the eddies which 
break off from the Stream (for example, that represented in fig. 30) carry 
enough momentum to produce lateral shearing stresses of the order which 
Munk has postulated. But this is true only if we regard Munk’s theory as 
applying to the climatological-mean Gulf Stream. 

Actually, of course, we are not satisfied with such a coarse statistical 
description of the Stream as the one to which the 100 km. scale of averaging 
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limits us. We want to describe more minute features of the Stream: a 5 or 
10 km. scale is more appropriate to our purpose. In this case, of course, the 
meanders and eddies cannot be regarded as turbulence, and the mean 
motion exhibits the structure observed by a single hydrographic section: 
we speak of this narrow filament of moving water as the instantaneous Gulf 
Stream. It is the Stream described in Chapter V. 

In attempting to frame a theory of the instantaneous Gulf Stream we do 
not have meanders and eddies at our disposal to provide large lateral 
stresses. As was mentioned in an earlier chapter (pp. 62 ff.), water-mass 
analysis indicates that there is not much mixing on a scale smaller than 
10 km. In an attempt to obtain some quantitative information on lateral 
shearing stresses, I recently (Stommel, 19556) analyzed a set of current 
observations made by Pillsbury in the Florida Straits, and found that the 
coefficient of eddy viscosity was probably less than 10° cm.?/sec., a value less 
than 2 per cent of the coefficient employed by Munk. Of course, these 
observations were not made in the Stream beyond the Straits, and the 
coefficient for that region is the one which we are really interested in 
knowing, but there are no long series of observations made there. Thus it 
seems difficult to discover enough turbulence of a scale smaller than 10 km. 
to provide important lateral-eddy stresses in the instantaneous Gulf Stream. 
We should emphasize here that apparently Munk was under the impression 
that he was dealing with the instantaneous Gulf Stream; his comparisons 
with hydrographic data indicate that. As I have suggested here, however, 
there is reason to suppose that his theory really applies to the climatological- 
mean Gulf Stream. We shall therefore attempt to frame a theory of the 
instantaneous Gulf Stream which does not take account of any important 
friction, but treats the Stream as essentially an inertial phenomenon. 

Before leaving this vexed question of the importance of lateral stresses 
we should note an argument proposed by Rossby (1936a) and Montgomery 
(1940) and recently questioned by Morgan (1956). Rossby asserted that 
since there is little bottom friction in the ocean, the torque of the wind 
stress applied to the surface of the sea can only be balanced by the torque 
of lateral shearing stresses around the coasts of the sea. Morgan (1956) has 
shown that there are other torques that need to be considered: a Coriolis 
torque, and a torque produced by differences of pressure at different parts 
of the coast. The relative importance of these various torques has not yet 
been ascertained; but much of the moral support for the indiscriminate use 
of lateral mixing and eddy coefficients has been removed. 

Let us now explore some ideas which may provide an explanation of the 
detailed structure of the Gulf Stream independent of friction. 
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THE POTENTIAL VORTICITY OF A LAYER 
OF UNIFORM DENSITY 


In order to include nonlinear dynamical terms and density stratification, it 
is helpful to introduce a quantity called potential vorticity, obtained in the 
following manner. Using the coérdinate systems employed in the preceding 
chapter, let us consider a homogeneous layer of fluid with density p, and 
thickness D, in which case the dynamical equations are 


dt?” pon ‘ (1) 
dv 1 Op 


where d/dt = (8/0t) + u(d/éx) + v(2/2y), and X and Y include both frictional 
driving and retarding forces, which we are not now interested in considering 
in detail. The Coriolis parameter f is regarded as a function of y. First, we 
eliminate the pressure gradients by cross-differentiation: 


Ou +5)= OY ox 


qirOrU+O(S+e =a a c) 


where ¢, the relative vorticity, is defined as = (dv/0x) — (eu/éy). The quantity 
f+(is called the absolute vorticity. The equation of continuity of the layer is 
dD Ou dv 
sss ay s\ bate abe) Gear 
dt ‘i (." +5) ‘ (4) 


Elimination of the horizontal divergence between the two equations results 
in the equation of potential vorticity, (f+ ¢)/D, 


dt w= D 5 (ae~ ay) (9) 


If there are no frictional forces, this equation simply states that the 
potential vorticity of a water column within the homogeneous layer cannot 
change as it moves from one place to another. Thus, as a column moves 
from one latitude to another, an adjustment must occur in the depth of the 
layer and in the relative vorticity in such a way as to keep the potential 
vorticity constant. 


THE POTENTIAL VORTICITY IN AN 
ISOPYCNAL LAYER 
Any actual Gulf Stream cross section exhibits a marked decrease in the 


thickness of layers (bounded by isotherms) in the high-velocity parts of the 
current. This vertical shrinking is particularly noticeable in the isothermal 
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layers warmer than 16°C. It seems dynamically significant that this 
vertical shrinking is just about the amount required for the potential 
vorticity of each isothermal layer to be uniform across the Stream, all the 
way from the inshore edge of the Stream, where vertical shrinking is at a 
maximum, to the Sargasso Sea. At the left-hand edge there is evidently a 
sharp discontinuity in potential vorticity. This is also true of the T-S 
correlations of the surface water. 

In order to demonstrate this remarkable uniformity (Stommel, 1955 a) of 
potential vorticity, let us consider the vertical thickness h of the layer 
bounded by the 17 and 19°C. isotherms in Worthington’s Gulf Stream 
section shown in fig. 33. We place the y-axis in the direction of flow of the 
stream, and direct the z-axis toward the right of the current, extending 
eastward into the Sargasso Sea. 

In the Sargasso Sea, on the right-hand side of the section, the depth of 
the isothermal layer is constant, and we denote it by hy. The velocity of the 
current vanishes there. Moreover, since the current is all in the y-direction, 
the relative vorticity is essentially given by dv/éx. The equation for uniform 


potential vorticity may then be expressed in the following way: / 
Ov 
wet (6) 
eT he 


The velocity v of the water in this isothermal layer may then be computed 
from the observed values of thickness h by direct integration: 


w@)={s(j--1) Of. (7) 


This computed velocity may now be compared with the geostrophically 
computed velocity. The results of such a comparison are shown in fig. 65. 
The fact that the velocities computed in these two very different and 
dynamically independent ways agree is good evidence that the potential 
vorticity in the Stream is nearly independent of cross-stream position. 


MODEL WITH UNIFORM POTENTIAL VORTICITY 


It is interesting to see how close a representation of the Gulf Stream can be 
obtained by using nothing but the principle of conservation of potential 
vorticity and the geostrophic approximation in a simple two-layer model. 
Taking the axes as before, we assume that a resting layer of uniform depth 
Dg, and density p,, on top of another resting layer of very great depth, and 
density p,, extends indefinitely in the positive x-direction. We now suppose 
that the interface between the two layers comes to the surface along the 
y-axis; thus D=0 at x=0. Therefore there must be a geostrophic current 
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Fig. 65. Vorticity analysis of Worthington’s section (19546; see figs. 33 and 
34 of the present study). The vertical thickness of the 17—19° layer is indicated 
in the upper graph. In the lower graph, the velocity computed by the vorticity 
formula (7) is shown by the solid curve; the geostrophic velocity as computed 
by Worthington, by the other curve. The fact that these two curves are fairly 
similar suggests that the potential vorticity in the upper layers of the Gulf 
Stream is fairly uniform across the Stream. For the dynamical significance of 
this, see the text. 


(Gulf Stream) in the upper layer just to the right of x =0, but, by hypothesis, 
it must not extend very far into the central Atlantic (that is, v=0 for 
large x). Since the lower layer is regarded as being at rest everywhere, the 
current velocity is given by the geostrophic relationship in terms of D: 


oD 
fo=g's (8) 
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where g’ = 9(P.—/)/P2. The configuration we have conjured up thus far does 
specify the total transport 7' of the stream, 


r=["vpae= SO (9) 
0 f2 

but does not specify the shape of the stream profile, the width of the stream, 
or whether countercurrents, and so forth, exist. The stream might be wide 
or narrow. Therefore some further physical statement must be introduced, 
to determine the stream completely. Suppose that we now require that the 
potential vorticity of the upper layer be the same for all z: 


ov 
f+ 
Ci oth 
7 =e: (10) 


Combining this requirement with the geostrophic equation, we are led to 
the following equation in D: 


53 (D—Dy): (11) 


where A?=g’ D,/f?. Therefore the cross-stream profile of depth and velocity 
is completely determined: 
D=D,(1—e-*), 


12 
v=(g'Do) e*". . 


If we now introduce numerical values D, = 800 m., and (p,—p,)/p,=2 x 10-3, 
the model gives a realistic transport, 7’'=64 x 10®°m.3/sec. The maximum 
velocity of the Stream is 4 m./sec. The width of the Stream, taken as the 
distance between the point of maximum velocity and that at which the 
velocity is reduced to 1/e of the maximum, is 40 km. The length, A, has been 
called by Rossby (1936a) the radius of deformation. Rossby introduced it 
in the presentation of his wake-stream theory, to explain countercurrents. 
Today it would seem that this term may actually be suitable for describing 
the Stream itself. 

Since this very crude model of uniform potential vorticity is independent 
of y, that is, of latitude, there is no provision for the growth of the stream, 
such as is observed in the real Gulf Stream between the Florida Straits and 
Cape Hatteras. In the actual Gulf Stream the water present at any latitude 
y’ is mostly drawn from lower latitudes, and hence one might reasonably 
ask why there should be any reason to expect the potential vorticity in the 
Stream at y’ to correspond to that outside the Stream at y’. The justification 
for applying the simple model described rests on the observed fact that in 
the central Atlantic between 10 and 35° N. latitude, the potential vorticity 
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of the top isothermal layers (> 10°C.) isremarkably constant. Suppose, for 
example, that the quantity D, is the thickness of the layer between the sea 
surface and the 10°C. isothermal surface. The depth of the 10° C. isotherm 
is shown in fig. 66. Since the relative vorticity in central oceanic regions is 
small, the reader may compute f/D, for the region between 10 and 35°N. 
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Fig. 66. Depth of the 10°C. isothermal surface in the western North 
Atlantic Ocean, according to Iselin (1936, fig. 47). Depths are given in meters, 


and convince himself that over this large part of the ocean the potential 
vorticity so defined is actually uniform and nearly constant. From a 
theoretical point of view it might be preferable to use density structure 
(see figs. 10-15). 

Dr George Morgan has pointed out to me that the easiest way to see the 
difference between the linear theories (the Munk theory, for example) and 
possible nonlinear theories is to interpret the potential-vorticity equation 
in the following way. The linear theory is steady, and essentially approxi- 
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mates the total derivative of the potential vorticity, (d/dt) [(€+f)/D], by 
fv/D. Thus € is regarded as small, and no account is taken of variations 
in D. The left-hand member of the potential-vorticity equation, formula (5) 
of this chapter, must balance against something, and hence all that remain 
are the various possible types of frictional terms. 

In the high-velocity regions of the Gulf Stream, the relative vorticity Cis 
negative and is not at all negligible as compared to f; inspection of the 
surface velocity profiles in fig. 32 will show the truth of this statement. 
Therefore it is quite conceivable that as a fluid column moves into the Gulf 
Stream, the changes in € and D can occur in such a way as to keep the 
potential vorticity constant and the conservation equation satisfied without 
the use of the frictional terms. Thus the theoretical model of a frictionless, 
nonlinear stream may describe the real Gulf Stream, the essential balance 
of terms within an isopycnal layer being of the following form: 


nse i eee at 


Wn sen. |p nab |=0. we 


A very interesting series of free, steady solutions of the motion in a 
frictionless, homogeneous ocean of constant depth has recently been 
obtained by Fofonoff (1954). The conclusion drawn from his work is that 
eastward-flowing currents cannot be broad, slow streams, but must be 
narrow and have high velocity and high relative vorticity. As yet, Fofonoff 
has not been able to obtain similar solutions for a two-layer ocean in which 
horizontal divergence is appreciable, but he has indicated that, in the North 
Atlantic, this will produce a westward, as well as a northward, intensifi- 
cation. When these solutions have been obtained we should expect the 
profiles of the strong currents to be similar to those obtained earlier in this 
chapter, in the consideration of conservation of potential vorticity. 
Fofonoff’s study is very similar to the uniform potential-vorticity model 
described above, and the models by Charney and Morgan described below. 
His paper (1954) was somewhat formal, however, and its applicability to the 
Gulf Stream becomes obvious only in the light of these and other studies. 

In formulating the physical model for a nonlinear steady-state theory of 
ocean currents we shall try to preserve the following characteristics. 

i. The ocean will be density-stratified. 

ii. The solution for central ocean transport will be the same as Sverdrup’s 
and Munk’s. 

iii. The narrow western stream (Gulf Stream) will be frictionless and 
convergent, and there will be a tendency toward development of a Fofonoff- 
type eastward flow in the region of the North Atlantic Drift. 

iv. All dissipation of the energy in the Stream will occur in a limited 
region toward the end, where meanders break up and multiple streams form. 

8 SGS 
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With respect to this region, the theory need not involve an assumption 
concerning lateral-eddy viscosity. The dissipation may be determined 
dynamically, as in a hydraulic jump in ordinary open-channel flow. 


THE REGION OF DECAY OF THE GULF STREAM 


After the Stream passes Cape Hatteras, it develops meanders, which increase 
in amplitude and eventually form detached eddies east of 60° W. longitude. 
The interpretation of observational data beyond this point is quite am- 
biguous, as has been mentioned in Chapter V. All indications are that the 
Stream breaks down into fragments and eddies which cannot be adequately 
described or surveyed by present means. Of all the regions of the Stream, 
this region of decay is the most difficult to depict and to understand. The 
region must be very extensive. The chart of the 10°C. isotherm depth 
(fig. 66) demonstrates that the potential vorticity of the upper layer is not 
even remotely uniform north of 35° N. latitude. The fragments, or, as they 
used to be called, ‘branches’, of the Gulf Stream still discernible in the 
North Atlantic Current (fig. 42) do not have the same potential vorticity, 
although there does appear to be a fair degree of uniformity within and to 
the south of any particular fragment. Apparently there are at least two 
complicated processes at work in this region which are not easily susceptible 
of simple theoretical treatment: 

i. The Stream becomes unstable and breaks up into fragments. Each 
fragment seems to carry an equal share of the transport. The potential 
vorticity of each fragment is uniform, but different from that of the parent 
Gulf Stream. These fragments can be traced at least halfway across the 
North Atlantic. They underlie the surface West Wind Drift. 

ii. Surface cooling and wind mixing profoundly modify the stratifi- 
cation of the surface layers, especially in the wintertime. Since, on the 
average, it takes more than a year for particles to pass entirely through 
the region of decay, it is reasonable to suppose that the potential 
vorticity of the waters in the region will be greatly changed by climatic 
factors. 

In the course of his fundamental studies of the atmospheric jet stream, 
Rossby and his colleagues at the University of Chicago have developed 
many interesting ideas which may have some bearing on the multiple- 
stream structures in the region of decay (1947). Neumann (1952) has also 
studied the region, and believes that the multiple streams are not decaying 
fragments of the Stream; he attributes them to persistent nonuniformities 
of the winds. My own impression is that wind patterns of the sort required 
do not persist long enough to produce such effects. 
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THE CRITICAL INTERNAL FROUDE NUMBER 


One of the important dimensionless numbers characterizing the flow of 
homogeneous frictionless fluids in channels with a free surface is the well- 
known Froude number, U?/gD, where U is the mean velocity of the fluid 
along the channel, g is the acceleration due to gravity, and D is the distance 
from the free surface to the bottom of the channel. According to ordinary 
hydraulic engineering practice, the flow is said to be subcritical, critical, or 
supercritical, depending upon whether the Froude number is less than, 
equal to, or greater than unity. In natural watercourses and engineering 
works such as canals and aqueducts, the flow is almost always subcritical 
except at certain control points (weirs and dams, for example), where the 
flow is locally critical. Whereas small perturbations in water level (long 
gravity waves) can move upstream in subcritical flow, they are brought to a 
standstill at points along the watercourse where the flow is critical. This 
explains the importance of points of critical flow in determining the water 
level upstream. Supercritical flow in nature is comparatively rare; it 
normally terminates abruptly in a hydraulic jump which converts the flow 
back to subcritical. 

Considering the great depth of the ocean, and the very moderate velocities 
of flow, it is clear that no ocean current even remotely approaches the 
ordinary critical Froude number. However, in a density-stratified fluid 
there are other ‘internal’ Froude numbers associated with the speed of 
propagation of long internal waves. In a system consisting of a very deep 
dense layer at rest, and a less dense surface layer of depth D moving at 
velocity U, the internal Froude number is U?/g’D, where g'=(Ap/p)g. In 
the stream of uniform potential vorticity discussed earlier in this chapter, 
in the theories of Morgan and Charney outlined later in this chapter, and in 
the simple theory of steady meanders presented in Chapter LX, parts of the 
flow are supercritical with respect to the internal Froude number. In the 
summer of 1955, in discussions at the Woods Hole Oceanographic Institu- 
tion, Dr Rossby pointed out that such regimes are likely to be unstable and 
probably do not exist in nature. Therefore we should not expect transverse 
profiles of velocity such as those given in equation (12) to hold for values of 
D/D, less than that for which v=,/(g'D). An easy calculation, for the 
stream of uniform vorticity, shows that this occurs at D=0-38 Dy. In view 
of the fact that this simple model of the Gulf Stream fits the data from 
observations so well, it is interesting to note that this is where the point of 
inflection in the depth of the 10°C. isotherm actually occurs (see fig. 33). 
Rossby (1951) has extended the concept of critical flow to flows with 
arbitrary vertical density structure. 

Up to the present, no one has succeeded in working out a detailed Froude 

8-2 


116 NONLINEAR THEORIES—INERTIAL 


theory which also takes account of the earth’s rotation. However, it does 
seem likely that the effect of the earth’s rotation will not radically alter the 
internal Froude criterion: the real Gulf Stream approaches critical internal 
Froude flow, and it is quite conceivable that as a result there are internal 
hydraulic jumps and other interesting small-scale phenomena, such as 
oblique shock fronts, along the left-hand, inshore, edge of the Stream. 


THE BOUNDARY-LAYER TECHNIQUE 


Although Munk’s first theoretical treatment of the wind-driven ocean 
currents involved complete analytical solutions for the entire ocean, he 
quickly realized that, since the viscous terms (the highest-order derivatives) 
were important only near the coasts, the problem could have been solved 
more simply by using Sverdrup’s (1947) simple first-order (in transport 
function) equations for the central parts of the ocean, and fitting a viscous 
boundary layer to the solution for the central regions at those places near 
the coasts where higher-order terms become significant. In the more com- 
plicated cases of triangular ocean basins this alternative procedure was 
adopted (Munk and Carrier, 1950). The ocean is thus divided into two 
regions: one, which we shall call the interior region, which includes all central 
regions of the ocean, and which is determined by the form of the wind stress; 
the other, the boundary-layer region, in which the higher-order terms 
are important. Actually, so far as the boundary layer is concerned, it does 
not matter what physical process produces the interior current regime. As 
long as some interior solution exists, be it wind-driven or thermodynamically 
produced, a linear viscous boundary layer can be fitted to it on the western 
coast. 

We shall take advantage of this natural separation of the problem into 
two parts in our formulation of the Gulf Stream as an inertial boundary 
layer. Furthermore, we shall limit the application of the theory to the 
region bounded by latitudes 10 and 35° N., the region of formation or growth 
of the Gulf Stream. North of 35° N. latitude the Stream breaks away from 
the coast, and we enter the region of decay. This region is not dealt with in 
either the Morgan theory or the Charney theory, both of which are dis- 
cussed in the next section. 


THE INERTIAL THEORIES OF MORGAN 
AND CHARNEY 


In the latter half of 1954 I made several futile attempts to construct a theory 
of an inertially governed Gulf Stream that would provide a more complete 
solution than the partial solutions described above. I confided my troubles 
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to DrJule G. Charney of the Institute for Advanced Study at Princeton, and 
to Dr George W. Morgan of Brown University, and in the spring of 1955 
these two investigators, independently, developed proper inertial theories 
of the Stream. Their theories lead to a nonlinear differential equation which 
can be integrated numerically. 

Let us first examine Morgan’s formulation. We consider an ocean bounded 
by meridional walls at x=0 and x=r, as in Chapter VII. For simplicity we 
also assume that the wind field over an interval of latitude 0 <y<s is given 
by the equation 


2 
r= —79(1- 4), T,=0. (14) 


According to the dynamical equation applying to the central part of the 
ocean [Chapter VII, equation (42); also Chapter XI, equation (5)] we can 
solve for the meridional transport per unit width, 1/,, in the interior: 


ous 2YTo 


=e= (15) 


On the east coast, =r, the zonal transport M,=—dy/0y vanishes; 
hence we have, from continuity, the following simple results for yy and W/,: 


= a2), (16) 
GR EO peg (17) 
a a Cee 


The advantage of the parabolic zonal wind-stress law, equation (14), is 
that it specifies a zonal transport at the western coast, «=0, independent of 
the latitude y. Solutions (15), (16), and (17) are independent of any assump- 
tion about the vertical density stratification. If we consider a homogeneous 
ocean, the variation in total depth due to slope of the free surface is negli- 
gible. If, on the other hand, we consider a two-layer system with density 
different in each layer, and the bottom layer at rest, we must also take into 
account the variation of the depth, D, of the top layer. In particular, the 
pertinent dynamical equation in the interior is the y-equation: 


On D 
~fm,=-95-(5). (18) 


Since, for the purpose of matching the boundary layer to the interior 
solution, we need only the value of D along the outer edge of the boundary 
layer (that is, for x=0), we can integrate equation (18) with respect to y, 
using M, at «=0 from equation (17) and taking into account the variation 
with latitude of the Coriolis parameter with latitude. We set f=fx+/xry. 
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In this way we obtain the following expression for use in the two-layer 
system: 


2M, 1 
D? (0, y) aa (fay +5 Fat?) +.D? (0, 0). (19) 


The quantity D (0, 0), the interior depth of the interface at the outer edge 
of the boundary layer, and at y=0, is a constant which we can choose 
arbitrarily so far as this theory is concerned. In this way we can form a 
variety of models. 

We now turn to the boundary layer itself. Morgan (1956) uses the 
following dynamical equations: 
gong @D* 
aa aeer (20) 
pr? 
“+ 9'D* = By"). (21) 


The asterisks are used to denote quantities defined within the boundary 
layer. Equation (20) states that the dynamical equation for force com- 
ponents normal to the coast is essentially geostrophic. Equation (21) is 
Bernoulli’s (nonlinear) equation for the upper layer. Next to the coast, u** 
is insignificant compared to v**. In the boundary layer the kinetic plus 
potential energies are a function B(y*) of the transport function only. 
Outside the boundary layer this is not so, because the wind stress has a long 
time to do work on the water there. It is convenient to transform these two 
equations into a form containing only y* and D*, by making use of the 
definition v* D* = 0y*/dx: 


yi g D* 
y= 5D +0), (22) 
oys* 
“= 2 D*[Biy*)—g' D*P”. (23) 
The quantity C(y) is easily evaluated along the outer edge of the boundary 


layer by matching y* and D* with interior solutions y and D at x=0 
[equations (16) and (19)]. In this way Morgan obtains the result: 


aa by’ M (0, 0) g' D*(0, 0) 
- | ii ai 


where we write M,(0, 0)=M,(0, y). When he substitutes this in equation 
(22) and solves for D**, Morgan obtains 


C(y) (24) 


_ Bry? M.(0, 0) 


g 


2f 
D*2= D*(0, 0) + 7 * (25) 
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Another relation between y* and D* can be obtained from the Bernoulli 
equation (21), because at the outer edge of the boundary layer v**—0, 


and hence g'D* = B(yr*). (26) 
Thus the function B(*) is determined by 
fy Bx oe 
Bhy(0, y)]=g'(DOO, y)l=9 | Dr, 0)+ 0 a] 5 ey 


When this value of B(7/) is substituted in equation (21) and D* is eliminated 
by equation (25), the following equation is obtained: 


(=) =24'| D0, pote "| 
Ox g 


U 


g 
2fx Br Pe 1/2 . 
«|| D0, O)+ a #4 er a | mv 
-| Dro. 0)+ + Tey |". 


The proper boundary condition is that x=0 be the streamline y*(0, y) =0. 
It is important to note that if £,=0, the quantity o*/ex =0 for all x. Thus, 
no boundary layer forms if the variation of the Coriolis parameter with 
latitude is neglected. Morgan proceeds to analyze the nonlinear differential 
equation (28) by introducing dimensionless variables, Z, 7, y*: 


r=st, y=sy, *=M,(0, 0) sp*; (29) 
and the parameters ¢ and 6 are defined as follows: 


D?(0, 0)g’ 


~ -M,(0, 0) Bs?” oe 
2 
pa fx, (31) 
fxs 
The differential equation (28) thus may be written 
oy* 2 263" 83 —_ = a 
9 * __ 772 
(e = ~ 200, 0)9 1/2 [e+ dyp* + 2yy* — y?] (32) 


x {let op* + prep? — [e+ dp* + 297h* — 97}"9}. 


For numerical integration, the involved numerical coefficient in equation 
(32) can be absorbed into the independent variable by the transformation 


ps8 1/2 a 
(among) ® a 
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Morgan (1956) gives his numerical results in several graphs, from which 
figs. 67 and 68 have been redrawn. Fig. 67 shows the variation of the 
dimensionless transport function 7/* with the dimensionless distance % from 
the coast =0 at y¥=1. The value chosen for /, in all these calculations was 
2 x 10-13/em. sec. The magnitudes chosen for other parameters are given in 
the legend to fig. 67. Except for the fact that the boundary of the southern 


Fig. 67. Dimensionless plot of mass-transport function for several different 
cases, according to Morgan (1956). 
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limit of the interior region of the model indicated by curve 2a is at the 
equator (0=0), curve 2a corresponds to a model very much like the real 
Gulf Stream. Morgan finds the width of the Stream to be 150 km. Curve 2¢ 
represents the same situation with d= 1, thatis, with the southern boundary 
at approximately 15° N. latitude. 
The shift of the southern bound- 
ary produces a slightly narrower 
stream, but the effect is quite 
small. The parameters chosen for 
curve lb represent a homogene- 
ous ocean with a depth equal to 
that of the upper layer of model 
2a. One sees that the stratifica- 
tion does not have much influence 
on the shape of the stream. Com- 
parison may also be made between 
curve 2a and curve la. Curve la 
represents a homogeneous ocean 
the depth of which is like that of 
the real ocean; this is the stream 
which would exist if the ocean 
were homogeneous. It is even 
narrower than the real Stream. 
A further indication of the 
insensibility of the Stream to 
density stratification shows up 
by comparing curves 2a and 2b. 
The density difference in the 
latter is twice that in the former; 0 02 04 06 08 40 
otherwise, all parameters are the 


Fig. 68. Dimensionless plot of the trans- 
same. port per unit width as a function of distance 


Curve 8 shows a very narrow from the western coast, according to Morgan 
stream, the narrowness being dae ee The numbers on the curves refer 
; : o the same choice of parameters as those 
to the increase in s. It applies shown in fig. 67. 
to inertial boundary layers that 
might be formed in connection with the pole-to-pole thermohaline 
circulation discussed in Chapter XI. 

In fig. 68 the dimensionless transport per unit width, d*/0z, is plotted 
against % for each of the cases discussed above. The transport tends to 
decrease monotonically with increasing x because of the decrease in velocity, 
but this is partly offset by the small D* in the two-layer models very close 
to x=0; thus in these cases the transport first increases and then decreases 
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with increasing Z. In this way the streams shown by curves 1b and 2a are 
of about the same width, but the latter is displaced somewhat away from 
the coast. 

Although Morgan’s curves indicate that stratification does not change the 
form of the stream in an inertial boundary layer, there is an important 
point that, as he has pointed out, needs emphasis. First, let us write equa- 
tion (25) in terms of the dimensionless variables: 


M0, 0)br8P gy 


D** = D*(0, 0) ea cae 2(0, 0) f*— 


Thus (0, 0) £;s?/g’ is a measure of the change of depth along a streamline, 
and ¢ is the ratio of the characteristic depth D*(0, 0) to the characteristic 
change in depth. For fixed 7 the depth is smallest at the coast. It vanishes 
when aan _ D**(0, 0)9' 
ae I~ M0, 0) fxs 

Thus Morgan shows that if e< 1, the solution cannot be valid for ¥>e, and 
the value ¥=e might be regarded as a latitude north of which a new regime 
of flow must occur. The region of decay, of meanders and eddies beyond 
Cape Hatteras, comes to mind. 

There are two further remarks that we can make concerning what, in 
nature, this limiting condition (35) may imply. First, it should be remem- 
bered that the warm water of the upper layer is produced in the interior by 
climatic influences acting upon the sea surface in subtropical and equatorial 
regions. Therefore it seems reasonable to suppose that the interface (or 
thermocline) ought to come to the surface in a latitude between the latitude 
of maximum curl of the wind stress and that of the northern limit of the 
subtropical gyre—between 35 and 50°N. latitude. In other words, the 
choice of D(0, 0) is arbitrary only so long as we do not include climatic 
influence in our models. Furthermore, if D(x, y) =0 does not coincide fairly 
well with y(x, y) =0 in the region of maximum westerly wind stress, water 
of the upper layer is not conserved in the upper layer. Of course, the real 
ocean is not a simple two-layer system, and we cannot demand complete 
conservation in reality, but the climatic influences acting upon the sea do 
not appear to be powerful enough to offset an excessive loss of surface water. 
Much of it must recirculate without major density modification. 

Secondly, the quantity D(0, 0) may be determined, to some extent, by the 
control action (in the hydraulic engineering sense) of the Gulf Stream itself. 
As the upper layer is made thinner and thinner, the velocities within it 
become higher and higher, until finally a minimum depth of that layer is 
reached beyond which the transport required by the interior wind-stress 
solution cannot be sustained [equation (35)]. It is at this very stage that 
the crude meander theory outlined in Chapter IX indicates the onset of 


(35) 
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meander instability. One is tempted to speculate, therefore, that the Gulf 
Stream automatically reaches the condition given by equation (35) at the 
latitude of maximum wind curl, and that it acts in much the same way as 
an ordinary open-channel transition (weir, overfall, and so on) to control 
the over-all depth of the thermocline in the North Atlantic Ocean. 

One of the most gratifying features of the inertial theory of the Gulf 
Stream is that it does not contain an arbitrary parameter such as the 
viscous boundary-layer theory does. A disconcerting feature is the absence 
of a well-defined countercurrent to the east of the Stream. Formally this is 
accounted for by the lower-order derivatives in the inertial boundary-layer 
theory. Too much emphasis should not be placed upon this apparent 
advantage of the viscous over the inertial theory, because there is a bit of 
confusion among various investigators over the countercurrent to the east 
of the Stream: there are two ‘countercurrents’ of widely different nature. 
The first kind is apparently associated with the warm core. It is evidently 
a dynamical necessity resulting from advection of warm water from lower 
latitudes (see fig. 33). The simple two-layer model treats the upper layer as 
uniform, so of course it cannot show this feature. Since this countercurrent 
has a width similar to that of the main Gulf Stream, I assume it is the 
countercurrent which the Munk theory attributes to lateral friction. 

The countercurrent of the second kind is a much deeper, broader feature 
of the circulation. It shows up somewhat in fig. 11, as indicated by the close 
spacing of the 800 and 1000 m. contours of the 27-0 o; surface along the 
30° N. latitude circle. It is much too wide to be part of the boundary layer, 
and there is no feature of the wind-stress distribution which can account 
for it in the interior. Perhaps this countercurrent of the second kind reveals 
some fundamental discrepancy between theory and reality; just what it is, 
I do not know. 

Whereas Morgan’s treatment (1956) shows quite clearly the role of the 
various physical parameters in the determination of the Stream profile, it 
is not aimed at giving a detailed representation of the actual Gulf Stream. 
Charney’s (1955) study constructs a two-layer ocean model with a choice of 
D and y on the outer edge of the boundary layer as near to what actually 
occurs in nature as is possible. The two studies thus supplement each other. 

Charney’s formulation makes use of the Bernoulli equation (21) and the 
equation of potential-vorticity conservation, equation (5). If we make use 
of the fact that in the boundary layer, | dv*/dx | > |@u*/dy |, and use the de- 
finition v* D* = dyy*/dx, the potential-vorticity equation is of the following 
form: 


7) (5 ) 
alageraaia 
CLINT Ce: = Fy). (36) 
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The arbitrary function F(y*) is determined by matching functions at the 
outer edge of the boundary layer with observed data. Here we make use of 
the fact that the higher-order terms in equations (21) and (36) vanish at the 
outer edge of the boundary layer. 


Ff _ ryys 
Sak); (37) 
g' D=B(p*). (38) 
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Fig. 69. Comparison of theoretical results and observational material, 
according to Charney’s (1955) analysis. a, calculated depth D* and mass- 
transport function y* of the upper layer; contours of D* (solid lines) are 
drawn for each 100 m., and streamlines (broken lines) for each 10 x 10® m.°/sec. 
b, schematic chart of the observed mean depths of the 10° C. isotherm. The 
scale of the abscissa is exaggerated so that details of the field will show up 
better. 


At the outer edge of the boundary layer D is prescribed by observation as 
a function of y. Since at this place 


epee UPR Bed B) 

Ee I 
we also obtain 7 as a function of y. Therefore, at the outer edge f/D and g'D 
are known functions of y, and hence of y. In this way Charney (1955) 
determines F(y) and B(y) along the outer edge of the boundary layer, and 


therefore at every point inside the boundary layer which is connected with 
the outer boundary by a transport line. 


NONLINEAR THEORIES—INERTIAL 125 


In order to make a close comparison between the results of theory and 
those of observation, Charney has chosen a parabolic form for the transport 
function along the outer edge of the boundary layer, and has introduced 
an arbitrary Florida Straits transport. The depth of the interface in the 
two-layer model is identified with the depth of the 10° C. isothermal surface. 
The results of Charney’s attempt to reproduce the flow pattern in the growth 
region of the Gulf Stream, between the Florida Straits and Cape Hatteras, 
are shown in fig. 69. Fig. 69, a, depicts the theoretically deduced contours of 
the thickness of the upper layer, D*, in intervals of 100 m. (solid lines); and 
also contours of the transport function, y*, at intervals of 10 x 10®m.3/see, 
(broken lines). The x-scale has been exaggerated by a factor of five. Fig. 69, }, 
is a schematization of the depths of the 10°C. isotherm taken from fig. 66. 
The agreement between the two charts is good, and, as Charney says, could 
doubtless be made even more complete by a numerical treatment involving 
actual coastal geometry, and real continuous density stratification along 
the outer edges of the boundary layer. At present such refinements do not 
seem worth while, however. One cannot help believing that a correct zero- 
order approximation has now been achieved. 


Chapter Nine 
MEANDERS IN THE STREAM 


Very little is really known about the meanders of the Gulf Stream. They 
have been observed only a few times; beyond their wavelength, the 
direction of propagation, and the fact that they may grow into detached 
eddies, nothing is known. There are no detailed data on their structure as 
a function of depth, especially below the 900 ft. level commonly reached by 
the bathythermograph. However, because they are essentially wavelike 
they are tempting theoretically, and attempts have been made to treat them 
by linearized perturbation equations. 


THE APPLICATION OF THE PERTURBATION METHOD 


The perturbation method has been eminently successful in the classical 
theories! of sound waves, of tides, of gravitational water waves, and of the 
Bénard-type thermal-convection cells. Most of these problems are peculiarly 
suited to the perturbation method because they involve very simple models: 
the fluid is at rest, or in uniform motion, and the coefficients of the pertur- 
bation equations are constants. Helmholtz (1868) first investigated the 
following unstable system. Two homogeneous currents move, one above the 
other, with uniform (but different) velocities, under the influence of gravity. 
The upper layer is less dense than the lower. The shear at the common 
boundary tends to produce instability; gravity tends to produce stability. 
Short waves on the common boundary are shown to be unstable; long waves 
are stable. 

Kelvin (1871) applied the perturbation technique to the problem of wind 


1 This historical section is after Queney (1950). 
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blowing over water and introduced capillarity as an additional stabilizing 
influence, which acts mostly on very short waves. Asa result, the system is 
stable for all wavelengths of perturbation if the wind is below a certain 
critical speed. With increasing wind speed a single wavelength becomes 
unstable. Kelvin proposed to apply this theory to the formation of wind 
waves, but observation shows that real wind waves are more complicated 
than those envisaged in this simple theory. Even this problem, which at 
first seems so easy and for which there is no lack of observational material, 
has never been successfully resolved. 

Meteorologists must deal with problems of a vastly more complicated 
nature when they attempt to apply the perturbation method to such 
phenomena as the development of waves and cyclones in the atmosphere. 
After J. Bjerknes and his co-workers discovered atmospheric fronts in 1919, 
they came to regard cyclones as unstable waves of the common boundary 
(inclined frontal surface) between two air masses with uniform (but different) 
densities and velocities. Tentative approximate solutions were given by 
V. Bjerknes, J. Bjerknes, Solberg, and Bergeron (1933). Owing to the fact 
that there are two stabilizing factors—gravity, operating now on short 
waves, and the earth’s rotation, acting on the longest waves—the system is 
dynamically stable so long as the wind shear is less than a critical value, but 
when this value is just exceeded, a dominant unstable wave appears. This 
was interpreted as a cyclone. 

The chief difficulty of this model is that it was soon discovered that 
fronts never exist outside of cyclones, a fact which made it seem unlikely 
that they could be the origin of cyclones. Some cyclones have no fronts 
associated with them at all. With the increase of meteorological data it has 
become increasing clear that cyclones are not the result of a degenerating 
mean zonal motion; indeed, Starr (1953) suggests that they are the cause of 
the mean motion. 

Modern meteorologists have abandoned the notion that large-scale 
atmospheric perturbations are necessarily associated with discontinuities 
(fronts), and have turned to other models, such as that of Rossby’s (1947) 
theory of long waves in the westerlies. Other present-day models involve 
vertical wind shears, and horizontal temperature gradients (see Fjortoft, 
1951, for a summary) of a form so remote from the structure of the Gulf 
Stream that they are not likely to be of help in understanding meanders. 

When we come to consider the Gulf Stream, we are not embarrassed by 
the difficulty involved in early cyclone theory: there is no doubt that the 
Gulf Stream is the primary phenomenon, and that the meanders are 
secondary. Hence the Gulf Stream may surely be regarded as the basic 
undisturbed current, on top of which the meanders appear as perturba- 
tions. 
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On the other hand, we have so little consecutive data on the develop- 
ment of a meander into an eddy that we are at a distinct disadvantage. For 
example: after an eddy breaks off, does the parent stream mend itself, or do 
the two ruptured ends remain broken, the severed fragment drifting away 
to the northeast, and the fresh end penetrating to the east to form a new 
extension of the Gulf Stream ? 


GULF STREAM MEANDERS 


The Gulf Stream meanders observed on the Multiple Ship Survey of June, 
1950 (Fuglister and Worthington, 1951), have inspired several theoretical 
studies. It has been suggested that they are analogous to the waves in the 
atmospheric jet stream, but both their small scale (100-400 km.) and their 
small velocities of propagation suggest that the variation of the Coriolis 
parameter with latitude is not a dominant factor in their dynamics. 
Haurwitz and Panofsky (1950) have constructed several models of currents 
in a homogeneous ocean with cross-stream velocity profiles similar to those 
observed in the true Gulf Stream, and have carried out an intricate pertur- 
bation analysis to show the existence of unstable waves with reasonable 
velocities of propagation. These waves are a result of the shearing instability 
of the Stream. 

In aseries of lectures at Woods Hole in 1954 on the stability of large-scale 
motions Dr Jule Charney showed how application of some meteorological 
studies might be made to Gulf Stream meanders, in particular, the work of 
Kuo and Phillips. Kuo (1949) has made an extensive study of the instability 
of single smooth-profiled streams in a barotropic atmosphere (one in which 
the pressure and density surfaces are parallel). The wavelength of maximum 
instability in Kuo’s theory is about 2-7 times the width of the stream, or 
about 140 km. Phillips (1951) has discussed the instability of a particular 
baroclinic model (pressure and density surfaces intersecting) consisting of 
two layers of equal thickness, one flowing over the other. For values of the 
parameters similar to those of the Gulf Stream, the maximum instability 
occurs for waves about 300 km. long, the rate of increase of amplitude, 
about three times in two weeks. Neither of these models resembles the real 
Gulf Stream very closely. Perhaps all we learn from them is that the Gulf 
Stream may be unstable, either barotropically or baroclinically, or both. 

It seemed to me (Stommel, 1953) that certain types of meanders might 
exist in which stratification and inertia are dynamically important. In 
order to illustrate such a system, meanders in a very wide current were 
studied. To simplify the analysis, the realistic cross-stream velocity profiles 
of Haurwitz and Panofsky (1950) were abandoned. The density stratifica- 
tion of the real ocean is approximated by a two-layer system. 
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A SIMPLE MEANDER THEORY FOR A VERY WIDE 
CURRENT IN A STRATIFIED OCEAN 


Let us suppose that the lower layer is very deep, and hence that the hori- 
zontal pressure gradients vanish in it at all times. In the undisturbed state 
a steady current U flows in the x-direction in the upper layer of thickness D. 
Associated with this current is a cross-stream pressure gradient of the 
following form: 

,oD 


fU=-q9 by” 


(1) 
We now suppose that there are small perturbations uw, v, in the velocity 
components, and h, in the elevation of the free surface, and that these 
quantities are independent of y, the cross-stream coérdinate. The pertur- 
bation equations may be written in the form 
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If the perturbations are all in the form e***-”) we obtain the frequency 
equation ik 
eur —pr—| pra?” D|a—p)+P=0. (6) 


where p=c/U=v/kU =p’ + ip", the real part p’ of which is the ratio of the 
velocity of propagation of the wave to the velocity of the current, and the 
imaginary part p” of which gives the instability of the wave motion. For 
the particular range of parameters involved, no one of these terms is small 
compared to the others. It is convenient to rewrite this equation in the 


form pio” Py" (6) 
where 
f2 \'8/f2+9'k2 D 
P=2( scar f? (7) 
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2 \-1/3 . 
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The roots of this equation are all real, provided P>3, in which case there 
are three types of stable wave present. If P <3 there is a region of unstable 
9 SGS 
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waves. These are of most interest to us because they are the only ones 
likely to grow large enough to be noticed on a ship survey. Examination of 
the coefficients of the frequency equation reveals that for U*<g’D all waves 
are stable. At U?=g’D a single wave number given by k=f/,/2 U becomes 
‘just unstable’, whereas all other wave numbers are stable. For slightly 
larger values of U? there is a narrow range of wave numbers about k=f/,/2 U 
in which waves are unstable. In the critical case of marginal stability the 
‘just unstable’ wave is stationary. 

One objection to applying this model to the meanders observed in the 
Gulf Stream is that the real Gulf Stream is not very wide. A more sophisti- 
cated theory would include lateral boundaries to the Stream and would 
provide for resting layers of water on each flank beyond the boundaries. 
Also, a perturbation theory applies only to waves of infinitesimal amplitude, 
whereas meanders often grow to large amplitude. Therefore it is important 
to regard this treatment of meanders as merely indicative of a possible 
mechanism for describing the meandering of a stratified current; but it is 
hard to see how it could be tested by actual observation of the crests and 
troughs of meanders, because the theory and observational techniques are 
both too crude to admit of a meaningful comparison. 

Application to the Gulf Stream.—Rossby (1951) has shown that the 
velocity of the Gulf Stream does in fact approach the critical value ./(g’D). 
This is also true for the stream of uniform potential vorticity (see Chapter 
VIII). We might expect the Stream to become progressively shallower 
downstream as a result of friction, and gradually to approach the critical 
condition. Because of the paucity of hydrographic sections of the Stream 
in any one year or season, it is necessary to construct a composite series of 
sections, in order to determine whether there is any noticeable change in 
the depth of the Stream along its axis. A number of sections all made in 
early June of several years have been assembled and recomputed: (i) one at 
Hatteras at about 74° W.; (ii) two more on the Montauk Point—Bermuda 
line surveyed in Iselin’s (1940) studies; (iii) one at 58° W.; and (iv) two Ice 
Patrol sections along the 50° W. meridian. The geostrophic transports at 
different depths were computed. In order to exhibit any changes in the 
depth of the Stream, the percentage of the total transport below certain 
selected depths was computed. The results are given in table 6; it is clear 
that there is no striking change in the depth of the Stream indicated by 
these sections. The inference to be drawn from this is that from Cape 
Hatteras to the tail of the Grand Banks the Stream is near to the critical 
velocity almost everywhere, and that unstable meanders might be expected 
anywhere. 

We may ask ourselves what the size of the meanders predicted by our 
two-layer meander theory might be expected to be. A surface layer 200 m. 


MEANDERS 131 
TABLE 6 
PERCENTAGE OF TOTAL TRANSPORT BENEATH SELECTED DEPTHS AT 
DIFFERENT JUNE HyDROGRAPHIC SECTIONS ALONG THE GULF STREAM 
(2000 Meter REFERENCE LEVEL) 
74° W. 68° W. 53° W 50° W 
Off Cape Montauk Pt.— SM ae Bvs 
Hattore Bermuda line Section Sections 
Depth in meters Stations 
Atlantis Atlantis Atlantis Ice Patrol 
4565-4570 | 2867—2871 | 3054-3058 | 2616-2620 | 2715-2719 | 4175-4184 
(ees eer. 5 oe, cus 100 100 100 100 100 100 
DOME Se cis cladee 94 92 92 92 91 92 
MOO Re os 3 fsb: 88 85 84 84 82 84 
OO ere sisters < 74 71 69 69 65 69 
OOF aarti ate’. 61 58 55 56 50 56 
ALO eri Ars tonne ce: fa ceve 49 46 43 44 37 45 
U0 wcieataeeeaSne nearer 30 38 23 26 19 28 
OU) teeta Seo eae 112 13 7 9 6 12 
OO Mareen sete es 4 7 2 3 2 4 
GUD: sa fee eae 1 1 0 1 0 1 
DOOOF ce ths me s.cctok ) 0 0 0 0 0 


Source: Stommel (1953, table 2). 


thick, moving at 200 cm./sec. with a density difference of Ap/p=2 x 10-°, is 
critical. The wavelength of the ‘just unstable’ perturbation corresponding 
to this choice of parameters is 180 km. All other wavelengths are stable 
and do not grow. It is interesting that this wavelength corresponds very 
closely to that of the large stationary meander observed to grow and detach 
into an eddy (see Chapter V), but other admissible choices of parameters 
would have led to poorer agreement. 

It is important to emphasize that the meander theory presented here is 
not complete or proven, but merely suggestive of a type of wave motion 
which may possibly dominate the dynamics of meanders. If the energy of 
meanders is absorbed from the potential energy of the cross-stream mass 
distribution rather than from the kinetic energy of the flow, the meanders 
are a part of the thermohaline circulation of the ocean. 

The adjustments which the water masses on either flank of the Gulf 
Stream must make during the passage of a series of meanders should also be 
an interesting subject of study. Recording pressure gauges set on the 
bottom of the ocean, or on sea mounts, might give indications of these 
motions. Unfortunately, the bathythermographic data on hand are not 
sufficient for any statistical study of fluctuations of properties in the slope 
waters and coastal waters. 

There is no certainty, of course, that the growth of the meanders observed 
in the Gulf Stream is necessarily due to an instability of the flow such as 

9-2 
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is conceived in the model of the two-layer meander theory. The meanders 
might conceivably be formed and forced to great amplitudes by the wind 
without the assistance of any instability. In the weather maps for the 
months of May and June, 1950, there is no marked wind pattern that could 
be directly forcing the meander pattern observed in the Multiple Ship 
Survey, and no violent storms over limited areas are noted. 

If we consider the eddy formation observed in the course of the Multiple 
Ship Survey of June, 1950, we note that the Stream is continuous and re- 
latively straight both upstream and downstream of the eddy at 60°W. 
(fig. 40). One wonders whether there was any marked variation in the 
component of surface wind stress along the current at 60°W. From the 
North Atlantic surface wind maps made up at La Guardia Field it is possible 
to obtain the daily ship reports for every 5° along the 40° N. circle of latitude 
from 70 to 50° W. If we then compute five-day means of the surface wind 
stress we find that there was a region of low downstream stress west of 
60° W., and a region of fairly high stress east of 60° W. This pattern per- 
sisted for almost thirty days before the survey of the eddy. The results are 
given in table 7, the selected coefficient of shear stress being 2 x 10-°. 


TABLE 7 


Mean CoMPonent oF Dainty WIND STRESS IN DOWNSTREAM DIRECTION 
FROM May 15, 1950, ro JUNE 19, 1950, on 40° N. LatiruDE 


Ibi@averi nisi Abed o aut mes 70° W. 65° W. 60° W. 55° W. 50° W. 45° W. 


Eastward computed wind ; d , : : ; 
Ce EEE } — 0:05 + 0-01 — 0-069 + 0:30 + 0-20 + 0-40 


& Position of eddy. 


It is natural to wonder whether the coincidence of the eddy with the 
abrupt change in wind stress is simply a chance coincidence, or expresses 
a cause-and-effect relationship between the eddy and a monthly mean 
feature of the wind stress. If the change in wind stress occurs in association 
with a change in mean current velocity (towed-electrode data on the 
Multiple Ship Survey of 1950 indicate that east of 60° W. the current 
velocities fell to about 50 per cent of the velocities to the west), the section 
of the Stream at 60° W. could act as an important block to passage of 
meanders. Meanders that were arrested at such a point in the Stream might 
easily grow into eddies. This is all conjecture, of course; other explanations 
are possible. For example, in Tolstoy’s bathymetric chart (see fig. 16) is 
a series of sea mounts, some as shallow as 1500 m., in the region where the 
Gulf Stream reaches 65 to 60° W. longitude. It is conceivable that these 
sea mounts, and other uncharted ones in the region, are the perturbing 
cause of the breakdown of the Stream in this area. Bottom features are 


MEANDERS 133 


known to exert an influence upon the path of some surface currents, as is 
shown, for example, in the case of the Antarctic Circumpolar Current, by 
Sverdrup, Johnson, and Fleming (1942). 


STABLE MEANDERS IN A STREAM OF 
UNIFORM POTENTIAL VORTICITY 


The following dynamical study treats the problem of stable meanders in a 
Stream of uniform absolute vorticity. 

We shall suppose, as before, that absolute vorticity is conserved in the 
upper layer, and hence that 


Ov 
la 2 
a ae (9) 


0 


= 


where h, is the value of h(x) at 2—>oo. 

We shall also suppose that the Stream is not straight, but slightly curved, 
with a radius of curvature Z. If # is much greater than the width of the 
Stream, and if we assign a positive sign to Z for cyclonic curvature, then 
the dynamical equation of the upper layer is 


fots=7 x (10) 


where g’=g(Ap/p). Eliminate h by cross-differentiation; and in this way 
we obtain 
we eda 


v+— 


7Re waa’ my) 


where c=,/(9’hy), and «=f/c. We regard the nonlinear term as small in 
comparison with the others; hence, we express v in series form: 


v4 on 
= — Sales 12 
v 0+ oR (aap (12) 
Therefore, the functions vp, v;, v2,... are determined by the following 
equations: 
og— aa ee (13) 
,+c w=av,, (14) 
Vgt+C—12u)0, =a 05, (15) 


Vg +¢-1 (2vyv, +02) =a-2 05, (16) 
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the solutions of which are: 


Vy =ce-2*, (17) 
y= gem, (18) 
v= $ (19) 
mace te, (20) 


The first approximation of h may be obtained from equation (10) by direct 


integration: 
= v v 
_ 7 {} (r+ 25) +3 dx+hg 


21 
gee (21) 
=> pap aia Be gest | 
The value of x for which h vanishes is given approximately by 
2 
= 22 
Oe 30K = 


Because the shape of the interface is not strikingly different in troughs and 
crests, it is not likely that it would be very easy to test the hypotheses 
entering into the theory by hydrographic sections across crests and troughs 
of meanders. A more plausible test is to compare the velocity in a crest and 
a trough at the same value of h. For purposes of illustration we construct a 


TABLE 8 


NUMERICAL EXAMPLE OF VELOCITY PROFILE IN MEANDERS 


h/ho vie (cyclonic)® v/e (anticyclonic)® 


Or Oeearoreroreis kc ociarers autres) sietouele se 0-958 1-060 
Cee netare nee ceest etree citattetars viele tere cehe x 0-862 0-950 
OrD Seavert SelGcutalede: sfesa’aistovegenenayeshs 0-772 0-829 
UERecob Molec onload Gmope oO Oe 0-680 0-720 
Yaar oc bonis wl sos eae © vieeiee 6) as 0-590 0-620 
(TAD Wp. O Din gd crore Und olor Ogsaoa ae 0-495 0-510 
OrGitsAtheie ote pico Sine omen oiaioro noe 0-395 0-408 
STR arte Pereiaatcua’s arrieie a alas 0-300 0-300 
Sipser kev ateamepeiel ctstch ote snr ysy tle ory < 0-200 | 0-200 
De Oia sane erene. swonensie she Dyevete a tases) aor 0-100 0-100 
We apararerevenet enone ave tere herent ont rene fa cote 0-000 0-000 


® Probable error, + 0-002. : 
> For h/hy values of 0:7 and higher, the values of v/e (cyclonic) and v/e (anticyclonic) are 
essentially equal, within graphical error. 
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table of h and v for a meander with radius of curvature Z=(2/3) x 10? em. 
and «= 10-7/em. The values of v are then plotted against h for both cyclonic 
and anticyclonic curvature. The results, scaled from the plot, are shown 
in table 8. The difference is most pronounced for small values of h/hy. The 
proper place for a test is probably at about h/h)=0-2. 

Since the total transport of the Stream in a crest is greater than that 
through a trough, the meander pattern moves, as a whole, downstream. The 
rate of this progression may be found by the following simple reasoning. 
First evaluate the total transport 7’: 


fl -| vhdx 
x=2/(3a?Z) 


ae 1 i 
=e Ol or 6a2Z |" 


When the values of parameters introduced above are used, the transport 
at the crest is greater than that through the trough by about 10 per cent 
of the average transport. Suppose that the shape of the meander is 


(23) 


x=2, cos (ky — vt). (24) 


At time t=0, the rate of increase of volume W of warm water between 
ky=0 and ky=7 due to the motion of the meander is 2x)h,(v/k). This, of 
course, must be equal to the excess of transport through a crest (ky=0) 
over that through a trough (ky=7), chy)/3a?Z. But at t=0 and y=0, the 
radius of curvature & is also given in terms of the equation 


fa 
FA= aan tol (25) 
hence 
vy ck 
== Gf [ (26) 


The rate of progress of a meander with a wavelength of about 300 km. 
is therefore about 5 cm./sec. The actual movement of meanders ought to 
be studied extensively by means of the air-borne radiation thermometer. 


Chapter Ten 


FLUCTUATIONS IN THE 
CURRENTS 


Many catastrophes of an economic kind, such as the failure of the rice crop 
in Japan, or of a certain fishery, or years of unusual numbers of icebergs in 
shipping lanes, are attributed to fluctuations in ocean currents. Very little 
is really known about such fluctuations. It takes years of careful and 
expensive observation to produce even a very crude description of them. 
The scientific programs of our oceanographic institutions are not geared to 
long-term problems of this kind; there is much pressure for novelty, much 
temptation to follow the latest fad, and a persistent though erroneous 
notion that all worth-while problems will eventually be solved by some 
simple, ingenious idea or clever gadget. A well-planned long-term survey 
designed to reveal fluctuations in ocean currents would be expensive and 
time-consuming. It might even fail, because of inadequacies of the tools 
we have at hand. But until this burdensome and not immediately reward- 
ing task is undertaken, our information about the fluctuations of ocean 
currents will always be fragmentary. 


THE DYNAMICS OF THE FLORIDA CURRENT 


Numerous studies of tide-gauge data have been made (Hela, 1952; Mont- 
gomery, 19386), to obtain information about seasonal fluctuations in the 
cross-stream slope of the sea surface; and then, by the geostrophic equa- 
tion, to compute seasonal fluctuations of the mean surface velocity of the 
Florida Current. In addition to the tide-gauge data already published 
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concerning Key West, Miami, and Cat Key, there are now available, for the 
first time, tide-gauge data for Havana. This information is obtained from 
a gauge set up by the United States Coast and Geodetic Survey in 1946. 
A leveling survey has been made between Key West and Miami (Mont- 
gomery, 1941a), but it is, of course, impossible to connect Havana or Cat 
Key to the same datum by ordinary means (Montgomery, 1947). Therefore 
we cannot obtain true differences in sea level across the Straits, but we 
can obtain fluctuations in the differences of sea level at two stations, as 
shown in fig. 70. At all four stations the maximum sea level occurs in 


O——O HAVANA—KEY WEST 
O4—A CAT KEY—MIAMI 
X——X HAVANA—CAT KEY 
O——O KEY WEST—MIAMI 


J F M A M J J A Ss te) N 19) 


Fig. 70. Annual changes in the difference of sea level between various pairs 
of stations on the Florida Straits. From Stommel (1953, fig. 1). 


September and October, apparently caused by the summer heating of the 
water. This average rise does not appear in the differences. 

One interesting feature of these differences is that the fluctuation in the 
Miami-—Cat Key differences is about twice that of the Key West—Havana 
differences. The maximum cross-stream slopes occur during July, when the 
flow is strongest (Hela, 1952), and the maximum downstream slope from 
Key West to Miami also occurs at this time. These relations of slopes to 
surface velocity are quite what might be expected from the most elemen- 
tary considerations of Bernoulli’s equation and the geostrophic equation 
(see Chapter III). But the fact that the fluctuations of the differences at 
the Miami—Cat Key section are larger than at the Key West—Havana 
section, and the fact that the downstream slope between Havana and Cat 
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Key is small in July, rather than largest, are at first sight rather startling 
and require explanation. 

Another noteworthy feature of the Florida Current is that the region of 
anticyclonic shear in the section off Miami is extensive. So far as direct 
velocity measurements from an anchored vessel are concerned, the data 


CM/S 


150 


50 


(0) 
MIAMI CAT KEY 
Fig. 71. Sample surface velocity profile of the Florida Current between 


Miami and Cat Key, according to Murray’s (1952) towed-electrode measure- 
ments. 


are scanty. The only observations at different depths are those made by 
Pillsbury (1891), and these show a rapid decrease of velocity with depth. 
The axis of maximum surface velocity is not in the center, but is dis- 
placed toward Miami; and this feature has been verified on many crossings 
of the Current off Miami by Murray (1952). Fig. 71 shows measurements 
made by towed electrodes on one of these crossings, but it should be em- 
phasized that the towed-electrode method is likely to give particularly 
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misleading readings of velocity in the Florida Straits. Nevertheless, the 
general fact is that a wide zone of anticyclonic vorticity (of approximately 
—[0-5+0-1]f, where f is the local Coriolis parameter) seems to be 
established and also requires explanation. It is, of course, impossible for 
the variation of the Coriolis parameter with latitude to cause such a shear 
over so short a distance as the length of the Florida Straits. 

At Miami the channel is only about one-half as deep and one-half as 
wide as at Key West. The change in depth can have no important in. 


CAT KEY 


HAVANA 


Fig. 72. Hypothetical position of the interface in the Florida Straits. The 
arrows indicate the magnitude of the current. In the actual Florida Straits 
there is a right-angle bend in the channel between Key West—Havana and 
Miami—Cat Key. Bottom topography is not indicated. From Stommel (1953, 
fig. 4). 


fluence on the flow, since the bottom water does not have an appreciable 
velocity at either section. The narrowing of the channel at Miami is very 
important hydrographically, for in order to pass through it the water is 
accelerated, and this requires a small drop in the level of the free surface 
from Key West to Miami. Since the lower layers are at rest (except for 
tides), the isopycnic surfaces must slope upward toward Miami to counter- 
act the axial pressure gradient in the upper accelerating layers. Because 
this slope must be some 500 times the drop of the free surface, it produces 
a marked decrease in the thickness of the surface layers as Miami is 
approached. By conservation of potential vorticity, this vertical shrinking 
produces an anticyclonic shear and magnifies the transverse geostrophic 
slope of the free surface. A crude quantitative analysis based on a 
two-layer model is presented below. Fig. 72 illustrates the average 
hydrographic structure of the Florida Current envisaged in this hypo- 
thesis. 
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Consider a layer of water of density p,, depth D, overlying a deep, 
resting layer of density p,. The upper layer has a velocity U at the Key 
West—Havana section («=0), and its vorticity is zero. Farther down- 
stream (x=Az), at the narrower Miami—Cat Key section, the Stream has 
accelerated by a drop in the free surface Ah. Because the lower layer is at 
rest, the interface must slope upward to offset the effect of the free surface. 
If the axial pressure gradient vanishes in the lower layer, then the following 
relation must hold: 


(2) OD _ ch (1) 


Pp, |} Ox Ou 


The linearized vorticity equation in this simple case may be written as 
follows, where ¢ is the vorticity of the upper layer: 


0 (E+f\ _ 
05 (*5)=0. (2) 


We assume that the vorticity vanishes at Key West; hence the vorticity 
at Miami is given approximately by 


Ah pr» 
eae D p.—py 


(3) 


The actual acceleration observed at Miami corresponds to a value of Ah of 
about —20 cm. (This is somewhat in excess of the —4-9 cm. obtained by 
leveling between Key West and Miami [Montgomery, 194la].) The 
density difference of the two layers is approximately (p,—/,)/p.= 2 x 10-*. 
The initial depth D of the current is roughly 250 m.; and this yields a value 
of the vorticity increase between Key West and Miami of A¢= —0-4 f, the 
observed value given above. 

If the high anticyclonic shear at Miami is actually due to the vertical 
shrinking of the upper layers, it should be possible, by careful hydro- 
graphic observations on the two sections, to try to detect the average rise 
of about 100 m. in the isopyenic surface at a depth of 100-300 m. It would 
be difficult and tedious to measure this hypothetical rise, because the rise 
may be partly masked by the large cross-stream slope of the isopyenic 
surfaces, and by variations in depth caused by the tides. Present data 
are insufficient to form a basis for a satisfactory test. 

The Miami tide gauge is most sensitive to fluctuations in the transport 
through the Florida Straits, because it lies at the western end of the 
shallowest section of the Stream. 
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IRREGULAR FLUCTUATIONS IN TRANSPORT OF 
THE FLORIDA CURRENT 


The most spectacular fluctuations of the Florida Current are to be seen in 
the periods of extraordinary flow revealed by the submarine-cable measure- 
ments throughout the months of December, 1952, and January and April, 
1953 (fig. 73). During these months the Florida Current maintained a flow 
of about 8 x 10® m./sec. in excess of its normal flow. Mr Jerome Namias 


1952 1953 


Fig. 73. Comparison of high-flow regimes in the Florida Current, as obtained 
from the submarine-cable measurements (curve A), with those computed from 
the five-day mean zonal pressure difference between 30 and 20° N. (curve B). 
The upper curve (A) shows excess of flow above mean flow. 


suggested to me that these periods of excessive flow might be associated 
with thirty-day persistent breakdowns of the North Atlantic wind circula- 
tion. After a number of trial analyses of weather maps, I decided to use 
the five-day mean pressure difference between 30 and 20° N., averaged all 
across the ocean, as an index of the strength of the trades north of the 
West Indies. These five-day mean zonal pressure differences (30° minus 20°) 
are plotted in fig. 73. 

It is clear that there are three sustained periods of breakdown of the 
Bermuda—Azores high, and that invariably they precede the periods of 
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high flow of the Florida Current by about a month. It is unlikely that the 
atmospheric-pressure disturbance is the direct cause of the anomalous 
flow, for two reasons. First, an adjustment to a pressure change would 
occur with the speed of a long gravitational surface wave—that is, in 
several hours, not several weeks. Secondly, the mass flux to compensate 
the surface-pressure drop would be of the order of 3 x 10" m.3, whereas the 
total excess flow of the Florida Current observed during each period is of 
the order of 3 x 10% m.°. It seems more likely that the anomalous currents 
are caused by changes in wind-stress distribution over the ocean. But even 
this influence is not direct, because the anomalous wind conditions asso- 
ciate low winds with high discharge of the Florida Current. However, if it 
is remembered that, because of the immense mass of water involved, the 
dynamic topography of the North Atlantic subtropical gyre is constrained 
to remain relatively unaffected by short-period (even thirty-day) wind 
changes, it is possible to offer a rationalization for this inverse relation in 
the following manner: All along the Lesser Antilles, at the entrance to the 
Caribbean Sea, the dynamic topography is such as to produce southwest- 
flowing currents in the upper 400 m. In the very surface layers the normal 
trades tend to produce a net northward drift (Montgomery, 1936a). The 
combined effect of these two tendencies is to produce an essentially west- 
ward flow of the upper 100 m. at 25° N., and a southerly component below 
this depth. The fact that Sargasso Sea water normally enters the Caribbean 
only below 100 m. and not at the surface (Parr, 1939a) lends support to 
the explanation just outlined. 

When the trades cease to blow at 25° N., however, the northward Ekman 
drift ceases abruptly (within 18 hr.), and even the surface waters move 
southwestward. Thus the flow into the Caribbean is increased. An order- 
of-magnitude estimate of this increase may be obtained from the following 
simple relation. In the case of a steady wind stress 7, acting in the 
x-direction (toward the east), the northward Ekman transport r/,, (see 
Ekman, 1905) is given by the relation 


rM,,.= —-— (4) 


where f is the Coriolis parameter, and r is the width of the ocean. Since T, 
is normally very nearly —1 dyne/cm.? at 25° N. latitude, f= 0-6 x 10-4/sec., 
and r=5x108cm., the northward Ekman transport across the 25° N. 
parallel is 8-3 x 10° m.3/sec. When the trades cease, this flux is added to 
the water entering the Caribbean. There is a distance of about 3000 km. 
through the Caribbean to be traversed before this increase reaches the 
Florida Straits. There is therefore the possibility of alag of about fifteen days 
before the first effect of the change in flow appears in the Florida Current. 
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A good check on these ideas could be obtained by means of potential 
measurements on a middle-Caribbean cable such as the one between 
Curacao and Santo Domingo. The time lag should be about half that 
observed in the Florida Straits. 

It is important to be aware that these remarks about causes of fluctua- 
tions of the Florida Current are only speculations and ideas, not proven 
theories in the usual physical sense. 

That such facile speculations may be treacherous is illustrated by certain 
occurrences that have been observed. In January, 1955, there was an 
exceptionally intense and widespread month-long breakdown of the 
Bermuda—Azores atmospheric high-pressure cell, and the northern fringes 
of the trades disappeared. This meteorological event was accompanied by 
a decrease in the transport of water through the Florida Straits—a singular 
instance of complete reversal of the simple explanation discussed above. 


SEASONAL FLUCTUATIONS OF THE GULF STREAM 
NORTHEAST OF CAPE HATTERAS 


In the previous chapters in which we discussed the details of single cross 
sections of the Gulf Stream we were on fairly solid observational ground. 
In Chapter IX, on meanders, the dangers of picking one of many different 
possible interpretations of the data, and ignoring the other interpretations, 
arose. There was an uncomfortable feeling of engaging more with a 
creation of our own minds than with a fact of nature. As we now proceed 
to discuss the fluctuations of the Gulf Stream, we are even more apt to find 
ourselves pursuing a will-o’-the-wisp. Although we do know something 
about seasonal changes in the shallow wind-drift currents in the upper 
100 m. of the ocean, and although we do have fairly adequate data, from 
tide gauges and the submarine cable, for a study of the Florida Current, we 
know virtually nothing about changes in the transport of the Gulf Stream 
itself. With these words of warning, we may now examine what little 
material there is at hand. 

There has never been a cruise which failed to find the Gulf Stream between 
the United States and Bermuda. It is quite clearly a permanent feature. 
Butare the position and transport of the Gulf Stream entirely independent of 
the seasonal fluctuations in the strength and position of the wind system ? 

The pilot charts which are based on hundreds of ship reports on file 
with the United States Navy Hydrographic Office indicate the presence of 
seasonal changes in the surface currents. Fuglister (1951a) has prepared 
a summary of these surface-current data in miles per day. They are pre- 
sented in fig. 74. The curves for the Florida Current off Miami, the Florida 
Current southwest of Cape Hatteras, and the Gulf Stream northeast of 
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Fig. 74. Fuglister’s plots of annual variation of surface currents in various 
areas of the North Atlantic. From Fuglister (1951a, fig. 2). 
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Hatteras, all show a pronounced maximum in early summer. The broken 
lines represent combined annual and semiannual components obtained by 
harmonic analysis. 

Two other types of data have been advanced to support the idea that 
there is a seasonal fluctuation in the Gulf Stream. Iselin (1940) computed 
the volume transport of the Gulf Stream from thirteen hydrographic 
sections made in the period 1937-1940, and believed that these showed an 
annual variation of about 15 x 108 m.3/sec. in the transport of the Stream. 
These data are plotted in heavy black dots in fig. 75. 
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Fig. 76. Annual variations of differences of sea level reduced to normal 
atmospheric pressure. The solid line is used for the Bermuda—Charleston 
section, the broken line for that from Key West to Miami. From Montgomery 
(19386, fig. 71). 


The monthly mean sea level from tide-gauge records should also be an 
indication of changes in transport. 

The geostrophic equation requires that the transverse slope of the free 
surface and surface velocity be proportional. A decrease of mean sea level 
at a station along the United States coast suggests an increase in Gulf 
Stream transport, but does not prove it. Montgomery (19386) has pointed 
out a number of corrections (such as allowance for seasonal variations in 
mean atmospheric pressure) which must be applied to the tide-gauge data. 

Records from Miami and Charleston are also plotted in fig. 75. Miami 
and Charleston are on the same side of the Stream, and hence the records 
from these two places do not, strictly speaking, serve as an indication of 
changes in the cross-stream slope. Montgomery (19386) has plotted a 
curve for changes in the cross-stream slope between Bermuda and 
Charleston, reproduced here as fig. 76. 
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The chief difficulties in utilizing tide-gauge data are: (i) it has been 
impossible, thus far, to make the proper allowance for the slope of the 
water surface produced by the direct stress of the local wind over the 
continental shelf; (ii) the cross-stream slope is related only to surface 
velocity, and, unless something definite is known about the vertical dis- 
tribution of velocity, does not give an unambiguous value of the total 
transport. 

One of the most striking features of fig. 75 is the small amount of 
hydrographic data available, and, worse yet, the concentration of the data 
toward the latter half of each year. The project was interrupted by the 
Second World War and was never resumed. The broken-line curve is 
drawn as a kind of fanciful interpolation, and the ordinate scales of the 
tide-gauge data are chosen simply to make the curves look alike. I think 
one is well advised, when confronted by a set of curves such as those in 
fig. 75, to regard them as merely suggestive, rather than to imagine that 
they prove anything. 


THEORETICAL STUDIES OF TRANSIENT CURRENTS 


Several attempts have been made to construct theoretical models that will 
help in understanding the transient state of the ocean circulation brought 
about by forced changes in the winds. For example, Ichiye (1951) has 
obtained an interesting solution for a wind-driven circulation produced by 
a periodic wind stress acting upon a homogeneous ocean within a rect- 
angular ocean. He uses the frictional boundary layer for the Gulf Stream. 
He finds that the dynamic topography on the western coast consists of a 
series of ‘waves’ of dynamic height moving toward the shore with in- 
creasing amplitude at the rate of one crest per year. Thus the annual 
fluctuation of transport of the Gulf Stream is associated with a shift in 
the Stream’s position. Veronis and Morgan (1955) have obtained a result 
very different from Ichiye’s. The difference apparently arises from their 
inclusion of the important time-variable term in the mass-continuity 
equation (which Ichiye takes as vanishing). Veronis and Morgan treat a 
homogeneous ocean, and find: that the mass transport of the Gulf Stream 
responds to variable wind systems of more than three months’ duration 
with lags of the order of nine days; and that there is no noticeable shift in 
the position of the Stream throughout the year. 

Both these theoretical studies treat the ocean as though it were homo- 
geneous, whereas, of course, the real ocean is stratified in layers of different 
density. In the subtropical regions of the Atlantic Ocean, most of the 
density change with depth occurs in the thermocline, as we have already 
seen (Chapter IV), and we get some approximation to the true density 
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stratification of the ocean by considering a two-layer ocean model. We 
must now ask ourselves, what influence does the density stratification of 
the ocean have upon its response to transient wind systems? Even very 
crude physical arguments seem to suggest that the density stratification 
of the ocean cannot respond rapidly to changes in the wind distribution. 
Examples of such arguments are as follows: 

a) There is an immense store of available potential energy in the deep 
warm-water mass in the Sargasso Sea, more than a thousand times the 
kinetic energy of all the currents of the North Atlantic. We may suppose 
that, on the average, the rate of dissipation of the kinetic energy of large- 
scale currents in the ocean is not more than the rate of work done by the 
stress of the wind on the ocean. (Much of the work done by the wind may 
be dissipated in waves or in vertical shear of the surface currents.) If the 
trade winds exert a stress of 1 dyne/cm.? on the sea surface, and the 
surface current there is roughly 20 cm./sec., the rate of work done on the 
ocean is roughly 20 ergs/cm.?/sec., and this sets a maximum average rate 
to dissipation of kinetic energy per unit area of the ocean currents by tur- 
bulence. We next compute the potential energy stored in the warm central 
core of the North Atlantic Ocean in excess of that which would exist were 
the thermocline level everywhere. In the latter case the thermocline 
would be at a uniform depth of about 550 m. In reality, the thermocline 
is at a depth of 900 m. in half the ocean, and 200 m. in the remainder. If 
the difference in density of the layers above and below the thermocline is 
taken as 3 x 10-* g./cm.3, the potential energy stored in the warm core is 
3 x 10° ergs/cm.?. The potential energy stored is thus sufficient to maintain 
the current system for at least 1700 days in the absence of the driving 
stress of the winds. The average kinetic energy of the North Atlantic 
Ocean is of the order of 4.x 106 ergs/em.?. It is no wonder, then, that the 
density structure of the Gulf Stream System is so uniform and constant 
a feature regardless of the vagaries of the weather. 

b) An alternative order-of-magnitude computation of the time constant 
for the decay of the wind-driven circulation following a complete cessation 
of the wind may be based on the time required for the volume of warm 
water in the Central North Atlantic Ocean (3000 km. x 6500 km. x 0-5 km.) 
to drain off at the full Gulf Stream rate of discharge (70 x 10° m.*/sec.), 
supposing all the Gulf Stream water to be lost to the Arctic seas. This 
period turns out to be 1600 days. 

c) Another estimate of the time constant is provided by a comparison 
of the total volume of warm water with the rate of convergence of surface 
- water due to the anticyclonic wind system: about 1000 days. 

These qualitative considerations suggest that the density structure of 
the ocean (and hence the current field as usually computed from hydro- 
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graphic-station data) does not respond completely to those fluctuations in 
the wind field which have periods of much less than a year. Since such 
arguments as the ones given under a-c might very well be misleading, 
Dr George Veronis and I (1956) have studied in detail the response of a 
two-density ocean to a variable applied wind stress. The model studied 
does not have boundaries, but the applied wind stress is taken as periodic 
in space. Friction does not appear to be an important factor for periods 
of a year or less, and inertial-gravitational wave motions are negligible for 
periods of more than several days. In order that the motion for a very long 
period shall approach the motion represented in the steady-state solutions 
by Sverdrup, Reid, and Munk (see Chapter VII), it is necessary to include 
f, the variation of Coriolis parameter with latitude. For very long periods 
(thirty years or more), the theoretical currents induced by the wind are 
confined entirely to the upper layer, as is very nearly true of wind-driven 
currents of the real ocean, and this surface intensification of currents pro- 
duced by a transient wind system does not occur without /. 

The essential point to make clear is that for periods of between a week 
and a year the density field in the ocean at mid-latitudes does not attain 
a full equilibrium response to the variable wind-stress pattern at any 
instant; but it does respond. Thus, in the North Atlantic, we cannot expect 
the depth of the main thermocline to adjust itself completely to the 
seasonal progress of the Bermuda—Azores anticyclonic wind system ; on the 
other hand, the main thermocline cannot be wholly insensitive to the 
seasonal wind changes, although it would be essentially insensitive to 
wind-stress changes with periods less than a week. The currents induced by 
wind-stress patterns with periods between a week and a year are not 
confined to either layer, but the vertically integrated transport is very 
nearly in equilibrium with the applied wind stress at any instant. Thus, so 
far as vertically integrated transport is concerned, the results obtained by 
Veronis and Morgan (1955) using a homogeneous model and a viscous 
Gulf Stream should be applicable even for a density-stratified ocean. How- 
ever, the wind-induced variations of oceanic circulation should not be 
accurately determinable by the traditional method of geostrophic computa- 
tion from the density field; and the currents are so small that they cannot 
be observed directly by current meter; and the horizontal displacements 
involved are too small to produce noticeable changes in the distribution of 
deep water masses. 

Apparently, in central oceanic areas (away from the equator) the chief 
directly observable effect of the purely wind-induced seasonal variation of 
oceanic circulation is the change in mean sea level, but even this is nearly 
obscured by seasonal changes associated with changes of local density 
arising from local heating, salinity variations, and so on. The island of 
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Bermuda is favorably placed near the western side of the North Atlantic 
Ocean, close to the region of maximum mean wind curl, but well away 
from the Gulf Stream. One should expect to find changes of sea level 
there which correspond closely to those described in the theory of Veronis 
and Morgan. 

Dr June Pattullo, of the Scripps Institution of Oceanography, has been 
engaged for several years in making a study of seasonal changes in mean 
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he 
Fig. 77. An analysis of annual variations of sea level at Bermuda. The four 
curves are described in the text. 


sea level throughout the world, and although the results of this valuable 
project have not yet been published,! she very kindly communicated to me 
the results of her analysis of the Bermuda sea level. Curves 1 and 2 of 
fig. 77 are plotted from Dr Pattullo’s computations. Curve 1 presents the 
monthly mean sea level computed on the basis of hydrographic-station 
data alone. It represents the change in sea level due to changes of local 
density structure and, according to theory, is not directly related to the 


1 Subsequently published by June Pattullo, Walter Munk, Roger Revelle, and 
Elizabeth Strong, in the Journal of Marine Research, 14 (1955):88-156, under the 
title, ‘The Seasonal Oscillation in Sea Level’. 
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wind-induced circulation. Curve 2 gives the actual tide-gauge data 
averaged for seventeen years and corrected for atmospheric pressure. The 
difference between curve 1 and curve 2 is plotted as curve 3. 

It is natural to inquire whether the difference (curve 3) can be due to the 
directly wind-induced changes of oceanic circulation. In order to make a 
rough, tentative order-of-magnitude calculation, I computed the mean 
monthly curl of the square of the zonal winds centered at 32°5 N. latitude 
across the Atlantic. This quantity should be directly proportional to the 
mean monthly curl of the wind stress in the latitude of Bermuda. Ac- 
cording to the theory of Veronis and Morgan (1955), the meridional 
transport in the central oceanic areas should be given by the relation 


__mean monthly transport _ mean monthly curl of squared zonal winds (5) 
~ mean yearly transport —— mean yearly curl of squared zonal winds ° 


We shall call this quantity NV. I was quite surprised to note that N varies 
from 0-6 in September and October to a maximum of 1-6 in December. 
This suggests very much larger variations in oceanic transport than those 
indicated by Iselin (1940) on the basis of hydrographic-station data. It is 
important to remember that the changes of transport which we are com- 
puting cannot appear in hydrographic data, since they are associated with 
the barotropic mode. The variation d in sea level due to these monthly 
variations in transport is given by the formula 


d= fl (N—1), (6) 
g 

where 7’ is the mean yearly transport in the central oceanic area of the 
North Atlantic at 32°5 N. latitude, and D is the total mean depth of the 
ocean. In plotting curve 4 I assumed that f=0-7 x 10-4/sec., 7’ =60 x 
10° m.3/sec., and D=4 km. The agreement of this very crude calculation 
with curve 3 is, to my mind, rather good, and suggests a verification of 
the theory of transient currents. However, a truly careful discussion of the 
problem must await the detailed results of Dr Pattullo’s studies of the 
ocean-wide distribution of changes in sea level. There are other physical 
obstacles to a convincing demonstration. For example, we do not know 
the ways in which the wind stress depends upon the air-sea temperature 
difference, and we do not have an adequate theory of the thermohaline 
circulation in central oceanic areas to use as a companion piece along with 
the theory of wind-driven circulation. These speculations about the data 
and about the applicability of the theoretical models will remain incon- 

clusive until really adequate data are available. 
Systematic long-duration series of oceanographic measurements are not 
obtained on a routine basis. The Western Union cable measurements are a 


152 FLUCTUATIONS 


solitary example of what is needed. In the future, after unattended 
instruments on buoys have been maintained for many years at various 
points in the oceans and as soon as long series of subsurface temperature, 
salinity, and velocity measurements have accumulated, it will be possible 
to discuss these matters more satisfactorily. 


LONG-PERIOD FLUCTUATIONS IN THE GULF STREAM 


There are not sufficient data available to permit us to discuss fluctuations 
in the Gulf Stream System lasting longer than one year. This is because 
early measurements were made without a knowledge of the existence of 
meanders and eddies. Unless allowance can be made for these short- 
period fluctuations, the long-period variations obtained by plotting the 
data are not statistically significant. 

The Kuroshio, which flows along the coast of Japan, and which appears 
superficially to be similar to the Gulf Stream, has exhibited a very remark- 
able change in the period 1919-1950. From 1919 to 1934 (Uda, 1938) the 
Kuroshio had a pattern of flow past Cape Shiono Misaki very much like 
that of the Gulf Stream off Cape Hatteras. Between 1935 and 1942 a large, 
semipermanent cold eddy, 200 km. in diameter, developed at the Cape, and 
the Kuroshio weakened and went completely around this cold mass. Since 
1942 the eddy has slowly decayed and the stream has returned to its 
original state. We have never observed anything like this in the Gulf 
Stream. 


Chapter Eleven 


ROLE OF THE THERMOHALINE 
CIRCULATION 


Although much effort has been directed toward detailed theoretical 
elucidation of the consequences of the action of wind stress in producing 
the oceanic circulation, there has not been a parallel development of the 
theory of the thermohaline process. The distribution of wind stress over the 
ocean surface is much better known than the distribution of net heating 
and cooling of the ocean; thus there is little hope, at present, of constructing 
as complete and satisfactory a theory of the thermohaline circulation as is 
possible for the wind-driven circulation. Nevertheless, it does seem desir- 
able to reconnoiter the problem and to try to construct, on a tentative 
status, a simple theory of the thermohaline circulation that can be applied 
to real ocean circulations. 

In this chapter, attention is first drawn to the fact that the now-familiar 
‘wind-curl’ equation introduced by Sverdrup (1947) is simply a statement 
of the fact that at every geographical position in the ocean the integrated 
horizontal divergence of the Ekman wind drift is compensated by a cor- 
responding convergence of the integrated geostrophic current. Next, a 
thermohaline process is specified in terms of a prescribed vertical velocity 
structure. At some particular mid-depth the vertical velocity has an 
extreme value; and this same depth is the level of no horizontal divergence. 
Two vorticity equations are formed by vertical integration over the layers 
above and below this level. It is shown that the thermohaline circulation 
appears as a kind of internal mode of motion, and that the level of no 
horizontal divergence is also the level of no meridional motion and, if 
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there is a finite interval of no horizontal divergence, this is a level of no 
motion. 

By applying these expressions to the central regions of the oceans, and 
using the principle of mass conservation, expressions are obtained for the 
integrated wind-produced and thermohaline-produced transports in the 
western boundary currents. (It is not necessary to consider the detailed 
structure of these currents, be they viscous or inertial.) A schematic model 
of a rectangular basin is examined, and an attempt is made to relate the 
theory in a more schematic form to actual measurements in the Atlantic 
Ocean. Several predictions (or more properly, deductions) are made: for 
example, that there is a great thermohaline countercurrent under the Gulf 
Stream. It is intended that these predictions shall stimulate direct measure- 
ments of deep ocean currents. The well-known discrepancy between the 
transports of the Gulf Stream and Kuroshio computed from the wind theory 
(Munk, 1950) and those determined from observation (dynamic computa- 
tion) is resolved. It is important to note that the Equatorial Current trans- 
ports in the theories of Sverdrup (1947) and Reid (19485) do not exhibit this 
discrepancy. Finally, some of the limitations and perplexities of the 
rectangular model are discussed. 


THE VORTICITY EQUATION FOR THE CENTRAL 
OCEAN, WIND STRESS ONLY 


The physical meaning of the ‘curl’ equation in the theory of the wind- 
driven oceanic circulation, as first propounded by Sverdrup (1947) and 
further developed by Munk (1950), can be most easily appreciated when 
explained in the following way. We write the linearized stationary equations 
of motion in the form 


Op OTe 

ore a ees 
ee (1) 

eee 

Spu= By Oe" 


Horizontal viscous stresses, inertial terms, and time dependence are neg- 
lected. The x-axis is positive toward the east, the y-axis is positive toward 
the north, and the z-axis is positive upward. The velocity components wu 
and v, the density p, the hydrostatic pressure p, and the vertical shearing- 
stress components 7, and 7, are regarded as functions of x, y, and z. The 
Coriolis parameter f is regarded as a function of y alone. 

The equations (1) are now integrated from the surface z =z, to some great 
constant depth z= —h, at which one supposes that the horizontal pressure 
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gradients and motions vanish. The components of mass transport per unit 
width are defined in the following manner: 


Zo Zo 
M,=| pudz, M,=| pvdz, (2) 
—h —h 


and a function P is introduced: 


2= | " pde. (3) 
—h 
The equations (1) in the integrated form are as follows: 
oP 
ne rere ae ie (4) 
oP sil 
io Tae ae 


The interchange of integration and differentiation in the pressure terms does 
not lead to any additional terms, because the lower limit is chosen at zero 
horizontal pressure gradients, and the terms introduced by the variation 
in the surface elevation (upper limit), —p(z)) (0z)/0x), and so forth, are 
demonstrably negligible (Munk, 1950). Cross-differentiation of these two 
equations, and use of the fact that the horizontal divergence of the in- 
tegrated mass transport vanishes, 


OMy  OMy _ 
Ge du. 


0, 


leads to the ‘curl’ or vorticity equation now so familiar in theories of wind- 


driven ocean circulation: aie 
fM,=curl 7 , (5) 


where 7 is the stress of the wind at the sea surface, and {=0f/dy. 

In asense, the simple mathematical manipulation of cross-differentiation 
tends to obscure the physical meaning of equation (5). Let us therefore 
begin again with the integrated equations (4) and break each mass-transport 
component into two parts: one set to represent the Ekman wind-drift 
transport components M,,, M,,, and the other the geostrophic transport 
components M,,, M,,: 


M,=M,,+M,,, (6) 
M,=M,,+M,,- 
The equations with which we must deal are the following: 
—fMye=T | 2=2,3 fle — Tybee (7) 
oP oP 
= =——} =——. 8 
f[M yg ox ? fM zg oy ( ) 
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Let us now form the horizontal divergence of the Ekman wind drift, 
div, M,, from equations (7): 


aie 7g ets Me deemed, (9) 
y 
This is a quantity which has often been contoured for various oceanic 
regions in the belief that it has something to do with upwelling. 
Because of the variation of the Coriolis parameter with latitude, all 
northward or southward geostrophic motions exhibit a horizontal diver- 


gence, div, YW, which can be obtained formally from equations (8): 


div, M,= aan 


In a stationary process the total horizontal divergence, div, M, must 
vanish, hence the sum of equations (9) and (10) vanishes, and we are led 
back to equation (5). The meaning is clearer, however, because we now see 
that what equation (5) expresses is the condition under which the divergence 
of the Ekman wind drift, produced by the wind, is compensated for by the 
divergence of the geostrophic flow. In Ekman’s original picture of the 
stationary wind-driven currents set up by the wind on a plane ocean with 
constant Coriolis parameter, the geostrophic flow is nondivergent, and hence 
he was forced to consider a divergent frictional layer at the bottom to com- 
pensate for the surface Ekman wind-drift convergence. Intense lateral 
friction could produce a divergence, too, but we rule it out by hypothesis. 
It is worth emphasizing that the role assigned to friction in the present 
hypothesis is comparatively insignificant: we suppose that it is dynamically 
important only in the Ekman wind-drift layer. This is a very different 
point of view from that espoused by Neumann, by Hidaka and the Japanese 
school, and by Stockmann and the Russian school. As we have seen, it is 
consistent with the original studies of Sverdrup, Reid, and Munk concerning 
the interior regions of the ocean; and it is possible that friction (in a climato- 
logical-mean sense) is important in certain regions of decay of the intense 
western currents (for example, in the detaching eddies in the Gulf Stream). 
The impossibility of stating definitely whether or not friction plays an 
important dynamical role in the interior of the ocean is one of the reasons 
why the present discussion is not, strictly speaking, a theory, but merely 
a hypothesis. 


(10) 
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THE VORTICITY EQUATION WITH A 
THERMOHALINE PROCESS 


Let us now attempt to introduce a simple thermohaline process into this 
model: For example, we can suppose that, in subtropical latitudes, there is 
a net vertical flux of heat through the surface of the sea tending to increase 
the temperature, and hence to decrease the density of the surface waters, 
and that the density field is actually maintained constant, in spite of this 
net heat flux, by virtue of a slow vertical mass flux pw bringing deep water 
from below into the surface layer. The vertical velocity w is certainly very 
small at the surface and at the bottom. At some intermediate depth z=z, 
it reaches an extreme value (maximum in the subtropics). Since, by mass 
continuity, 


a] 4) 4) 
gpa ge) (11) 


the level of extreme vertical velocity, z=z,;, also corresponds to the level of 
no horizontal divergence. In the layers above and below this level there 
are net divergences, and we shall suppose that they are balanced by 
geostrophic divergences in these layers. We now divide the ocean into two 
such layers, indicating the layer above the level of no horizontal divergence 
by the subscript 1 and the layer below by the subscript 2. We then proceed 
to form the integrated vorticity equations for each layer. First, the follow- 
ing definitions are introduced: 


0 0 
Mis4=| pudz; M,,=| pvdz 
i 7) 


is 3 (12) 
Ma=| pudz; M,.=| pvdz 
D D 


where D is the depth of the bottom of the ocean. We cross-differentiate the 
equations of motion (1), obtaining 


bors] (om) +5, (oo) |= (Gt 2). (13) 
If we now assume that vertical shearing stresses are not important below 
the depth of the Ekman spiral we may take the mid-depth flow as being 
essentially frictionless at z=z;. Furthermore, for simplicity, we shall take 
the depth of the bottom, z=D, as a constant, in order to avoid, at this 
stage, the rather superfluous complications of carrying terms of the form 
(pu) |,-p (2D/x), and so forth. The integrated forms of the vorticity equa- 
tions of the two layers are, then: 


£M,,=—f(pw),+curl 7, (14) 
$M yp=f pw). (15) 
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These equations are of the same form as equation (5), but contain an 
additional term because of the vertical mass flux associated with the ther- 
mohaline circulation. The sum of the two equations is exactly the same as 
equation (5); and hence we see that the introduction of the thermohaline 
circulation in this manner does not produce any net transport over both 
layers in the interior of the ocean, but does result in a kind of internal mode 
of circulation. As we shall see, there is every reason to suppose that the 
flow induced in each layer may be of the same order of magnitude as that 
produced by the wind stress. Since there is so much uncertainty in choosing 
the depth of no motion for a reference level in dynamic computations from 
actual hydrographic data, it is important to form a clear idea of the role of 
this internal thermohaline mode. Indeed, one of the theses of this study is 
that the apparent discrepancy between the theoretical Gulf Stream trans- 
port and that computed ‘dynamically’ (Munk, 1950) is simply a result 
of a confusion of this sort. In a later section this will be considered in 
some detail. 

In order to obtain some idea of the magnitude of the velocity w; required 
to produce a transport equal to that produced by wind stress, we note that 


Munk (1950) gives the value curl 7 = —0-7x 108 g./sec./sec. at 35° N. latitude 
in the North Atlantic Ocean. An upward flux of (ow); =0-85 x 10-4 em./sec. 
(roughly 8-5 cm./day or 30 m./year) produces an equivalent transport. This 
is not an exorbitantly high vertical velocity. It corresponds closely to 
values premised by Riley (1951). It may be compared with what little is 
known of global heat-budget considerations. For example, suppose that in 
subtropical regions there is a net downward heat flux of 150 g.cal./day/em.? 
into the sea surface, and that this heat is used to heat water flowing up from 
below. The average temperature increase is in the neighborhood of 15°C., 
hence the vertical flow is 10 cm./day. These figures should be regarded as 
maximum heat fluxes near the equator: they would require that the ocean 
transport more of the equatorial heat surplus toward the poles than does 
the atmosphere. 

From the point of view of equation (13) the level of no horizontal 
divergence is also the level of no meridional motion. Moreover, if there is a 
small but finite interval of depth in which there is no horizontal divergence, 
this interval will also be a layer of no motion and no vertical shear of geo- 
strophic motion. This may be in part an explanation for Defant’s (1941) 
intuitive choice of the level of no motion in the Atlantic Ocean. In all that 
follows I have been forced to make the assumption that the level of no 
horizontal divergence is of finite thickness, in order to make it a level of no 
motion as well. 

It is perhaps worth while to reémphasize at this point that all we have 
done in this section is to introduce new fields of horizontal divergence 
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associated with a hypothetical thermohaline process, in addition to an 
Ekman wind-drift divergence at the top of the ocean. These additional 
divergences are compensated for in each layer by the divergence of the 
meridional geostrophic flow. 


THE TRANSPORT OF THE WESTERN CURRENTS 


The role of the western currents in the ocean is to maintain continuity of 
mass in the ocean basins bounded by coasts. Generally, the solutions of the 
interior vorticity equations lead to an accumulation of water at certain 
latitudes and in certain layers, and conservation of mass can only be restored 
by rapid western currents in which inertial or viscous forces become of 
dominating importance. In order to compute the transports in the upper, 
G,, and lower, G',, layers of the western currents, we need merely make use 
of the interior solution for the north transport component and the equation 
of mass continuity (11) in its integrated forms: 


4) O 
Be ee My, =(pw);, (16) 
ra) 4) 
On Bt ag ae Pe (17) 


Terms from the differentiation of the limits of integration do not appear, 
since z=z, is a level of no motion. Let us now consider an ocean bounded at 
y=0 by the z-axis and by meridional walls at x=0, r. The boundary con- 
dition at x=r is taken as M,,=M,,=0, following Sverdrup (1947). The 
quantities M,, and M,,. are known from equations (15), and hence we may 
find M,,, M2: 


r 0 
M(x, n=| eer dx, (18) 
r é 


In general we will be led to the result M,,+0 and M,,+0 at «=0. These 
fluxes must be regarded as absorbed by the boundary layer or western 
current solution, hence the western current transports are given by 


Gy(y) = in| Ir pwdydx, (20) 


ee (21) 
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where 


Fe j | *(fpw—eutlt ) de, (22) 


1)=; } “fpwar. (23) 


WIND-DRIVEN AND THERMOHALINE CIRCULATIONS 
IN A RECTANGULAR OCEAN 


In order to portray some of the chief features of the thermohaline circu- 
lation, the analytic solution for a simple rectangular basin has been ob- 
tained. Coasts are placed at x=0, «=r, and y=0, but the ocean is left open 
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Fig. 78. A rectangular ocean with an idealized thermohaline process. The 
vertically integrated transports per unit width in the interior (or central) 
regions of the ocean are shown on the 2, y plane at the right of the figure. The 
broken-line arrows show magnitude and direction of the vertically integrated 
transport per unit width due to wind stress alone. When the idealized vertical 
mass flux is introduced, new patterns of flow appear; the heavy solid arrows 
indicate the transport per unit width of the upper layer, the fine solid arrows 
indicate that of the lower layer. To the left are shown the total transport 
functions of the western currents. 


on the northern edge. A wind stress of the form 7,=—7,cosly is applied 
to the surface, and a vertical flux of pw; = pw, cos ly is fixed at the level of no 
horizontal divergence. The values used in computing the transports ex- 
hibited in fig. 78 were: r=7/l= 6000 km.; 7) =3 dynes/cm.? (this fairly large 
figure for the surface stress is used to bring the curl at y=7/2/ up to Munk’s 
value of 0-7 x 10-8 g./sec./see. and to include the 1-3 meridional factor); 
W,)= 10-4 em./sec.; 8=2 x 10-13/em./sec.; and f=fy. Detail of the western 
currents is not shown, but the total transport of the western current in each 
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layer is shown as a function of latitude by the curves on the left-hand side 
of the diagram. The features of the circulation in this idealized basin, and 
with this very arbitrary choice of distribution of vertical velocity, cannot, 
of course, be applied directly to any real oceanic basin, but they do give an 
indication of the way in which the introduction of the thermohaline process 
modifies the current pattern in the central ocean regions. The broken-line 
arrows indicate direction and amount of transport in the upper layer due to 
wind stress alone, which are similar to those given by the Munk theory. 
They also represent the total transport of both layers in the thermohaline 
model, because, as we have seen above, the transports of the thermohaline 
process are an internal mode of motion, in which the transport of one layer 
is compensated by a counterflow in the other layer. On the left-hand side 
of the figure the transport of the upper layer due to wind alone, or alter- 
natively the transport of both layers (G,+G,) in the thermohaline case, in 
the western current, is indicated by the broken-line curve. Because of the 
choice of parameters, the maximum transport, 36 x 10° m.3/sec., matches 
that computed by Munk (1950) for the Gulf Stream. 

The introduction of vertical mass flux at an arbitrary level of no hori- 
zontal divergence completely alters the picture of the transports, so far as 
individual layers are concerned. In the upper layer (heavy solid arrows and 
curves), this modification is revealed as: (i) a general strengthening of 
transport over the southern half of the basin, (ii) a turning toward the north 
in the northern half, (iii) an increase in the transport in the western current 
in low latitudes, (iv) an abrupt cutoff of the western current somewhat south 
of the maximum westerlies, and (v) a reversal of the western current in the 
northern quarter of the ocean, contrary to the wind-driven current. 

In the lower layer (light solid arrows and curves): (i) strong southwest- 
ward transports occur in the northern half of the basin; (ii) this flow con- 
centrates toward the west into a powerful countercurrent underneath the 
surface western current; (iii) at the latitude of no vertical flux there is no 
meridional transport in the interior of the ocean, and thus the northern half 
of the lower layer is completely cut off from the southern half over the entire 
interior region of the ocean, except for the connection through the counter- 
flow underneath the western current; (iv) the flow in the southern half is 
northeastward and slow; and (v) the water in the lower layer, in the 
southern half, is relatively stagnant, and more than 200 years old (that is, 
since contact with the surface). 

It is interesting to note that the behavior of the thermohaline circulation 
implied by this model in the southern half of the interior of the ocean is 
inverse to the action of a simple Hadley-type convective circulation: 
warmed water moves equatorward, cold water poleward. The northern half 
of the interior solution behaves like a Hadley cell. 

II SGS 
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TRANSPORT CALCULATIONS FROM HYDROGRAPHIC 
DATA FROM THE ATLANTIC OCEAN 


The rectangular ocean basin discussed in the above section is not meant to 
be a true representation of any real ocean basin. It does, however, illustrate 
certain dynamical principles which we will now try to use in a rough, 
tentative way in a discussion of the circulation of the North Atlantic Ocean. 
In order to make the discussion as quantitative as possible, the geostrophic 
transports across various sections in the Atlantic are computed, and an 
elementary method for shifting the level of no motion is introduced. There 
is, of course, no positively certain way to determine the depth of no motion. 
A variety of arguments will therefore be introduced in this connection, none 
of which can be regarded as really satisfactory. 

In the past two years Worthington has embarked upon a program of 
repeating all old (1930-1940) Atlantis hydrographic sections. Moreover, 
he has made his casts to the very bottom, and thus there is now becoming 
available, for the first time, a thorough deep network of sections across the 
western North Atlantic from which an adequate study of the geostrophic 
transport can be made. As Worthington has pointed out in a recent lecture 
at Honolulu, the feature of principal dynamical interest which these new 
sections reveal is that there are horizontal density gradients under the Gulf 
Stream at all depths sufficient to introduce sizable contributions to geo- 
strophic transports as determined by standard ‘dynamical’ computation. 

As a result, the net geostrophic transport computed for each section 
depends strongly upon the choice of the level of no motion, and a change of 
this level—for example, from 2000 m. to the bottom at 5000 m.—makes a 
very large change in the net transport. In other words, the suppositions 
introduced in making the integrations leading to equation (4), namely, that 
the horizontal pressure gradients and motions vanish for all depths beneath 
a certain great depth, are not valid. Nevertheless, in any particular section 
there may be a depth at which (though not below which) horizontal pressure 
gradients parallel to the section and motions normal to the section vanish. 

This uncertainty in the depth of no motion makes it much more difficult 
to compare the results of theory and observation than might be inferred 
from Iselin’s analysis of early observations (1936, 1940). Thus Munk’s (1950) 
comparison of his own theoretically deduced Gulf Stream transport of 
35 x 108 m.3/sec. with the ‘observed’ transport of 74 x 10° m.°%/sec. is fairly 
meaningless. The latter figure was computed by Iselin with an assumed 
depth of no motion at 2000 m., neglecting counterflows above 2000 m. on 
the sides of the main current, and also neglecting counterflows beneath the 
Stream. An even worse discrepancy can be obtained, for example, by using 
Worthington’s new Gulf Stream sections, excluding countercurrents, and 


THERMOHALINE FEATURES 163 


taking the depth of no motion at 5000 m. In this case the net ‘observed’ 
transport is 123 x 10® m.3/sec. (Worthington, 19546). A consistent use of 
the bottom as the depth of no motion has caused Worthington serious 
trouble with the conservation of mass, and he has hoped to circumvent this 
trouble by supposing that there are large-scale nongeostrophic flows in the 
ocean, particularly between Bermuda and the West Indies. I cannot offer 
any positive proof to refute Worthington’s position on this question, but 
it is my opinion! that these difficulties may be mostly resolved in terms of 
the thermohaline circulation, that the depth of no motion in the western 
North Atlantic is, roughly, 1600 m. (in agreement with the estimate by 
Defant, 1941), and that the actual net transport of the Gulf Stream System 
is much more nearly equal to Munk’s theoretical value of 35 x 10® m.3/sec. 
than previous estimates neglecting the deep countercurrent underneath the 
Gulf Stream have indicated. A crucial test of this explanation will be to 
observe by direct means the direction of flow at various depths beneath the 
Gulf Stream. A deep counterflow is a necessary, but not sufficient, feature 
to validate this theory. I therefore emphasize the point that this hypothesis 
of the thermohaline circulation can be struck a mortal blow should direct 
observation disprove the existence of a deep countercurrent on the western 
side of the North Atlantic Ocean.* 

In order to obtain a rough idea of what the transports in the actual Gulf 
Stream might be, the following crude analysis is presented. 


1 T have tried to adhere, as much as possible, and perhaps rather slavishly, to the 
notion that the flow in the deep water is steady and geostrophic, and that the 
observed data which we have are not much distorted by purely local effects or short- 
period fluctuations. Recently other students of the subject have begun to doubt 
both of these assumptions, but the degree to which they are true is not yet known. 
Even granting these assumptions, there are other possible interpretations of the data 
which are quite different from the one I describe here. The conception of the current 
systems presented here is not a unique solution, nor, indeed, are the schemes proposed 
by other writers. I think the multiplicity of possible interpretations needs emphasis 
in order to serve as a warning to the reader. The reason I have used such simple 
assumptions is that otherwise there are so many degrees of freedom that almost any 
interpretation becomes possible. 

2 As the manuscript for this book was completed in mid-1955, it has not been 
possible to include more recent developments in the text; however, in March and 
early April of 1957 a preliminary set of direct measurements of currents under the 
Gulf Stream (using Swallow’s new neutrally buoyant floats) has been made by J. C. 
Swallow and L. V. Worthington and reported in a letter to Nature magazine 179 
(June 8, 1957):1183-1184, and since these measurements have an important bearing 
upon the opinions expressed in the text, I insert an extract from the letter as a footnote 
added in proof: 

‘In choosing the most suitable part of the Gulf Stream system in which to work, 
various factors were considered. The surface velocities in the Stream off the American 
continent usually exceed 200 cm./sec., and it was felt that such strong currents would 

11-2 
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A HYPOTHETICAL ANALYSIS OF THE GEOSTROPHIC 
TRANSPORTS IN THE ATLANTIC OCEAN 


Let us define a function 7'(z, d), the transport per unit depth across a hydro- 
graphic section, as a function of depth, z, and of the reference level (level of 
no motion), z=d. If S(z) is the ratio of the width of the section at depth z 
to that at the surface, taken as completely determined by bathymetry, then 
the transport per unit depth, 7'(z, d), can be expressed in terms of the same 
function referred to the bottom z=6 as the level of no motion in the following 
manner: T ( d, b) 
T(z, d)=T(z, b)— S@) 


The total vertically integrated transport across the section, 7 (d), is a 
function of reference level only: 


TF (d)= | ” 12, d) de. (25) 
b 


S(z). (24) 


make it difficult to keep track of the deep floats. At Stommel’s suggestion, an area 
off Cape Romain, South Carolina, was chosen. Here the shallow (less than 800 m.) 
Florida Current flows over the Blake Plateau, while strong pressure gradients are 
found in the deep water farther offshore. Farther north, towards Cape Hatteras, the 
violent shallow gradients are superimposed on the deep ones. To the south the deep 
gradients dwindle away in a manner not yet understood. On these counts it was 
decided to place the floats as close as was practicable to the junction of the surface 
Stream with the deep water. 

‘The bulk of the work was carried out in about lat. 33° N., and between long. 
75° 30’ W. and 76° W. Excellent Loran coverage exists in this area. Two ships took 
part in this work, which was a joint venture of the National Institute of Oceanography 
and the Woods Hole Oceanographic Institution. The current measurements were made 
by the R.R.S. DISCOVERY II, while the R.V. ATLANTIS occupied hydrographic 
stations, in which serial observations of temperature and salinity were made in order 
to provide a nearly synoptic picture of the deep pressure conditions. 

‘Nine floats were followed, of which seven were in deep south-going water. The 
measurements lasted for periods of 1-4 days, with some overlaps when more than one 
float was being followed. Three floats at 2,500 metres moved in directions between 
south and south-west with mean velocities between 2-6 and 9:5 cm./sec., and four 
floats at 2,800 metres depth moved almost due south with velocities of 9-7—17-4 em./ 
sec. Additional evidence for a south-going deep current was obtained by A. S. 
Laughton, who took underwater photographs of the deflexion of a ball suspended on 
a string, only 50 cm. above the sea floor, in a depth of 3,200 metres. A southward 
movement of about 5 em./sec. was found at that depth. 

‘The ATLANTIS hydrographic stations were made at right angles to the drift of 
each float. The spacing of the stations, eighty-eight in number, did not exceed ten 
miles in the vicinity of the floats. In one case a rectangular pattern of eight stations 
not more than three miles apart was laid. There are some irregularities in the slopes 
of the deep isobars; but the average conditions suggest that the level of no motion in 
this area most probably lies between 1,500 and 2,000 metres, if the southerly move- 
ment of the deep water is to be accounted for by the geostrophic equation.’ 
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In most sections the function S(z) is nearly unity, except in the vicinity 
of the bottom of the section; thus a graph of 7'(z, d) for any particular choice 
of reference level z=d can be approximately transformed to the function 
T(z, da’) for another reference level z=d' by simply shifting the origin of 
function by a constant amount, 


T(z, d') = T(z, d)—T(d’, d), (26) 


so that the level of no motion can be made to move up and down without 
changing the shape of the curve. Thus, in fig. 79, the approximate con- 
‘sequences of a change of reference level on the vertical distribution of 
transport per unit depth can easily be visualized by shifting the origin of 
the abscissa to left or right so that the curve will intersect the ordinate at 
the desired reference level. 

In fig. 79, one of the sections at which we can most certainly assign a 
reference level is section 9 across the Florida Straits. We know from the 
low salinity of the water at the depth of 600-700 m. immediately north of 
the Florida Straits that the reference level must be at the bottom, or below 
it. Direct current measurements and the results of electrical-potential 
recordings between Key West and Havana indicate that the mean flow is 
nearly 26 x 10° m.3/sec.; hence, unless the level of no motion is just about at 
700 m. (the bottom) a serious discrepancy would exist between this estimate 
and the vertically integrated geostrophic transport through the section. 
We therefore take the transport-per-unit-depth curve for section 9, shown 
in fig. 79, as a starting point in our analysis of the geostrophic transport of 
the Gulf Stream in particular and the Atlantic Ocean in general. If the 
flow in the western North Atlantic is geostrophic, then the total mass trans- 
port through sections 8 and 9 must be very nearly the same as that through 
section 6 or section 7. So far as I can determine, this requirement is best 
met by the choice of reference levels indicated in fig. 79. In order to demon- 
strate the degree to which the total mass conservation is achieved by this 
choice, I have prepared the more detailed graph shown here as fig. 80, using 
some of Worthington’s newest data. 

It is reassuring to note that not only do the total transports 7 (2) balance 
approximately, but also the transports at each depth 7'(z, d) balance, very 
roughly, when the depth of no motion is chosen in the interval of 1200- 
1600 m. Evidently we cannot expect precise agreement. The shifting of the 
level of no motion uniformly over the whole section may be too crude; there 
may be differences in the level of no motion at various points across each 
section. In any case, the choice of the interval of 1200-1600 m. leads to 
much smaller violations of mass continuity than the choice of the bottom 
as reference level. 

We can extend the argument of total mass conservation to other 
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sections, but, because of vertical flow, we cannot expect the individual 
transport per unit depth to be conserved over large areas. Thus, for 
example, we may reasonably expect the vertically integrated transport 
to the north across section 7 to be equal to that to the south across 
section 3, because the Atlantic Ocean is essentially closed toward the 


north. 
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Fig. 79. Geostrophic transport per unit depth, 7(z,d), across various 
sections of the Atlantic Ocean. The total transport across various latitude 
circles due to Ekman wind drift alone is indicated at the lower right corner; 
the depth distribution is estimated. Western boundary currents are the major 
features in sections 6, 7, 8, 9, 10, 11, 12, and 15. The other sections represent 
interior (or central) oceanic regions. 

The broad arrow-point symbol indicates the depth to which actual data 
extend, below which the curve may be extrapolated; the long, flat rectangle 
indicates the depth of the solid bottom. The stations used in the respective 
sections were as follows (A= Atlantis; C= Caryn; Ch=Challenger, 1932 cruises; 
D= Discovery; and M= Meteor): 


Sec. Sta. Sec. Sta. Sec. Sta. 
1 Chl, 12 8 A 5313, 5328 14 M 188, 201 
2 Ch 1287, 1293 9° A 51155, 5162 15 M 165, 167 
3 A 4628, 4637 LON) FC2202; 218 16 M 168, 180 
4 A 5307; MVE-3 ie Avb2o2, 0204, 17 « M Isis sas 
5 A 5203, 5210 5235, 5262 184 
6 A 5180, 5202 12 M 289, 293 18 M 75, 87 
7 A 5295, 5312 13. ~M 218, 289 19 D 1156, 1165 
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Sections 3, 4, and 5 are horse-latitude sections where winds are small but 
the curl of the wind stress is at a maximum. They include parts of the 
interior that are well removed from the western boundary currents. This is 
an important region of the ocean from a theoretical point of view. We should 
establish here whether the interior solution really does correspond to what 
is expected from the theory of wind-driven currents. It is not an easy 
section to construct with confidence: no single section is available across the 
interior of the North Atlantic with many deep-water data (the Dana depth 
determinations are unreliable). Moreover, 
two stations such as those used in section Syms Pies eo ae 
3 are really not sufficient below the Mid- 
Atlantic Ridge: the pressure gradients in 
the two basins may be different. Section 4 
is composite. Below 2000 m. it is most 
unsatisfactory. Section 5 is made from 
Worthington’s reliable new deep-water 
data, but it does not include the entire 
width of the North Atlantic. It is com- 
puted from the bottom. The conservation 
of mass requires that the reference level 
be placed at about 1200m., and it is 
interesting to note that this depth roughly 
coincides with the level at which the great 
salt pall from the Mediterranean seems to 
hang. From the fact that this highly saline 
water extends almost due westward in 
a tongue issuing from the Straits of 
Gibraltar it may be shown that at depths : 
of between 1000 and 2000 m. the value Fig. 80. Detail of the trans- 
of the meridional component of T(z, b) Pe eer eecee eee 
must be less than 0-2 x 10®m.3/100m./sec. own in fig. 79. 

Were it greater, the wedge of Mediterra- 


10 12 


DEPTH IN METERS 


nean Water would not point due west—there would be a perceptible tilt. 
A choice of about 1600 m. is indicated by the curves in fig. 79. In section 5, 
for which good deep-water data are available, this does not imply large 
transports of deep water and bottom water. The total vertically inte- 
grated transport across section 3 appears to be of the order of 


30-40 x 10° m.3/sec., 


a figure not inconsistent with the interior (central) solution given by the 
Munk theory. 
Let us now proceed to discuss other sections in fig. 79. The northernmost 
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section, section 1, goes from shore to shore, and hence should include interior 
and boundary currents. The 1932 Challenger stations used extended only 
to a depth of 2000 m. If 2000 m. is used as the depth of no motion there is a 
net flux above 2000 m. of roughly 30x 10° m.3/sec. to the north across 
section 1. By shifting the level of no motion to about 1500 m., or even 
shallower, it is possible to obtain a balance of mass flux across the section 
(as depicted in fig. 79). The flux due to Ekman wind drift across section 1 is 
shown, approximately, at the bottom of fig. 79. It is too small to play any 
important role in the net balance across section 1. Our reasoning forces us 
to admit a deep southward flow across section 1 between 1500 m. and the 
bottom. A similar pattern of flow appears to occur across section 2, but, 
although the data for this section extend to 3000 m., the section does not 
extend all the way to the western coast. 

In terms of the foregoing hypotheses, the deep currents, which at sections 
1 and 2 are distributed over much of the interior, are compressed into a 
narrow western boundary current at sections 6 and 7 and do not exist in the 
interior at sections 3, 4, and 5. Because of the shallow sill in the Florida 
Straits (indeed, even over the Blake Plateau), the deep countercurrent 
shown below 1600 m. on sections 6 and 7 splits away (sections 8 and 10) 
from that part of the surface current (sections 9 and 11) which passes 
through the Caribbean. Because the curl of the wind stress is nearly uniform 
between 20 and 40°N., the surface currents across sections 9 and 1] are 
more properly regarded as wind-driven than are the surface currents across 
sections 8 and 10. By hypothesis, the schematic model of the rectangular 
ocean described above causes this deep current to be absorbed by vertical 
motion into the surface layer of the interior of the subtropical North 
Atlantic. In the real ocean it appears that only part of the deep flow across 
sections 6, 7, 8, and 10 can be so absorbed, and that most of it must 
flow across the equator into the South Atlantic Ocean, where it appears on 
both section 12, beneath the Guiana Current, and section 15, beneath the 
Brazil Current. The deep currents of the South Atlantic are mostly confined 
to narrow western boundary currents. Sections 12, 13, and 14 are within 
10° of the equator, where there is some reason to be slightly suspicious of the 
geostrophic calculations. Also, the flow of the Ekman wind drift cannot be 
ignored at these latitudes. It is interesting to note that the flow of Antarctic 
Intermediate Water appears to be stronger across section 13 than across 
section 12; one wonders whether this indicates vertical velocity near 1000 m. 
or a surprisingly deep influence of curl of the wind stress. The transport 
across interior section 14 seems remarkably low, especially when one notes 
by what a narrow margin the geostrophic transport misses being canceled 
by the Ekman wind drift. 

We now come to sections 15 and 16, which are the counterpart in the 
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South Atlantic of sections 7 and 3 in the North Atlantic.* Sections 3 and 16 
are both situated in subtropical regions of high atmospheric pressure: the 
regions of maximum curl of the mean wind stress between the mid-latitude 
westerlies and the trades. Sections 3 and 16 show about the same equator- 
ward transport; that across section 16 is, perhaps, a little deeper. Each is 
very much what would be expected from the wind-stress theory alone. It is 
only when we compare the two western current sections 7 and 15 that we 
see a striking difference between hemispheres. Although section 15 is very 
complicated, the depth of the level of no motion may be ascertained fairly 
well from water-mass considerations. The choice shown in fig. 79 indicates 
three such levels between different water masses. It is consistent with the 
descriptions by previous investigators. The difference in appearance be- 
tween section 7 and section 15 is so remarkable as to call for further 
description and elaboration. In both sections the total vertically integrated 
transport poleward must be equal to the total wind drift equatorward 
across sections 3 and 16 respectively, from mass-conservation arguments. 
The wind-driven component of section 7 and that of section 15 must be 
essentially the same. The marked difference between the function for the 
Gulf Stream and that of the Brazil Current must therefore be due to the 
different way in which the thermohaline western boundary currents enter. 
In the Gulf Stream the thermohaline circulation reénforces the poleward 
surface motion; in the Brazil Current the situation is quite the contrary. 
The separation of the wind and thermohaline contributions is shown 
schematically in fig. 81. 

In summary, we see, therefore, that the wind-driven component of the 
circulation is much the same in the two hemispheres, but that in the western 
boundary currents the thermohaline circulation plays quite different roles 
in the Gulf Stream and in the Brazil Current, making the former appear 
much more strongly on charts of the topography of surfaces, equal pressure, 
or density. 

For purposes of comparison, the function 7'(z, d) is drawn also for the 
transport through a north-south Discovery section across the Antarctic 
Circumpolar Current, section 19, 4500 m. being taken as the reference level. 
The trifling Benguela Current is shown in section 17. 

Transport curves for a rough two-layer division of the circulation of 
the Atlantic Ocean are indicated in fig. 82,a@ and 6. Each line represents 
about 10 x 10® m.3/sec. 


3 Since this chapter was written, a much more precise analysis of deep water 
transports in the South Atlantic Ocean has been published by Georg Wiist: in Papers 
in Marine Biology and Oceanography, Supplement to Vol. 3 of Deep-Sea Research, 
1955, entitled ‘Stromgeschwindigkeiten im Tiefen- und Bodenwasser des Atlantischen 
Ozeans auf Grund dynamischer Berechnung der Meteor-Profile der Deutschen 
Atlantischen Expedition 1925/27’, pp. 373-397. 
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WEAKNESSES AND DIFFICULTIES IN THE 
PRESENT THERMOHALINE HYPOTHESIS 


The chief weakness of the treatment outlined here, as I see it, is the hypo- 
thetical nature of the quantity pw;. How shall we measure it directly? If 
frictional stresses in deep water are really negligible, then much information 
concerning the field of w might be obtained from an analysis of central 
oceanic hydrographic data, since 
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INTERMEDIATE 
Fig. 81. Schematic diagram of the possible explanation of the very different 
transport-per-unit-depth curves for the Gulf Stream and the Brazil Current. 


but the available very deep data in the North Atlantic are quite insufficient 
for such a study. 

Bottom topography must be of considerable importance in the thermo- 
haline circulation, in which horizontal pressure gradients do not vanish at 
the bottom. Major topographic features like the Mid-Atlantic Ridge must 
be of great importance; along the eastern rise of this ridge we might very 
well find a deep southward current (of the boundary-layer type) similar to 
the one under the Gulf Stream. Where bottom water is confined to deep 
basins the thermohaline circulation may not be of direct importance, and 
the circulation within the basin may actually be driven by small frictional 
stresses exerted upon it by the thermohaline circulation above sill depth, 
in much the same manner as in an ordinary wind-driven circulation. 

Finally, the vertical distribution of vertical velocity must be of great 
importance in determining the vertical temperature and salinity structure 
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in the ocean. Were the eddy diffusion coefficients independent of depth, the 
depth of maximum w might nearly correspond to the bottom of the main 
thermocline. 

In summary, then, the line of argument and interpretation of data pre- 
sented here, although by no means conclusive, do suggest: (i) that the 
interior of the subtropical and equatorial Atlantic Ocean is essentially wind- 
driven, from a dynamical point of view; (ii) that in subpolar regions thermo- 
haline processes are important; (iii) that in the western boundary currents, 
effects due to the mass-conservation requirements imposed by both the 
wind stress and the thermohaline process are important, even in subtropical 
and equatorial latitudes; (iv) that the portion of the surface-layer western 
current which is driven by wind stress is intimately connected with Munk’s 
gyres, or a cell-like structure of wind currents; (v) that the thermohaline 
portion of the surface-layer western current flows from 40°S. latitude to 
40°N. latitude almost unchanged in transport; (vi) that the deep circu- 
lation at all latitudes is essentially thermohaline, but in subtropical and 
equatorial latitudes is confined to a narrow western current of the (perhaps 
inertial) boundary-layer type; and (vii) that, therefore, the mean vertical 
flux of mass in the interior of subtropical regions at depths of 1000-2000 m. 
is probably considerably less than that used in the theoretical rectangular 
ocean (less than 10 em./day).* 


4 Since the preparation of the manuscript for this book, I have written and 
published a review article which supplements the ideas presented in this chapter: 
‘A Survey of Ocean Current Theory’, Deep-Sea Research, 4 (1957):149-184. Also, 
I have attempted to make an analysis of hydrographic data in the North Atlantic 
according to equation (27) under the title: ‘On the Determination of the Depth of 
No Meridional Motion’, Deep-Sea Research, 3 (1956): 273-278. 


Chapter Twelve 
CONCLUDING REMARKS 


The central North Atlantic Ocean is covered to a depth of about 700 m. 
with a layer of warm water, slowly drifting toward the southwest under the 
combined influence of the stress of the wind and the earth’s rotation. This 
central body of warm water, called the Sargasso Sea, does not reach the 
western coasts of the Atlantic Ocean (the eastern coast of North America), 
but is bounded some distance offshore by the Gulf Stream—a narrow, 
intense, northeastward-flowing current which returns to the north again 
the southward-driven Sargasso Sea water that has passed through the 
Caribbean and has turned through the Florida Straits. 

The Gulf Stream flows along the western boundary of the warm Sargasso 
Sea surface water. As the Stream turns toward the east, off the Grand 
Banks, it acts as a kind of dynamic barrier, or dam, which, by virtue of 
Coriolis forces, restrains the warm Sargasso water from overflowing the 

‘colder northern water of the North Atlantic. The water in the Stream is 
not significantly different in temperature from the large mass of warm 
water which lies to the right of its direction of motion. The Gulf Stream is 
not an ocean river of hot water. The intensity of flow of the Stream, the 
Stream’s direction, and its temperature are not primary climatic factors in 
determining the climate of Europe; but the role which it plays in deter- 
mining the northern boundaries and average temperature structure of the 
Sargasso Sea must be of critical climatic importance. However, we do not 
yet know how this role is physically connected with such observables as 
total transport. 

Certain remarks concerning the very nature of oceanographic research 
itself are made in a later section of this final chapter. 
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JUSTIFICATION FOR TREATING THE GULF 
STREAM SYSTEM AS AN ENTITY 


Perhaps it has occurred to the reader to wonder whether there is any 
justification for treating the Gulf Stream System as a physical system in 
any sense separable from the general circulation of the North Atlantic 
Ocean. Ifthe Gulf Stream were a river of very warm water flowing through 
a colder environment, there would be little doubt that it should be con- 
sidered a distinct physical phenomenon. As we have seen in the preceding 
chapters, the most clear-cut feature of the entire system is that it is really 
a flow along the very edge of the juncture of a mass of cold water and a 
mass of warm water. Therefore, one is led to inquire whether the primary 
physical phenomena to discuss are not the origin and nature of the two 
contrasting water masses, and whether it is not proper to regard the Gulf 
Stream current system as merely a secondary feature associated with their 
zone of contact. 

From a historical point of view, the Gulf Stream has always been the 
primary feature of interest in the western North Atlantic, and the efforts 
of the first scientific oceanographic explorations were concentrated on it 
because it was so clear-cut an entity. ‘It’ could be crossed in a day; ‘its’ 
total water transport could be determined; ‘its’ changes in position could 
be plotted; ‘it’ had an edge, a core of high velocity, a countercurrent, and 
so on. Compared to other regions of the ocean, it seemed one where 
techniques of measurement, such as the indirect velocity determinations 
based on the geostrophic approximation, could be made with a reasonable 
degree of precision. When one considers the vast central oceanic areas, or 
the eastern sides of oceans, where the current features are very faint and 
indistinct—almost beyond the reach of standard means of measurement— 
it is no wonder that oceanographers on the western sides of oceans often 
consider themselves especially fortunate to possess such strong current 
systems near their coasts and, as a result, have spent much more time 
cruising in these strong currents than in obtaining detailed information 
about the great water masses on either side. 

Also, from the point of view of theoretical studies, the Gulf Stream 
System has a distinct reality of its own. The outstanding feature of these 
theories (discussed in Chapters VII and VIII) is that the transport solu- 
tions are independent of any assumptions regarding the density structure. 
The Gulf Stream would exist, according to these theories, whether it is 
driven by wind stress or by a thermohaline process. It would exist whether 
the ocean were half as deep as it is, or ten times as deep. The velocities 
would be different, but the pattern of the transport lines would be the 
same. The Gulf Stream System emerges from these theories as a boundary 
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phenomenon occurring at a certain boundary of the oceans where high- 
order derivatives in the governing equation (elsewhere negligible) are 
locally important. In this sense, the Gulf Stream ‘exists’ as a physical 
entity in the same sense as the boundary layer next to an aerofoil exists: 
it is a recognizable part of a larger physical system. 

The idea of a western boundary current—whether it be viscous or 
inertial—is of basic importance because it simplifies the physics and 
mathematics of the central regions of the ocean. The higher-order terms in 
the vorticity equations are restricted to a narrow coastal region, and the 
remainder of the ocean can be studied with very simple equations (as we 
have already seen in Chapter XI). This possibility of emancipating the 
interior solution from boundary conditions at the west may very well be 
the most important lesson to be learned from the study of the Gulf Stream. 


CLIMATE AND THE GULF STREAM 


There is scarcely any more firmly rooted idea in the mind of the layman 
than the notion that the Gulf Stream keeps the European climate warm. 
So long as it was believed that the Gulf Stream was a river of warm water, 
this idea did make sense. It is no longer possible to be so certain of the 
direct climatological influence of the Gulf Stream, for it now seems that it 
is not so much the Stream itself that is important, as the position and 
temperature of the large mass of warm water on its right-hand flank. 
Any argument which purports to relate climate to the transport of the 
Gulf Stream must include some account of the physical relations between 
the Stream and the density (and thermal) structure of the water on either 
side. In fact, Iselin (1940) went so far as to speculate that warming of the 
European climate might actually be least during periods of increasing 
transport of the Gulf Stream. This idea was based on the hypothesis that 
the processes which produce the warm surface masses of the Central 
Atlantic Water are more or less constant in time, and hence, owing to the 
geostrophic relation, an increasing transport of the North Atlantic Gyre 
must be accompanied by simultaneous deepening of the thermocline in the 
Sargasso Sea and a radial shrinking of the current system. In this way, 
warm surface water would be withdrawn from the north, and the European 
climate might be colder. Conversely, a weakening of the transport of the 
current system would, Iselin reasoned, be accompanied by a rising thermo- 
cline throughout the Sargasso Sea, and the excess warm water would force 
the current system radially outward and farther northward, and might 
even succeed in discharging quantities of surface water to high latitudes, 
thus warming the European climate. There is no convincing evidence to 
prove that this (or any other) sequence of events actually takes place, nor 
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can we even discuss the hypotheses very satisfactorily; in particular, we 
know very little about the rates of formation of water masses. But these 
arguments are of special interest, because they confound the popular idea 
of the role of the Gulf Stream in European climate. For all we know, the 
European climate might actually be warmer if the direction of rotation of 
the North Atlantic Eddy were reversed. 

During the past thirty years there has been an increase of about 2° in 
the surface temperatures of the Norwegian Sea, and also a decrease of 
perhaps 50 m. in the depth of the 10° C. isotherm throughout the Sargasso 
Sea. (The statistical significance of this change of position of the isotherm 
is doubtful.) These two isolated bits of information suggest that the North 
Atlantic surface circulation is slowing down and that the Gulf Stream 
transport is consequently decreasing, and that this, in turn, results in 
greater warming of the coast of Europe. This indication favors Iselin’s 
proposal. 

The theory of wind-driven ocean currents suggests that the effect of 
diminished winds would make the North Atlantic circulation shallower in 
the middle, but would not change its radius. Toward the periphery of the 
circulation the isotherms might become deeper. However, it is important 
to remember that the theory, in its present state of development, deals only 
with integrated velocities, and hence that there is no clear theoretical 
indication of the ways in which temperature, heat transport, and surface- 
current velocity might vary with the wind. Moreover, the climate and 
winds over the ocean are not unaffected by the state of ocean currents. 
As a result, the theory of wind-driven currents permits us to draw in- 
ferences about the change of currents resulting from only very slight 
changes of the winds. It seems to me quite reckless to assert, for example, 
that some current of the Carboniferous Glacial Period, flowing in an ocean 
which does not even exist now, had a greater or lesser transport than the 
present-day Gulf Stream. The more these hypothetical oceans of the past 
depart from the configuration and climate of the oceans of today, the less 
we can hope to explain and describe their behavior.! When the disparity 
between supposed conditions of the past and observed conditions of the 
present becomes as great as that envisaged in Chamberlain’s theory of 
the reversal of the deep circulation due to the sinking of highly saline 
water at the equator, the unhappy physical oceanographer can make no 
constructive comment. 


1 Many examples of purely verbal theories that invoke physical processes with 
which our present primitive theoretical models cannot yet cope are given in Chapter IIT 
of C. E. P. Brooks’s book Climate Through the Ages (1949). Another example, 
published since this text was written, is ‘A Theory of Ice Ages’, by M. Ewing and 
W.L. Donn, Science, 123 (1956): 1061-1066. 
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SOME REMARKS OF A POLEMICAL NATURE 


By far the greatest part of physical oceanographic knowledge was accumu- 
lated by the following process: 


First, a grand cruise, or expedition, brings back data obtained at many 
hydrographic stations. 

Secondly, extensive plots, graphs, and tabulations of the data are made 
and published for the benefit of future generations (e.g., in Bulletin 
Hydrographique). 

Thirdly, certain of the more striking features of the data plots are noted. 

Fourthly, some plausible hypotheses are advanced to explain them. 


At this stage the procedure had usually exhausted the energy of those in- 
volved, and almost always the funds, and the study usually stopped. Ina 
few projects, such as the Discovery Antarctic studies, the Scripps Marine 
Life Program, and the Woods Hole Gulf Stream studies, attempts were 
made repeatedly (third step) to bring the striking features into better focus, 
the first and second steps being repeated in various combinations and per- 
mutations, with the goal of presenting the main features of the findings in 
great detail. Thus, much of the effort in the Gulf Stream work has been 
devoted to attempts to detect changes in the total transport of the Stream, 
to surveys of the shapes of meanders, and to determinations of the 
‘width’ of the Stream and of the transverse surface velocity profile. The 
plausible hypotheses (fourth step) are modified as needed, to avoid conflict 
with the observations. Sometimes the choice of hypothesis is easy and the 
modification slight, because a variety of different plausible hypotheses all 
fit the observations. 

These first four steps are absolutely necessary, of course; without them 
we would know nothing at all about the ocean as it really is. But for the 
full development of the science there must be one additional step: 


Fifthly, the plausible hypotheses must be tested by specially designed 
observations; in this way theories can be rejected or accepted, or may be 
modified to become acceptable. 


As an example of an untested hypothesis one may take any particular 
choice of the ‘depth of no motion’ for dynamic computations of the Gulf 
Stream (see Chapter XI). A knowledge of the truth or falsity of this 
hypothesis will become indispensable once the theoretical models begin to 
approach reality. To make a satisfactory test will exercise our ingenuity. 
The surface currents are too strong to permit the making of accurate direct 
current measurements from an anchored ship. Oceanographers may have 
to devise and employ some new observational technique in order to make 
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this important test, but the great obstacles should not blind us to the 
necessity of undertaking it. Perhaps the most promising technique for 
direct deep-current measurements underneath the Gulf Stream is a method 
which is now being refined by Dr J. C. Swallow at the British National 
Institute of Oceanography and was recently described at a London con- 
ference on the deep ocean circulation. A small float, so ballasted as to 
reach and remain at a fixed predetermined depth, carries a battery- 
powered ultrasonic noisemaker. A ship with a directional hydrophone 
array can track the drift of the float, and hence obtain direct measures of 
the drift of deep currents. Whether this technique can be refined to 
measure currents as accurately (say to 0-5 cm./sec.) as will be needed to 
determine the level of no motion below the Gulf Stream, remains to be 
seen; but it is important that this test should be made.? 

Another example of an untested hypothesis is the entire conception that 
large-scale lateral mixing occurs in the ocean. As yet, there does not seem 
to be any obvious way to test this idea. 

I should like to make it clear, finally, that I am not belittling the survey 
type of oceanography, nor even purely theoretical speculation. I am 
pleading that more attention be given to a difficult middle ground: the 
testing of hypotheses. I have not explored this middle ground very 
thoroughly, and the few examples given in this book may not even be the 
important ones; but perhaps they are illustrative of the point of view in 
which attention is directed not toward a purely descriptive art, nor toward 
analytical refinements of idealized oceans, but toward an understanding of 
the physical processes which control the hydrodynamics of oceanic circu- 
lation. Too much of the theory of oceanography has depended upon purely 
hypothetical physical processes. Many of the hypotheses suggested have 
a peculiar dreamlike quality, and it behooves us to submit them to especial 
scrutiny and to test them by observation. 


2 For a report of very recent work of this kind see Chapter XI, footnote 2. 


The definitions of only the chief mathematical and physical symbols used 
in this book are indicated below. For a subscript or mark appended to a 
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symbol and used only once or twice, no special listing is given. 


SYMBOLS FROM THE ROMAN ALPHABET 


average coefficient of horizontal-eddy viscosity 
north-south dimension of a rectangular basin 
internal long-wave velocity, ./(9’ho) 


undisturbed depth of a homogeneous ocean, or of the top 
layer of a two-layer ocean 


body force representing wind stress at surface; also, in 
Chapter VIII, a function of y 


Coriolis parameter 


a function of y in Chapter VIII; total transport of 
western current in Chapter XI 


acceleration due to gravity 
reduced gravity, equal to g(Ap/p) 


intensity of the earth’s magnetic field; subscripts denote 
components 


perturbed depth 
horizontal coefficient of eddy viscosity 


vertical coefficient of eddy viscosity 
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wave number 

length of cable between towed electrodes 

horizontal components of vertically integrated mass 
transport per unit width 

pressure 

an artifice, coefficient of friction on the bottom of a 
homogeneous ocean, assuming a linear dependence on 
velocity 

radius of curvature of a streamline 

east-west dimension of a rectangular basin 

north-south dimension of wind system 

transport per unit depth 

vertically integrated transport 


the T-S diagram, with temperature as ordinate, salinity 
as abscissa; used in water-mass analysis 


mean velocity, used to linearize inertial terms in the per- 
turbation theory 


horizontal component of velocity in the x-direction 
horizontal component of velocity in the y-direction 


horizontal codrdinate, usually meaning (in this book) 
directed eastward 


horizontal codrdinate, usually meaning (in this book) 
directed northward 


vertical codrdinate, directed upward 


SYMBOLS FROM THE GREEK ALPHABET 


1/A 

meridional gradient of the Coriolis parameter 
drag coefficient of wind on sea surface 
relative vorticity, vertical component 

radius of deformation 

wavelength unit (microns) of radiation 
meander frequency 

density of water 

density of air 

10° (p—1) 

a special measure of density, defined on p. 31 
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shearing stress across a horizontal surface 


stream function, or mass-transport function, depending 
on the context 


SPECIAL SYMBOLS 


parts per thousand, unit of salinity 
star, used to denote quantities in inertial boundary layer 


arrow, used over a quantity to indicate a vector 
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HYDROGRAPHIC STATION DATA 


All hydrographic station data obtained in the Gulf Stream, consisting of 
precise temperature and salinity measurements, have been published at 
Copenhagen since 1910 by the Conseil permanent international pour 
l’exploration de la mer in its annual Bulletin hydrographique, for the period 
beginning with the year 1908. Because of the great care taken in editing 
the Bulletin, there is an interval of about five years between the date of 
observation and the date of publication. 

A complete index, including charts, of all the hydrographic data from 
the Atlantic Ocean prior to 1932 is contained in International Aspects of 
Oceanography, by T. Wayland Vaughan (1937). 

Two papers by C. O’D. Iselin (1936, 1940) may be used as an index for 
Atlantis stations made in the Gulf Stream between 1930 and 1940; the 
original data are to be found in Bulletin hydrographique. 

Since the Second World War there have been relatively few hydro- 
graphic sections of the Gulf Stream; specifically: 


1946, Sept. Atlantis stations 44144423 (no salinities), a Chesa- 
peake Bay—Bermuda section. 

1947, June. Altlantis stations 4565-4570, a section off Cape 
Hatteras. 

1947, Dec. Atlantis stations 4623-4627, a section along 65° W. 
longitude. 

1949, July-Sept. Atlantis stations 4760-4832, an areal survey 
near the Grand Banks. 
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1950, June. International Ice Patrol stations 4175-4184, a 
section near 50° W. longitude (see Soule, 1950). 

1950, Oct._Nov. Atlantis stations 4842-4882, three sections at 
68° W. longitude (see Worthington, 19545). 

1953, Aug. Atlantis stations 5081-5125, three sections at about 
73, 70, and 61° W. longitude respectively. 

1954, Sept. Atlantis stations 5176-5202, along 65° W. longitude, 
a very deep section. 

1955, June. Allantis stations 5295-5312, Chesapeake Bay to 
Bermuda, very deep. 

1955, June-July. Atlantis stations 5344-5360, east-southeast of 
Cape Romain, S.C., very deep. 

1956, Nov. Atlantis stations 5416-5444, south of Grand Banks on 
50° W. longitude, very deep. 


The data from all these hydrographic stations have been submitted to 
the Bulletin hydrographique, and have been or will eventually be published. 


BATHYTHERMOGRAPH DATA 


The data from bathythermograph measurements are so voluminous that 
they have never been published in full. The information is available only 
in raw, unreduced form in the card files of the Woods Hole Oceanographic 
Institution, and in the form of rough running sections drawn aboard ship. 
A few sections have been published (see Iselin and Fuglister, 1948). 
8. L. Strack has made an analysis of fifty-eight bathythermograph crossings 
of the Stream to study the relation of temperature at the surface to that 
at various depths as a function of season, but this study has not been 
published. 


DISCUSSIONS AND SUMMARIES OF ATLANTIC 
OCEAN DATA 


In addition to the sources mentioned above there are several general 
studies, atlases, chart collections, and so on, containing much useful in- 
formation about the Atlantic Ocean, but not concerned specifically with the 
Gulf Stream. 

1. Charts of surface temperature, salinity, and density, according to 
season, are given in Giinther Bohnecke’s (1936) atlas, consisting of seventy- 
four plates and entitled Atlas zw: Temperatur, Salzgehalt und Dichte an der 
Oberfliiche des Atlantischen Ozeans, in Deutsche Ailantische Expedition auf 
dem Forschungs- und Vermessungsschiff Meteor, 1925-1927, Wissenschaft- 
liche Ergebnisse. (Band 5—Atlas.) 
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2. Charts of distribution of mean temperature, salinity, and density at 
various depths are included in an atlas of 103 plates prepared by Georg 
Wiist and Albert Defant (1936), Atlas zur Schichtung und Zirkulation des 
Atlantischen Ozeans, which is also in the multivolume work on the 
scientific results of the Meteor Expedition (Band 6—Atlas). 

3. For charts of the distribution of properties on various o, surfaces 
in the southern North Atlantic Ocean, see R. B. Montgomery (1938a). 

4. For a study of the distribution of various chemical constituents of 
the Atlantic, including an attempt to explain the distribution quantita- 
tively, see G. A. Riley (1951). 

5. Charts of mean winds, charts of climatic factors, and similar source 
material are to be found in Atlas of Climatic Charts of the Oceans, prepared 
by the United States Weather Bureau (1938). 

6. Charts showing surface winds, surface currents, and features of 
navigational interest are available in Atlas of Pilot Charts, Atlantic Ocean, 
prepared by the United States Navy, Hydrographic Office (1950). 

7. For a general geography of the Atlantic Ocean, in a popular vein 
and yet informative and fairly comprehensive, see Gerhard Schott’s 
Geographie des Atlantischen Ozeans, published at Hamburg in three 
editions: 1912, 1926, and 1944. It is the third edition that has been 
utilized in the present study. 
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TABLE A 


SHorT TABLE OF S1IGMA-T VALUES 


Values of co: 


Temperature, 
in degrees Salinity, in °/,. 
Centigrade 
34 35 36 37 
0 27-32 28°13 28-93 = 
3 27°11 27-91 28-70 —— 
6 26:78 27°57 28-36 —— 
9 26:36 27°14 27-92 sos 
Ie 25°83 26-61 27°38 28-16 
15 25-22 25-99 26°76 27-53 
18 24-53 25-29 26°05 26°82 
21 23°75 24-51 25°27 26:04 
24 22-91 23-66 24-42 25:18 
vA7f 21-99 22-74. 23°49 24-25 
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